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ABSTRACT 

This paper discusses the electrical and magnetic design and analysis of a permanent 
magnet generation module for naval applications. Numerous design issues are addressed and 
several issues are raised about the potential improvements a PM generation system can offer. A 
proposed 16 MW PM generation module design is presented along with a detailed design 
methodology. 

Eighty different machines and power conversion modules are sized, designed, and 
analyzed with a final design selected. Specifically, sizing and detailed machine design and 
analysis is performed examining the effects of numerous parameters including number of phases, 
number of poles, magnetic geometry, machine dimensions, and material types. Analytical 
models are developed to study rotor losses caused by stator winding time and space harmonics 
and slot space harmonics. Power electronics and conversion modules to connect the high-speed 
generator to a DC distribution system are designed and analyzed. In depth simulation of the 
eighty complete systems is performed using the software programs MATLAB (Version 12.0, 
Mathworks) and PSIM (Version 6.0, Powersim, Inc.). 

The 16 MW permanent magnet generation module, consisting of the generator and 
associated power electronics, provides an excellent alternative to traditional wound rotor 
synchronous machines. The final design offers significant reductions in both weight and 
volume. Specifically, it is estimated that the PM generation module has a 7x reduction in 
volume and a lOx reduction in weight compared to similarly rated wound rotor systems. These 
reductions can provide flexibility to naval architects since power, weight, and volume are 
integral parts of the design and construction processes. However, further study is necessary to 
verify the PM generation modules thermal, structural, and mechanical performance. 
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Chapter 1 Introduction 


1.1 Purpose 

The purpose of this thesis is to design and analyze a permanent magnet generator and 
power module for naval applications. When deciding whether to implement an electrical 
technology or component onto a naval vessel, the size, weight, and cost are the major factors for 
successful integration. Significant performance improvements and cost reduction of power 
electronics, coupled with the availability and decreasing cost of high energy permanent magnet 
(PM) materials makes PM generators attractive for naval usage. These machines offer numerous 
desirable features, including light weight, small size, simple mechanical construction, easy 
maintenance, good reliability, and high efficiency [1]. 

Before analysis of a generator can begin, it must be properly designed for typical naval 
power requirements. This involves sizing the generator along with designing the associated 
power electronics for connecting the machine to the distribution system. A specific concern 
associated with PM generators is possible inefficiencies and excessive heating; in particular rotor 
losses caused by space and time harmonics during the energy conversion processes. The 
optimum machine design is one that delivers the required power through a matching process 
between the generator and the power electronic converter [2]. 

1.2 Problem 

The Navy’s commitment to develop an integrated electric power system for the next 
generation warships offers the expectation of using the installed generation capacity to power 
ship propulsion, advanced weapons components, and high power combat control systems [3]. As 
these electrical loads increase, it becomes increasingly important to efficiently utilize installed 
power as well as develop smaller, effective power generation systems. Navy ships are extremely 
high performance systems and therefore power and weight considerations are integral parts of 
the design process. 

The life cycle of a navy ship is on average 2-3 times longer than a commercial ship and 
therefore navy ships undergo excessive modernization and upgrades throughout their service life. 
Many of the newer components have significantly higher power requirements than the originals 


13 





putting a burden on the power generation system. Therefore, the Navy is moving toward 
designing all electric ships with integrated power systems (IPS) and increased power generation 
that can be efficiently managed to meet future demands. 

The effective integration of electrical power in future naval ships requires the 
development of technologies that ensure volume and mass reduction in critical mechanisms. 
Military ships require higher power density components, impose more stringent acoustic and EM 
signature requirements, and subject systems to harsher environments than commercial 
applications [4]. Rotating generators, coupled with prime movers, need to be lighter in weight 
and higher in power density. High-speed PM generators provide a substantial reduction in size 
and weight making them a logical choice for naval applications. Currently, the Navy has not 
designed or built a high power (megawatt) PM generator and therefore the need exists. 

In conjunction with constructing a PM generator, a DC bus architecture is one of the 
preferred schemes for the future [5]. DC power distribution systems can offer a size and weight 
reduction compared to high-power AC systems [3]. With a DC bus distribution, the PM 
generator can be optimized independent of producing 60 Hz frequency as required in the past. 
However, conversion of the high-frequency generator AC output to DC requires power 
electronics. Rectification of AC to DC presents the problem of creating harmonics in the input 
current which are then reflected back onto the generator causing rotor losses. In addition, the 
generator produces space harmonics which also produce losses in the rotor. Therefore, the PM 
generator and power electronics module (PEM) need to de designed and optimized to deliver 
constant DC power while minimizing machine losses. 

1.3 Background 

1.3.1 History 

The Navy has designed and built electric ships since the early part of the 20* century. 

The original advantages perceived for electric ships, superior performance, reduced manning, 
arrangement flexibility, and fuel efficiency, are still relevant today [6]. In the early part of the 
20* century, diesel-electric submarines, small surface ships, and some battleships and carriers 
had electric propulsion. By the late 1940s, mechanical drive systems became popular because of 
improvements in metallurgy and manufacturing. However, the capability of mechanical 
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transmissions reached its limits in the late 20 century spawning a renewed interest in electric 
drive and integrated power systems for military applications. 

For the last thirty years, the commercial industry has designed and operated integrated 
power systems (IPS). An IPS is a ship architectural paradigm in which the ship’s power and 
propulsion are provided hy a common electrical distribution system instead of having a separate 
mechanical drive for propulsion. For military applications, an IPS provides numerous benefits: 

• Decreased life cycle costs because of increased fuel economy and efficiency (a Navy ship 
with IPS may consume 10-25% less fuel than a similar ship with mechanical drive [7]) 

• Increased ship design and arrangement flexibility since the ship is not limited to having a 
long mechanical shaft line 

• Reduced system complexity 

• Higher degree of modular design using power components 

• Broad industrial base for implementing IPS design 

• High levels of automation and control 

• Increased power available for non-propulsion uses since a Navy ship spends a large 
portion of its time operating at low propulsion levels (approximately 95% of the time) 

• Increased stealth, survivability, and payload 

Since an IPS provides power for both the ship loads and propulsion, larger generation capability 
is required. As the Navy builds its new ships with IPS architectures, compact, high-power 
generation systems must be examined to help facilitate implementation of the new designs. 

1.3.2 Power Generation & Distribution 

Almost all naval core power generators to date are air-cooled 50/60 Hz machines that are 
military derivatives of commercial generators and are therefore size and weight excessive [8]. A 
typical turbine generator system is shown in Figure 1 and Table 1 contains nominal 
characteristics for several current commercial and naval generators. 
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Figure 1: Typical Turbine Generator System 


Table 1: Examples of Current Generator Characteristics 
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1 

3 

21 

3600 

4.7 

4 

3.4 

50.0 

2 

3 

36.5 

3600 

6.2 

4.37 

3.76 

63.6 

3 

3 

26 

3600 

6.2 

3.56 

4.5 

68.1 

4 

3 

25 

3600 

5.18 

3.1 

4.15 

57.3 


With the advent of high-power, cost-effective power electronics, it is no longer necessary to 
generate power at 50/60 Hz so generators can be optimized independent of frequency. High¬ 
speed, power-dense generators become the logical choice for naval purposes. Permanent magnet 
machines are ideal for this high-speed application due to their simple structure and high power 
density [9]. 

Since the generator can now be designed to produce higher frequencies, distribution 
architectures are not limited to being 60 Hz. Either a high frequency AC system or a DC system 
can be designed through the use of power electronics modules, with the DC distribution being 
preferred because of its advantages in size and weight. In this thesis, a high-power DC zonal 
architecture is assumed using solid state converters to generate AC where needed. Each zone is 
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electronically isolated from the other zones with automatic fault detection and reconfiguration to 
provide continuous power during damaged conditions. The PM generator sets and power 
electronic conversion modules serve as the backbone thereby providing a reliable power system 
for navy ships. 

1.4 Scope 

The scope of this thesis is limited to the PM generator and associated power electronics 
AC-DC conversion module. The following is accomplished: 

• Determine the electrical power requirements for a Navy IPS ship in order to properly size 
the generator 

• Compare typical wound rotor machine design to a permanent magnet design to determine 
applicability for IPS applications 

• Conduct material analysis and selection for the generator design 

• Perform initial PM generator detailed design 

• Design the power electronics conversion module to perform high-power AC-DC 
conversion 

• Conduct detailed analysis of rotor losses of the PM generator, in particular those caused 
by time and space harmonics 

• Perform numerous iterations of machine and power electronics designs to develop 
optimized generation scheme 


17 


Page Intentionally Left Blank 


18 







Chapter 2 Power Requirements and Machine Selection 

2.1 Machine & Module Requirements 

On board Navy ships, electricity is used to provide power to virtually all components, 
including mission systems, support systems, combat systems, and communications systems. In 
addition, as ships continue to be upgraded and modernized, more power is needed for newer 
combat systems and weapons components. Most current naval platforms have some form of 
mechanical propulsion system with separate ship service electrical generators supplying the 
ship’s power. With an IPS ship, the ship’s generators provide power for propulsion and the 
ship’s service loads, and through proper utilization, power is efficiently managed. 

To properly size the PM generator, the power requirements must be identified and 
therefore a typical load list is developed for an EPS naval ship and is included in Appendix A. 
The overall power requirement for the generator is 16 MW. Since size and weight are important 
factors and the generator can be optimized independent of frequency, high-speed operation and 
maximum power-density are desired. Therefore, the highest possible speed is selected while 
ensuring the PM generator is compatible with both gas turbines and steam turbines. 

Traditionally, gas turbines run at much higher speeds than steam turbines causing the 
steam turbines to be more limiting. From information collected from the Elliot Turbomachinery 
Company, Inc., 16,000 RPM is approximately the highest speed steam turbine that can 
reasonably be constmcted at the megawatt power level [13]. Therefore, to provide a degree of 
conservatism, 13,000 RPM is selected for the nominal design speed for the PM generator. 

The power electronics module (PEM) converts the AC voltage from the generator to 
700 VDC. Overall, the PEM and generator must be designed so that losses suffered by the 
permanent magnets on the generator rotor are minimal. Table 2 lists the general requirements for 
the entire system. 

2.2 Machine Selection 

Military ships require high power density components and improved acoustic and 
electromagnetic signature requirements while subjecting systems to harsh environments [14]. It 
is therefore important to ensure the power generation system is capable and efficient. 
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Table 2: General Machine & Module Requirements 


1 .., . ’ : VIS; • 

Of 

Generator Power 

16 MW 

Generator Speed 

13,000 RPM 

PEM Output Voltage 

700±5 VDC 

PEM Output Ripple 

0.7 VDC (@ 16 MW) 

Generator Rotor Losses 

Minimal 


2.2.1 Permanent Magnet versus Wound Rotor 

Reducing the size and weight of ship’s turbine generator sets offers significant 
advantages to naval architects. Replacing older generators with lightweight ones could make it 
possible to decrease the size of some generator sets by as much as 50% [15]. PM generators 
therefore become an attractive alternative compared to wound rotor machines because of the 
availability and decreasing cost of high energy PM materials along with improved power 
electronics. 

A wound rotor generator normally consists of armature windings on a stationary stator 
frame with field windings on an inner rotor. The rotor is turned by a prime mover, usually a gas 
or steam turbine, and current is supplied to the field windings through brushes or a brushless 
exciter. As the current-carrying field windings rotate past the stator windings, current is 
produced in the stator windings through Faraday’s Law. An example of a wound rotor machine 
is shown in Figure 2 [8] and Figure 3 [16]. 



Figure 2: Example of Wound Rotor Generator 





















Figure 3: Cross Section of Wound Rotor Generator 


Wound rotor generators have been the backbone of power generation for the U. S. Navy because 
they are a proven technology that is battle tested. They offer steady voltage regulation using the 
field windings, a large air gap for producing the rotor magnetic flux, low fault currents, and high 
power capabilities. 

However, the machines tend to be complex, weight excessive, and they require field 
windings which limit design alternatives. There are several drivers which cause these problems 
to occur. First, to generate the necessary magnetic flux levels, wound rotor generators have large 
pole pitches to support the required field windings. These pole pitch windings in turn require 
larger end turns and thick back iron to support the magnetic flux, both of which contribute to 
increased size. Second, because of the winding losses in the rotor, large cooling systems can be 
required thus increasing the number of support components. 

High speed generators offer a reduction in machine size and weight because as a 
machine’s speed increases, its size decreases for a given output power. The PM generator is 
ideal for high-speed applications because of its simple structure and high power density [9]. In a 
PM generator, the rotor field windings are replaced by permanent magnets which do not require 
additional excitation. As the permanent magnets are rotated by the prime mover, current is 
produced in the stator windings. An example of a PM generator is shown in Figure 4 [8]. 
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Figure 4: Example of PM Generator 

PM generators offer several advantages: they have no rotor windings so they are less 
complicated; they have high efficiencies; the gap field flux is not dependent on large pole pitches 
so the machine requires less back iron and can have a greater number of smaller poles; and they 
usually require smaller and fewer support systems. Assuming the same flux density and 
circumferential arc, doubling the number of poles produces the same radial flux but requires half 
the stator core thickness, as shown in Figure 5. 



Figure 5: Flux vs. Number of Poles 
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However, PM generators also have some disadvantages. They do not possess field 
excitation control and therefore voltage regulation can be problematic. This can be corrected by 
using external voltage control such as large capacitor banks or power electronics, as well as 
choosing the turns on the stator winding properly to produce the anticipated required nominal 
voltage. Additionally, since the permanent magnet fields cannot be turned off, there exists the 
risk of excessive currents in the event of an internal fault. This problem can also be solved 
through the design of the turbine governor and controller or dynamic braking. Overall, the 
advantages of the PM generator over the traditional wound rotor generator make it a better 
alternative for high-speed navy applications. A summary comparing the different designs is 
given in Table 3. 


Table 3: Comparison of Wound Rotor and PM Generators 





Wound Rotor 

Steady voltage regulation with 
field windings 

Weight excessive 

High power capabilities 

Large size 

Large air gap for flux 

Rotor windings & associated 
losses 

Low fault currents 

Large support systems 

Proven, robust design 


Permanent Magnet 

Less complicated 

Lack of inherent voltage 
regulation 

Reduced size and weight 

Potential fault currents 

High efficiency 

Magnet losses 

No excitation supply or field 
windings 


High speed applicability 



2.2.2 Type of Permanent Magnet Machine 

There are numerous layout possibilities for permanent magnet machines and only the 
most common are discussed here. These include radial flux inner rotor, radial flux outer rotor, 
and axial flux designs. In most PM machines, flux crosses from the rotor to the stator in the 
radial direction [17]. The first t 5 q)e, the radial flux inner rotor design, is the closest configuration 
to the classical AC synchronous generator. An example of this design is shown in Figure 6 [17]. 
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Figure 6: Example of Inner Rotor PM Machine 

In this type of machine, the windings are placed on the stator, either in slots or in a 
slotless ring, and the magnets are surface mounted on the rotor or buried in the rotor. Buried 
magnet designs often result in rotors that are larger than equivalent surface-magnet machines 
with high-energy magnets [18]. Buried magnet machines can also have significant structural 
issues in high-power applications [19]. When the magnets are surface mounted and the machine 
is operated at high speed, the magnets are often secured with a retaining device made of either 
alloy steel or carbon-fiber. Overall, the inner rotor machine possesses high torque/power 
capability and good heat conduction and cooling properties making it ideal for high-speed, 
higher-power applications [18]. 

The radial flux outer rotor machines are commonly used in hard disk drives, small 
computer ventilation fans, and some blowers. This type of design is very efficient, low-cost, 
easy to manufacture, and applicable for low-power applications [18]. It is the opposite of the 
inner rotor because the stator in on the inside with the rotor and magnets on the outside. A cross 
section of an outer rotor machine is shown in Figure 7 [17]. 
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Figure 7: Example of Outer Rotor PM Machine 


The axial flux machine is significantly different than the previous two because flux flows 
in the axial direction vice radial direction and the windings are oriented radially vice axially (see 
Figure 8 for an example diagram [20]). 



The main advantages of this design are their low cost, flat shape, and smooth rotation. However, 
if axial-flux machines are operated at high speeds (above 1000 RPM), eddy-current losses and 
heating can become excessive [18]. Also, stator construction is difficult because it must be 
laminated circumferentially. An example of this design is the turntable for a record player. 
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Overall, because of its inherent advantages in heat removal and cooling, the abundance of 
manufacturing capabilities, and its high-power, high-speed applicability, the radial flux inner 
rotor with surface mounted magnets is selected for the 16 MW PM generator design. 
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Chapter 3 Material Selection and Machine Initial Design 


3.1 Material Selection 

One of the key considerations during the electromagnetic, structural, and thermal design of 
a permanent-magnet machine is the selection of the magnet, stator, and rotor materials [21]. 
Machine output, heat rise, weight, and cost are a few of the characteristics which are directly 
influenced hy selection of the machine materials [22]. 

3.1.1 Permanent Magnets 

The size and performance of high-speed PM generators depend on the permanent magnet 
material properties [9]. The magnets must be selected to provide the necessary air gap magnetic 
field and ample coercive force to compensate for possible damaging effects while minimizing 
the volume of material because of cost and weight considerations [23]. 

Ferromagnetic materials are the most common substances used in the construction of 
machines and their properties are normally described using B-H curves and hysteresis loops. 
These curves represent an average material characteristic that reflects the non-linear property of 
the permeability of the material but ignores the multi-valued properties [17]. An example of a B- 
H curve is shown in Figure 9. 



Figure 9: Example of B-H Curve 
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Several basic magnetic properties are of critical importance for the permanent magnets in a PM 
machine: 

• Remnant Flux Density (Br): It is the value of the flux density remaining after 
magnetization and it directly influences the air gap flux and magnet sizes. 

• Coercivity (He): It is the value of magnetizing field needed to reduce the flux density in 
the magnet to zero and it gives a first order estimate of a magnet’s resistance to 
demagnetization. 

• Energy Product (BHmax): It is the maximum energy product of the magnet and it is 
inversely proportional to the total magnet volume required. 

• Recoil Permeability (prec): It is the gradient of the B-H curve and it gives the magnet’s 
ability to return to its initial magnetization after subjected to damaging forces. If the 
magnet goes below Hk, then it will recoil along a lower line resulting in a lower magnet 
flux density. 

• Load Line: It is a line drawn from the origin to the magnet operating point on the 
hysteresis curve (Bn,). The magnitude of the slope of the load line is the permeance 
coefficient. 

Permanent magnet materials come in many varieties and the four most common types for 
machine applications are Alnico, Ferrites, SmCo material, and NdFeB material. Table 4 and 
Figure 10 show the characteristics and typical B-H curves for these materials [18]. 


Table 4: Magnet Material Properties 








Coercivity (He) 

kA/m 

40- 130 

180-400 

800-1500 

800- 1900 

Recoil Permeability 

(Prec) 


1.9-7 

1.05-1.15 



Energy Product (BHmax) 

klW 

20-100 

24-36 

140-220 

180-320 

Maximum Temperature 


500-550 

250 

250-350 

100 - 200 

Br Temperature 
Coefficient 

%/’c 

-0.01 to -0.02 

-0.2 

-0.05 

-0.08 to -0.15 


The rare-earth magnets, SmCo and NdFeB, have become more popular for high performance 
applications because of their greater power density, high coercivity, high flux densities, and 
linearity of the demagnetization curves [24]. 
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Alnico — Ferrite ■” SmCo .NdFeB 



Figure 10: Typical Magnet B-H Curves 

Between the two rare-earth permanent magnets, NdFeB is preferred because it is cheaper and 
more readily available. It does possess some adverse characteristics such as moderate corrosion 
and lower resistance to temperature effects, but these can be controlled using surface treatments 
and proper cooling [25]. Therefore, NdFeB magnets are selected for use in the PM generator 
with the conservatively assumed values listed in Table 5. 


Table 5: Selected Magnet Properties 





Remanence (Br) 


1.2 

Coercivity (He) 

kA/m 

900 

Recoil Permeability (prec) 


1.05 

Energy Product (BHmax) 


260 

Maximum Temperature 


180 

Resistivity 

pfl/m 

1.43 
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3.1.2 Stator and Rotor Material 

The type of material chosen for the stator and rotor is important because it impacts the 
machine losses and efficiency. The rotor is usually built from the same material as the stator for 
ease of construction but is can be made of any economical steel provided it is strong enough for 
the given function [18]. No one material is optimum for every application and the normal 
criteria for selection are cost, permeability, core losses, and saturation flux. It is important that 
the material act as a flux guide and absorb the minimum amount of magnetomotive force (MMF) 
so that the flux is concentrated in the air gap. In addition, the material should minimize core 
losses including hysteresis and eddy current losses. 

High-quality, non-oriented, electrical grade lamination steels are typically used in most 
machines because the laminations help minimize losses. The four main materials are low carbon 
steels, silicon (Si) steels, nickel (Ni) alloy steels, and cobalt (Co) alloy steels. Low carbon steels 
are the lowest cost and are used in high volume applications where high core losses are 
acceptable. Silicon steels usually have 3% silicon which increases the resistivity to reduce eddy 
current losses. They are selected and specified based on core loss, with each grade (Ml9, M27, 
M36, and M43) having higher core losses and lower cost [22]. The lamination thickness is a 
tradeoff between cost and performance and the most common sizes are 0.014 in, 0.0185 in, and 
0.025 in (29 gauge, 26 gauge, and 24 gauge). 

Nickel alloys are either 49% or 80% nickel and they have lower losses than the silicon 
steel but are much more expensive. In addition, they require careful handling and not suited for 
high flux density environments (above 0.8 T) because of saturation. The cobalt alloys are only 
used in extremely high-performance situations such as military aircraft and space applications 
because of the high cost. Table 6 summarizes the different stator materials and the Ml9, 29- 
gauge electrical silicon steel is selected for the PM generator because it is economical, its thin 
laminations minimize losses, and it has a saturation flux density of about 1.8 T [2], [18], [22]. 
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Table 6: Laminated Steel Properties 








Low Carbon 
Steel 

Fair 

Good 

Good 

Best 

0.5 

Si Steel 

Good 

Good 

Fair 

Good 

1.0 

Thin Si Steel 

Better 

Good 

Fair 

Fair 

10.0 


Good 

Fair 



12.0 

mmmm 

Better 

Low 

Best 


15.0 


Good 

Best 

Good 


45.0 


3.2 Machine Design Parameters 
3.2.1 Stator Mechanical Design 

The stator is an important part of the machine because it serves as the main structural 
component, it provides the housing for the armature windings, and it completes the flux path for 
the magnetic circuit. The main consideration in the mechanical design of the stator is whether to 
make it slotted or slotless. A slotless stator has the armature windings located in the air gap of 
the machine as shown in Figure 11 [19]. 



Figure 11: Slotless Stator Design 

One of the advantages of the slotless construction is unique winding layouts are possible to meet 
specific performance goals. Another advantage is that the space available for the armature 
windings increases by a factor of about two since there are no stator teeth. This produces lower 
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conductor losses since less current flows in each winding. The flux density is reduced, however, 
because the effective air gap is much larger since the windings are in the air gap. Overall, there 
exists a higher electrical loading and a lower magnetic loading. 

One disadvantage of the slotless design is there are no good conduction paths to remove 
the heat generated from the windings. This reduces the allowable current density in the windings 
and lowers the power output. Another disadvantage is that the windings are directly exposed to 
the rotating flux which raises the possibility of additional eddy-current loss in the conductors and 
further losses due to circulating currents in the windings [18]. Overall, the performance of a 
slotless stator is almost always lower than that of an equivalent slotted stator design and 
therefore slotless stators do not appear often in high-power applications [17]. 

Slotted stators are the traditional stator design and consist of openings around the stator 
for the armature windings as shown in Figure 12 [26]. The openings provide rigid housings for 
the conductors and associated insulation. 



Figure 12: Slotted Stator Design 

Stator slots vary in size and shape with the most common configurations being rectangular or 
trapezoidal. In this paper, the slots are assumed to be approximately rectangular as shown in 
Figure 13 and contain form-wound windings so that the depression width is the same as the slot 
top width. 
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Figure 13: Stator Slot Geometry 

Slotting is used because is provides a good way to achieve a narrow air gap length while keeping 
the winding conductors close to the magnets to maximize the flux linkage. The slots also greatly 
increase the surface contact area between the windings and stator steel providing a path of low 
thermal resistance for good heat conduction which is important for keeping the windings and 
magnets cool. 

The resulting narrow air gap from the slots makes the permeance greater and therefore 
the air gap flux density greater producing a more powerful machine. In addition, the depression 
in the slot tops help control parasitic losses in the rotor by improving the uniformity of the air 
gap field. The limits of the size of the slots are twofold: the magnetic teeth must be able to carry 
the air-gap flux without saturating and the slots must be large enough to support the necessary 
current density in the windings. Typical limits for stator current density are shown in Table 7 
and in this paper it is assumed that the limit on current density (J) is 3000 A/cm [18]. 


Table 7: Stator Current Densities 


1 '";, '-‘/liiMiMMbod 


Natural Convection 


Fan Cooled 

800-1200 

Liquid Cooled 

2300 - 3200 


The disadvantages of the slots are that cogging torque may be a problem and it can be 
costly to insert the windings if proper construction techniques are not used. Overall, however, 
slotted designs are preferred in high-power applications and therefore a slotted stator is selected 
for the 16 MW PM generator. 
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The number of slots in the machine is usually a function of the number of phases and 
windings and can vary based on the application. The initial design of the generator assumes a 
three-phase machine but the number of phases will be examined to determine an optimum 
design. In order to allow for high-power operation and the possibility of a high number of 
phases, 36 slots is chosen for the initial generator design. 


3.2.2 Rotor Mechanical Design 

For high-speed applications, the rotor aspect ratio, defined as length-to-diameter (L/D), is 
a critical parameter. If it is relatively low, then the rotor has high stiffness and good dynamics 
but a large diameter which increases the weight and makes magnet retention extremely difficult. 
Additionally, the centrifugal force on the surface-mounted magnets is directly proportional to the 
rotor diameter so the rotor radial size must not be excessive. 

Permanent magnet machines offer flexibility in selecting pole sizes which allows for 
smaller diameters. They are therefore ideal for high-speed applications because they can have 
higher L/D ratios. This is because they do not have rotor field windings which have end turns 
necessitating big pole pitches and large diameters. A normal L/D ratio for a wound rotor 
machine is 0.5 - 1.0 compared to 1 - 3 for a PM machine [27]. Staying close to these ranges 
usually provides a first order estimate of satisfactory machine dynamic performance and 
acceptably low vibrations or oscillations. 

The rotor radius and the rotational speed also determine the tip speed of the machine 
which is the surface velocity of the rotor (as defined by Eqn 3-1). 

''tip = 

where ©m = angular speed (rad/sec) 

R = rotor radius (m) 


Eqn 3-1 

For most rotating machines, the upper limit on tip speed is between 100 - 250 m/s depending on 
the design. For surface magnet PM machines, retaining sleeves are sometimes used to help keep 
the magnets in place and allow for higher speeds. These sleeves can be constructed from alloy 
steel, carbon fiber, or other materials. The metal sleeves usually provide increased mechanical 
performance but have eddy current losses. 
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The carbon fiber and graphite composite sleeves have high strength-to-weight ratios 
which produce a thin sleeve and the sleeve’s lower conductivities yield reduced eddy current 
losses. However, the carbon fiber and graphite composite sleeves have lower temperature 
ratings and lower thermal conductivities making heat removal and increased cooling for the 
magnets and sleeve important issues [28]. Overall, the use of a retaining sleeve is necessary for 
the 16 MW generator since it is operating at high-speed and this allows the maximum tip speed 
limit at the rotor surface to be 200 m/s. The actual material for the retaining sleeve is examined 
later when detailed rotor loss analysis is performed. 


3.2.3 Number of Poles and Magnet Pole Design 

The optimum number of poles is a complex function depending on a number of factors 
including the magnet material, the speed of rotation, the desired output frequency, and the 
mechanical assembly of the rotor. An even number of poles is always used because this provides 
a balanced rotational design. As the number of poles increases, the individual pole pitch goes 
down which reduces the amount of stator back iron needed to support the magnetic flux. In 
addition, for a given power/torque, as the pole number rises, the required magnet volume 
decreases. 

Assuming a constant mechanical rotation speed, the generated electrical frequency is 
proportional to the number of poles as shown in Eqn 3-2. 

N-(2p) = 120-f 

where N = speed (RPM) 

p = number of pole pairs 
f = electrical frequency (Hz) 


Eqn 3-2 

If a PM generator is going to be the source for a DC bus through a rectifier system, a high pole 
number is desirable because as the electrical frequency increases, support components such as 
filter capacitors and inductors can be much smaller. Therefore, for a given rotational speed, one 
cheap and efficient solution is to have a higher number of pole pairs and frequency [27]. 
However, as the frequency increases, higher stator losses result because core losses are 
proportional to frequency squared. In addition, as the pole number gets larger, the number of 
slots per pole per phase decreases and can cause the voltage waveforms to become less 
sinusoidal so all factors must be considered. 
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The pole arc of the magnets can also be varied. Magnets seldom span the full pole pitch 
because the flux at the transition between north and south poles leaks between poles without 
linking the coils in the stator. The gaps between the poles usually contain non-magnet pieces, 
such as soft-iron, so that no flux crosses over the air gap between magnets. A full pole arc is 0me 
= 180°E and produces a full voltage waveform but has increased harmonic content. As the pole 
arc is reduced (up to 20 - 30 %) and those areas are filled in with soft-iron pieces, the resulting 
flux waveform is more sinusoidal and has fewer harmonics and therefore lower rotor losses [29]. 

The magnet poles are sometimes skewed to reduce cogging torque and smooth out 
variations in air gap reluctance, flux, and voltage waveforms. Skewing of the magnets occurs 
axially along the length of the rotor to provide a constant rotational torque and prevent pole 
pieces from exactly lining up with stator teeth. A skew factor is used to account for this effect 
and is shown in Eqn 3-3. 


hn = 



where ©g = skew angle, radE 
n = harmonic number 


Eqn 3-3 

As the pole number is increased, the stator conductors-per-pole decreases so that the per- 
unit inductance and synchronous reactance decreases with higher pole number. This can 
sometimes result in improved performance of the machine since the reactance is lower. Overall, 
the initial 16MW generator has 6 poles but this is examined later to determine an optimal design. 

3.2.4 Magnetic Dimensions 

The primary magnetic dimensions that affect a PM machine are the air gap and the 
magnet height. These two parameters play a major role in determining the air gap magnetic 
field, the air gap flux density, and the induced voltage in the machine. To a first order 
approximation, the air-gap flux density (Bg) can be represented by Eqn 3-4 [30]. 

The radial air gap is usually made as small as possible to maximize the air gap flux 
density, minimize the flux leakage, and produce a lower reluctance value since the air gap 
constitutes the largest part of the machine permeance/reluctance. However, the use of rare-earth 
permanent magnets (NdFeB or SmCo) with their higher flux density and coercive force permit 
some flexibility in the size of the air gap. 
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®g= 


‘m 


+ g 


-B, 


m 

where = magnet height (mm) 
g = air gap (mm) 

= magnet remnant flux density (T) 

■| 

Eqn 3-4 

Once the permanent magnet material is selected, the desired air gap flux density and 
induced voltage help determine the magnet height needed. If the magnet height is too large, the 
air gap flux density might be significant enough to cause the stator core material to saturate 
which reduces machine performance. The goal is to use the minimal amount of magnet material 
to achieve the desired effect because this reduces the size and weight of the machine and 
decreases the magnet material cost. Also, losses in the magnets can be reduced by using smaller 
magnets. In order to provide uniform magnetic fields, the magnet height is usually larger than 
the air gap by a factor of 5 - 10. 

3.2.5 Number of Phases 

In general, the number of phases affects a machine’s power, current, and voltage ratings 
as shown in Eqn 3-5. If the power is fixed, then as the number of phases increases, the phase 
voltage and/or current decreases, assuming the total number of turns is constant. 

|P+jQl =qVI 

where P = real power (W) 

Q = reactive power (VAR) 
q = number of phases 
V = RMS phase voltage (V) 

I = RMS current (A) 

Eqn 3-5 

Most motors and generators are three-phase machines because it is the industry standard, it is the 
most common form of power, and it is the lowest number of phases that produces balaneed 
torque with out pulsations in rotating machines. However, higher utilizations in generators can 
be achieved with higher phase numbers especially if the generator is connected through power 
electronics to a DC bus distribution. This is because the higher number of phases produces 
lower ripple in the DC bus voltage. 
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However, the AC line current harmonics are more substantial in increased phase 
machines because the triple-n harmonics are higher order as the phase number increases. For 
example, a 3-phase machine suppresses harmonics of order 3n, a 5-phase machine eliminates 
order 5n, and a 7-phase machine removes order 7n. Therefore, in higher phase machines, a 
greater number of large harmonics result in the AC line current. Also, as the number of phases 
increases, the phase inductances and reactances change since there are a greater number of 
windings influencing each other. 

Most machines are usually designed with the phases balanced meaning that they have an 
evenly-spaced phase distribution around the stator of the machine. This produces voltage 
waveforms that are identical in shape from phase to phase but differ by a phase offset angle. In 
order to initially size the PM generator, it is assumed to have three phases but this will be 
optimized later in conjunction with the power electronics module. It is also assumed that the 
phases are always balanced. 

3.2.6 Slots per Pole per Phase 

The number of slots per pole per phase (m) is an extremely important design parameter 
when considering generator design and it is calculated using Eqn 3-6. It is used to help 
determine the relationship and interactions between the rotor poles and the stator windings as 
well as shape the generated back voltage of the machine. When m is an integer, the machine is 
an integral slot machine and when m has a fractional part, it is a fractional slot machine. 



2pq 

where Ng = number of slots 

p = pole pairs 
q = number of phases 

Eqn 3-6 

In an integral slot machine, the back EMFs of all of the coils making up a phase winding 
are in phase with each other and add up so that the final voltage amplitude is the direct sum of 
the individual coil voltages. In a fractional slot machine, the back EMF of all of the coils are not 
in phase so the net voltage has a different shape than the individual winding voltages. Varying 
the number of slots/pole/phase is one method used to produce a more sinusoidal voltage 
waveform and reduce the harmonics generated by the machine. 
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3.2.7 Stator Windings 


The stator windings are the location where the generator voltage is induced due to the 
time varying magnetic flux caused by the permanent magnets on the rotor. In a slotted machine, 
the winding arrangement is used to help shape the back voltage to produce a more sinusoidal 
waveform. The windings can be distributed by three methods: pitch, skew, or 
breadth/distribution. 

The pitch of a winding (a) refers to the angular displacement between the sides of a coil, 
usually expressed in electrical degrees or radians. When the individual coil pitch differs from 
180^ E, the winding is said to be short-pitched or fractional-pitched. This causes angular 
segments where the back voltage is zero because the flux linkage is constant and can help 
produce a sinusoidal waveform when multiple coils are connected. It also has the advantage of 
lowering the coil resistance and making the stator end windings more manageable. 

Windings in the stator can also be skewed axially along the length of the machine. This 
requires the stator slots to be more intricately designed which complicates the mechanical 
construction of large machines. Therefore, since the generator being designed is a large, high- 
power machine, skewing of the stator windings is not used but skewing of the rotor is employed. 

The breadth of a stator winding results from the coils occupying a distribution or range of 
slots within a phase belt. A stator winding normally consists of several coils each separated by 
an electrical angle y. The distribution of the coils causes each to link the rotor flux slightly out of 
phase with each other so when they are added together, they produce a more sinusoidal 
waveform. 

Within each stator slot, there are geometric size constraints which determine how many 
conductors can be placed in a slot. In smaller machines, coils are composed of round insulated 
wires that are placed in the stator slot along with insulation material. A slot fill factor (X®) is used 
to determine how much of the slot cross-sectional area is occupied by winding material as shown 
in Eqn 3-7. 

^ _ WindingArea 
^ TotalSlotArea 


Eqn 3-7 

In larger machines, form-wound windings are used for ease of construction and for better 
performance. A sketch of what a form-wound winding looks like is shown in Figure 14. 
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Typically, machines contain two coil sides per slot making the winding a double-layer design 
[17], [18], [19]. Overall, slot fill factors vary in value from 0.30 - 0.70, depending on the 
number and size of the conductors in the slots as well as the amount of labor utilized. In this 
paper, a slot fill factor of 0.50 is assumed. 








Figure 14: Example of Form-Wound Winding 

In conductors that carry high-frequency currents, skin effect can become an issue and 
affect the operation of the machine. Skin effect is caused by eddy currents in the windings 
themselves due to the changing magnetic field. These eddy currents force the current flowing in 
the conductor to crowd to the outer edges of the conductor. This in turn causes the current to 
flow through a smaller cross-sectional area and increase the resistance of the conductor. 
However, the generator under design is expected to operate at less than 2 kHz and for 
frequencies below 12 kHz, Rac/Rdc < 1-01 so skin effect can be neglected [32]. 

Within a phase, stator windings can be connected in wye or delta patterns as well as 
series or parallel. Almost all machines use series, wye-connected windings because they provide 
the safest alternative. This is because in a delta or parallel connection, the back EMFs can 
produce circulating currents which can result in addition losses, heating, or damage. Therefore, 
wye series connected windings are selected for use in the designs in this paper. 
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3.3 Machine Calculated Parameters 
3.3.1 Basic Model 

Since the machine is assumed to balanced, parameters can be determined on a per-phase 
basis and then applied to all of the phases. Each phase of the machine can therefore be modeled 
as shown in Figure 15. 



Figure 15: Per Phase Model 

The armature resistance (Ra) is the resistance of the windings of the machine and it is usually 
relatively small. The synchronous inductance (Ls) of the machine comes from the inductance of 
the windings and is composed of the air gap inductance, the slot leakage inductance, and the end- 
turn inductance. The back voltage (Ea) is produced through the flux linkage in the windings 
from the rotating magnetic field in the machine. Lastly, Va is the terminal voltage and is found 
using basic circuit analysis once the other parameters are known. 

3.3.2 Winding Resistances 

The stator coils in the machine are made of copper and therefore have some resistance to 
the current flow. This resistance of the copper phase windings is calculated using Eqn 3-8. 


where I = length of conductor 
a = winding conductivity 
A = winding cross-sectional area 

Eqn 3-8 

The length of the conductor comes from the windings traveling twice the length of the machine 
and twice around the end turns of the machine. It is assumed that the end turns follow roughly a 
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circular path from one side of the machine to the other where the radius of the circle is the 
distance to one half the stator slot height. The cross-sectional area of the conductor is obtained 
from the slot area and slot fill factor as shown in Eqn 3-9, assuming form-wound windings. 

A-^ 

2-Nc 

where Ag = slot area 

Ng = turns per coil 


Eqn 3-9 

3.3.3 Winding & Magnet Factors 

As discussed in section 3.2.7, windings are normally not full-pitched or concentrated but 
rather are short-pitched and have breadth associated with them. To account for these effects, a 
winding factor (kw) is utilized which is the ratio of flux linked by an actual winding to the flux 
linked by a full-pitch, concentrated winding having the same number of turns. The winding 
factor is the product of a pitch factor (kp) and a breadth/distribution factor (kb) as shown in Eqn 
3-10. 

Ic “ Ic • Ic 1 
^pn 

Eqn 3-10 

The pitch factor accounts for the windings spanning a electrical degrees vice spanning a 
full 180® E as shown in Figure 16 [26]. 



Figure 16: Short-Pitch Coil 
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The pitch factor is the ratio of the flux produced by a short-pitch coil to the flux produced by a 
full-pitch coil. Short-pitching is an important means for eliminating harmonics and improving 
the power quality of the machine. The pitch factor can be derived with the final result shown in 
Eqn 3-11. 



n a ^ 


n-7i ^ 

= sin 


•sin 

.. 


V 2 ^ 


1 2 j 


where n = harmonic number 


Eqn 3-11 

The breadth factor explains the effect of the windings occupying a distribution or range 
of slots within a phase belt. A phase winding normally consists of numerous coils connected 
together linking flux slightly out of phase with each other as shown in Figure 17 [26]. 



The breadth factor can be derived either magnetically or geometrically to obtain Eqn 3-12. 


sm 


nm-y 


^bn “ 


^ 2 ; 


msin 


'^ny 


where n = harmonic number 

m = slots per pole per phase 
y = coil electrical angle 


Eqn 3-12 

In addition to estimating different winding effects, the geometry of the magnetic air gap 
must be represented. Field methods are utilized along with vector potential analysis to develop 
expressions that account for different magnetic gap geometries. Reference [19] contains detailed 
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derivations of the magnetic gap factor (kgn) for several magnet and slot configurations. The 
equation for the slotted stator, surface magnet configuration is shown in Eqn 3-13. 


, np-1 




_ r.2"P 


np 


np + 1 


• R 


np+1 j^^np+ 1 ^ 


np 


np 




where Rg = outer magnetic boundary Rg = outer boundary of magnet 

Rj = inner magnetic boundary R^ = inner boundary of magnet 


Eqn 3-13 


3.3.4 Flux and Voltage 

The primary significance of the magnetic flux linkage in a machine is that it induces 
voltage across a winding whenever the flux varies with time as explained through Faraday’s 
Law. The first step in the process is to determine the air gap flux density. The flux from the 
magnet poles crosses the air gap to the stator windings but some flux leaks along the way and 
this is accounted for using a leakage factor (Ki ~ 0.95 for surface magnets). In addition, the flux 
path is normally dominated by the air gap reluctance since the reluctance of the stator steel is 
much less than that in the air gap. However, a reluctance factor (Kr ~ 1.05 for surface magnets) 
is used to compensate for the small effects of the steel reluctance on the air gap flux. 

The presence of the slots in the stator also affects the air gap flux density because of the 
difference in permeance caused by the slots. The flux crossing the air gap in a slot region travels 
farther before reaching the highly permeable stator back iron. Carter’s coefficient (Kc) is used to 
account for this effect [17]. The air gap flux density is also affected by the magnet geometry in 
the air gap as previously described by Eqn 3-13. Since the magnet poles rotate north/south, the 
air gap flux density shape can be approximated as shown in Figure 18. This can be represented 
as a Fourier series using only odd components because of half-wave symmetry as shown in Eqn 
3-18. Overall, the air gap flux density is calculated using Eqn 3-14 through Eqn 3-18. 
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Figure 18: Air Gap Flux Density 


where Wg = average slot width 
Wj = tooth width 

'Cs = Ws + w, 

Eqn 3-14 

where = effective air gap 
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Eqn 3-15 
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where PC = permeance coefficient 

= flux concentration factor (A 


Eqn 3-16 
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where = recoil permeability 
= remnant flux density 


Eqn 3-17 
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0^ = magnet physical angle 
n = harmonic number 
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Eqn 3-18 

Once the flux density is known, the flux must be calculated. Given a machine uniform in 
the axial direction, the flux linked by a single, full-pitched coil which spans an angle from 0 to 
7i/p is represented by Eqn 3-19. Assuming Baux is sinusoidally distributed, the peak flux for this 
ideal coil is given by Eqn 3-20. 


O = 


n 

'•— 
P 
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where B = radial flux through coll 



flux 

P 


Eqn 3-19 


Eqn 3-20 

Given that there are Na coils in a stator phase winding and including all the real winding effects, 
the total flux linkage is shown in Eqn 3-21. Through Faraday’s Law, the back EMF for the 
machine is given by Eqn 3-22. 


oo 

?l(0) = y,k •sin(np0) 


n = 1 
n odd 


where 




Eqn 3-21 
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Ea= ^ Vjj Sin(np0) 


n = 1 

n odd 


where V„ = 


Eqn 3-22 

With permanent magnet excitation, the field cannot be controlled as in a wound-rotor 
machine so the number of turns in the stator phase windings must be chosen so that the machine 
EMF is close to the nominal system voltage. Another option is to use power electronics to 
convert the machine EMF to the system voltage thereby providing steady regulation. 

Another effect that is sometimes an issue is armature reaction. It is caused by current 
flowing in the stator windings which creates a magnetic field that tends to distort the magnetic 
field established by the permanent magnets. For surface-magnet machines, because the magnet 
recoil permeability is approximately one and the magnet height is large compared to the air gap, 
the armature reaction flux density is small. As long as the stator teeth are not highly saturated 
due to the permanent magnets, armature reaction is negligible [17], [18]. 


3.3.5 Machine Inductances 

In a slotted permanent magnet machine, there are three distinct components of 
inductance: air gap inductance, slot leakage inductance, and end-turn inductance. The most 
accurate means for calculating these parameters is finite element analysis but analytical methods 
provide almost as good results (within a couple of percent) and are used in this paper [33]. 

The air gap inductance is usually the largest portion of the total inductance and it is due 
to the interaction of the stator windings with the flux crossing the air gap. To calculate the air 
gap inductance, a full-pitch, concentrated winding carrying a current I is initially examined 
which leads to an air gap flux density shown in Eqn 3-23. 


oo 

®flux= X 


n = 1 

n odd 


where 
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Eqn 3-23 


47 







When this concept is expanded to polyphase windings with balanced operation, the air gap flux 
density becomes Eqn 3-24. 


oo 

®flux= Xi ®n Sin(np0) 


n = 1 

n odd 


where 




q 4 Pq 
2 nTt (g + hj„) 2p 


Eqn 3-24 

The flux can be found using equation Eqn 3-19 and the total flux linkage is X, = NaO. With all 
real winding effects included, the air gap inductance is then given by Eqn 3-25. 



q 4 P 0‘^s’^st'^a '^wn 


Eqn 3-25 

In addition to the air gap, the coil currents generate a magnetic field that crosses from one 
side of the slot to the other producing a slot leakage inductance. For calculating the slot leakage 
inductance, it is assumed that the slot is rectangular with slot depressions (Figure 13) which 
results in a slot permeance per unit length shown in Eqn 3-26 [17], [18], [19]. 


Pemi = 


1 hs 


3 w 


St 


w. 


Eqn 3-26 

Assuming m slots per pole per phase and a standard double layer winding, it can be shown that 
the slot leakage inductance is given by Eqn 3-27 through Eqn 3-29 [19]. 


Lgg = 2pLgjPerm- 


4.N/(m- N,p) + 2.N,p.Np 


(self) 


Eqn 3-27 


L^rn ~ 2- p • Lgj ■ Perm N^p • (mutual) 


Eqn 3-28 
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Eqn 3-29 

The end turn induetance is the smallest of the three components. It is created by the 
magnetic field that surrounds a coil after it leaves one slot and before it enters another slot. 

Since it is extremely difficult to accurately determine because of complex winding patterns, a 
rough approximation is used. It is assumed that the end turns are semi-circular with a radius 
equal to one-half the mean coil pitch. Using reference [17], the total end turn inductance per 
phase is shown in Eqn 3-30. 
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Eqn 3-30 

The total inductance for the phase is the sum of the three inductances, ignoring other small 
factors.. 


^ “ ^ag ^lot ^ ®0'^ 


Eqn 3-31 

\ 

3.3.6 Basic Losses 

Losses in a machine consist of core losses, conductor losses, friction and windage losses, 
and rotor losses. Rotor losses include magnet losses and retaining can losses, they require 
detailed waveform and harmonic analysis, and therefore are discussed later in Chapter 5. 

3.3.6.1 Core Losses 

High-speed generator stator core losses (per weight) can be greater than normal machines 
because of the higher frequencies. These losses are minimized by using laminated steels in the 
stator construction as discussed in 3.1.2 and by not generating frequencies that are too high. 

Core losses consist of hysteresis and eddy current losses. Hysteresis loss results from the 
steel not wanting to change magnetic state. As the flux density varies, the material traverses the 
B-H curve and energy is lost. Eddy current loss is also caused by the variation in flux density. 


49 


Electrical currents are induced within the ferromagnetic material and circulate dissipating power 
because of the resistivity of the material. 

Because there are usually various imperfections in materials, the best way to approximate 
core losses is to use empirical loss data. If the flux density is estimated for each part of a 
machine and the mass of the steel calculated, empirical core loss data can be used to estimate the 
total losses. Empirical data for M-19, 29 gauge material is obtained as shown in Figure 19. 
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Figure 19: Core Loss Data 

An exponential curve fit is then applied to the data to obtain an equation for estimating the core 
losses (Eqn 3-32) [9], [35], [36]. 

£ 15 £ f 


Eqn 3-32 


The values for the bases and exponents in Eqn 3-32 are listed in Table 8. 
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Table 8: Core Loss Parameters 




Base Power (Po) 

36.79 W/lb 

Base Flux Density (Bo) 

LOT 

Flux Density Exponent (sb) 

2.12 

Base Frequency (fo) 

1000 Hz 

Frequency Exponent (sf) 

1.68 


3.3.6.2 Conductor Losses 


Conductor losses arise from the current flowing in the stator windings. The resistance of 
the windings is calculated using Eqn 3-8 from section 3.3.2. The conductor losses are then found 
using the traditional power equation for a resistance (Eqn 3-33). 


P. = qi. R, 


Eqn 3-33 

3.3.6.3 Friction & Windage Losses 

For rotors operating at high-speed, friction and windage in air can cause losses which 
result in inefficiency and heat production. These losses are calculated using the power necessary 
to overcome the drag resistance of a rotating cylinder as given by Eqn 3-34 [32]. 


3 4 

Pwind “^f^ Pair ® R '^t 

where Cj = friction coefficient 
Pgir = density of air 


Eqn 3-34 

The friction coefficient depends on numerous factors such as surface roughness and flow region. 
Since the air gap is a small annulus and the rotor is spinning at high speed, it is assumed that the 
air in the gap is in the turbulent region. Therefore, the coefficient of friction can be 
approximated by Eqn 3-35 [34]. 

C/-0.0725 Rey-^'® 
where Rey = Reynold's Number 


Eqn 3-35 
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3.4 Machine Sizing Methods 


3.4.1 Basic Sizing Method 

Whenever a machine is being designed, it is important to perform some back-of-the- 
envelope calculations to gain insight into initial sizing estimates. Most generators are 
constrained by two competing design parameters, air gap magnetic shear stress and rotor tip 
speed. Air gap magnetic shear stress (t) is the magnetic shear force developed per unit gap area 
and is constrained by magnetic design and thermal management [8]. It is proportional to the 
product of the surface current density and magnetic flux density as shown in Eqn 3-36. 

T o<= Kz-Bg 

where f = sheaf stress (psi) 

Kz = surface current density 
Bg = air gap flux density 


Eqn 3-36 

Typical values for air gap shear stress for different types of generators are shown in Table 9 [8], 
[11], [16], [18], [31]. For the basic sizing calculations, 15 psi is assumed since the generator is a 
large liquid-cooled machine. 


Table 9: Air Gap Shear Stress Values 




Small Air-Cooled 

1 - 5 psi 

Large Air-Cooled 

■ihbieeisishhh 

Large Liquid-Cooled 


High-Temperature Superconducting 

■■EDBSIsSHHH 


Rotor tip speed is discussed previously in section 3.2.2 and given by Eqn 3-1 with the 
assumed limit of 200 m/s. The fundamental machine power equation is utilized to derive the 
rotor radius and stack length of the machine (Eqn 3-37). 

P=2-7tRLgt'cv^ip 

where R = rotor radius 
Lg( = stack length 


Eqn 3-37 
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In order to simplify the equation, the L/D ratio is substituted in for Lsf Using the air gap shear 
stress, the rotor tip speed limit, and the power rating of the machine (16 MW), the basic power 
equation is iterated to obtain a final L/D ratio, rotor radius, and stack length while matching the 
desired rotational speed of the machine (13,000 RPM). Using an air gap flux density of 0.8 T, a 
pole pair value of 3, a slot height of 15 mm, and a slot fill fraction of 0.5, the frequency and 
current density of the machine are found. The detailed MATLAB code is contained in Appendix 
B with the results shown below. 


Basic Machine Design 

Input Parameters: 

Power = 16000.0 kW 

Shear Stress = 

15.0 psi 

L/D Ratio = 

2.85 

Tip Speed = 

200.0 m/s 

Pole Pairs = 

3.0 

Air Gap Bg = 

0.80 T 

Output: 

Rotor Radius = 

0.147 m 

Stack Length = 

0.838 m 

Speed = 

13000 RPM 

Frequency = 

650.0 Hz 

Ja 

1757.67 A/cm2 




3.4.2 Detailed Sizing Method One 

Once basic sizing of the machine is complete, in-depth analysis is conducted to ascertain 
the overall performance and scale of the 16 MW generator. Two detailed sizing methods are 
developed using MATLAB code. To help gain an understanding for current PM machine 
designs, numerous articles, transactions, and proceedings are examined and a database is 
developed containing comprehensive specifications for over 20 different PM machines (see 
Appendix C). This database is used to help develop input parameters for the first method (Table 
10), although many of these values are obtained after numerous iterations through the code. 


Table 10: Input Parameters for Sizing Method 1 


Jl^i^iiefer 


if -:^rani!efiBr 


Required Power 

16 MW 

Number of Phases 

3 

Rotational Speed 

13,000 RPM 

Number of Slots 

36 

Power Factor Angle 

0 deg 

Slots Short-pitched 

1 

Rotor Radius 

0.147 m 

Peripheral Tooth Fraction 

0.5 

Magnet Height 

25 mm 

Slot Depth 

25 mm 

Stack Length 

0.838 m 

Slot Depression Depth 

0.5 mm 

Pole Pairs (p) 

3 

Slot Depression Width 

N/A 

Magnet Br 

1.2 T 

Stator Back Iron Ratio 

0.7 

Magnet Angle 


Turns per Coil 

1 

Magnet Skew Angle 


Slot Fill Fraction 

0.5 
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Air Gap _ 

Steel Density 
Magnet Density 


4 mm 
7700 kg/m 
7400 kg/m 


Winding Conductivity 

Conductor Density 


6.0 X 10^ S/m 
8900 kg/m" 


Once the input parameters are entered, the first step in sizing is to generate the geometry 


of the machine. This includes determining the number of slots per pole per phase (Eqn 3-5), 
number of armature turns (Eqn 3-38), tooth width, slot dimensions, stator back iron dimensions, 
coil pitch, and winding end turn geometry. 


= 2-pmN^ 

where N^, = Turns per coil 

Ng assumes each slot has 2 half coils 


Eqn 3-38 

Next, the electrical frequency and rotor surface speed are determined using Eqn 3-land Eqn 
3-39. 


- PJl 
60 

where 


CO=2-7tf 0)=p(0jjj 

CO = electrical frequency (rad/sec) 
cogi = mechanical frequency (rad/sec) 


Eqn 3-39 

Winding, skew, and magnetic gap factors are then estimated as discussed in sections 3.2.3 and 
3.3.3. The magnitude of the air gap magnetic flux density (Bg) is determined accounting for 
slots, varying reluctances, and flux leakage per Eqn 3-14 through Eqn 3-17. 

The magnetic flux and back voltage magnitudes can then be calculated. The fundamental 
component of the magnetic flux links the stator windings to create useful voltage. Therefore, 
only the fundamental components of Eqn 3-21 and Eqn 3-22 are used to determine the internal 
voltage of the generator as shown in Eqn 3-40 and Eqn 3-41. 


B 1 = —B - k^j-sin 

71 ^ ^ 




P% 

V 2 y 


Eqn 3-40 
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Ea - 



P 


Eqn 3-41 

Utilizing equations from section 3.3.5, the machine inductances and reactances are 
established. The lengths, volumes, and masses of components and the overall generator are 
calculated using basic geometric equations and relationships. A 15% service mass fraction is 
added to the total mass estimate to account for the additional services associated with large 
liquid-cooled machines [31]. Once the mass of each of the stator parts is known, the core losses 
are estimated in accordance with section 3.3.6.1. 

The terminal voltage and current of the machine must then be calculated accounting for 
conductor losses and windage losses. The vector relationship (Figure 20) between terminal 
voltage (Va), internal voltage (Ea), and the synchronous reactance voltage drop is utilized to 
obtain Eqn 3-42. The armature resistance is usually ignored because it is much smaller than the 
synchronous reactance. 



Figure 20: Voltage Vector Relationship 


(Xs-Ia-cosvf - Xg la-simi; 

Eqn 3-42 
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The machine efficiency (r|) is then easily obtained. Appendix D contains the detailed MATLAB 

code and the results are shown below. 

PM Machine Design, Version 1: Surface Magnet, Slotted Stator 
Machine Size: 


Machine Diameter = 

0.472 m 

Machine Length = 

1.003 m 

Rotor radius = 

0.147 m 

Active length = 

0.838 m 

Slot Avg Width = 

16.493 mm 

Slot Height = 

25.000 mm 

Back Iron Thick = 
Machine Ratings: 

34.300 mm 

Tooth Width = 

15.403 mm 

Power Rating = 

16000.0 kW 

Speed = 

13000 RPM 

Va (RMS) = 

2341 V 

Current = 

2267.5 A 

Ea (RMS) = 

2925 V 

Arm Resistance = 

0.00526 ohm 

Synch Reactance = 

0.768 ohm 

Synch Induct = 

0.188 mH 

Stator Cur Den = 

2199.7 A/cm2 Tip Speed = 

200 m/s 

Efficiency = 

0.992 

Power Factor = 

1.000 

Phases = 

Stator Parameters: 

3 

Frequency = 

650.0 Hz 

Number of Slots = 

36 

Num Arm Turns = 

12 

Breadth Factor = 

0.966 

Pitch Factor = 

0.966 

Tooth Flux Den = 
Slots/pole/phase = 

Rotor Parameters: 

1.59 T 

2.00 

Back Iron = 

1.14T 

Magnet Height = 

25.00 mm 

Magnet Angle = 

50.0 degm 

Air gap = 

4.00 mm 

Pole Pairs = 

3 

Magnet Remanence = 

1.20 T 

Aig Gap Bg = 

0.80 T 

Magnet Factor = 

Machine Losses: 

0.949 

Skew Factor = 

0.989 

Core Loss = 

11.6kW 

Armature Loss = 

81.2 kW 

Windage Loss = 

Machine Weights: 

30.7 kW 

Rotor Loss = 

TBDkW 

Core = 

396.49 kg 

Shaft = 

438.05 kg 

Magnet = 

64.74 kg 

Armature = 

89.58 kg 

Services = 

148.33 kg 

Total = 

1137.18 kg 


3.4.3 Detailed Sizing Method Two 

In order to provide a check on the methodology of the previous sizing procedure, a 
second MATLAB code is constructed. The second method is developed using a combination of 
processes from references [2], [17], [18], [19], [26], and [37]. 
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Table 11: Input Parameters for Sizing Method 2 



Value 




16^ 

Current Density (Ja) 

■rWlllMBBII 

Rotational Speed 

13,000 RPM 

Number of Phases 

3 

Power Factor Angle 

0 deg 

Slots/Pole/Phase (m) 

2 

Maximum Vtip 

200 m/s 

Slots Short-pitched 

1 

Bsat 

1.65 T 

Avg Slot Width (Ws) 

16.0 mm 

Pole Pairs (p) 

3 

Slot Depth 

25 mm 

Magnet Br 

1.2 T 

Slot Depression Depth 

0.5 mm 

Permeance Coefficient (PC) 

5.74 

Slot Depression Width 

N/A 

Magnet Skew Angle 


Turns per Coil 

1 

Air Gap 

4 mm 

Slot Fill Fraction 

0.5 



Winding Conductivity 

gfifUaggM 



Conductor Density 



Many of the same equations used in the first method are utilized but the process differs in some 
of the input parameters and calculations. The second sizing method inputs are shown in Table 
11 with different entries from method 1 being Vtip, Bsat, Ws, m, Ja, and PC. 

First, given the maximum Vtip (200 m/s) and the rotational speed (13,000 RPM), the 
electrical frequency and rotor radius are computed. Next, the winding and skew factors are 
determined similar to the first sizing method. The magnet dimensions (h^ and Gm) are then 
determined along with the air gap flux density. Using Bsat, the PC, and Eqn 3-43 through Eqn 
3-46, the magnet arc width and height are iterated until the tooth width (wt) equals the average 
slot width (Ws). This differs from the first method where the peripheral tooth fraction is used to 
set the tooth width equal to the slot top width (Wst). 


where = magnet pitch coverage coefficient 

Eqn 3-43 

^ m“ S 

Eqn 3-44 
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Eqn 3-46 

The permeance coefficient (PC) normally varies between 5-15 with higher PCs only 
used for high performance space and aircraft applications [18], [37]. A PC value of 4 - 6 is 
typical for large PM machines where the air gap flux density is approximately 60 - 80% of the 
remnant flux density [17]. The PC is input as 5.74 to coincide with the PC calculated in the first 
sizing method. 

Once the magnet dimensions are known, the geometry of the machine is generated using 
similar equations as the first method except for the stator core back iron depth. In the first 
method, the input variable “stator back iron ratio” is used whereas in the second method, the 
saturation flux density is applied as shown in Eqn 3-47. 

_ 7tDgem Bg 
"" 4.p 


Eqn 3-47 

Utilizing a power balancing procedure, the stack length, terminal voltage, and terminal 
current are iterated to obtain a complete design. Core losses, conductor losses, and windage 
losses are calculated in accordance with sections 3.3.6.1, 3.3.6.2, and 3.3.6.3. The power 
crossing the air gap of the machine is determined using Eqn 3-48 [2], [37]. 

where kg = electrical power waveform factor 
kj = current waveform factor 


Eqn 3-48 
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The electrical power waveform factor is defined in Eqn 3-49 and the current waveform factor is 
shown in Eqn 3-50. These factors depend on the hack EMF and current waveforms generated in 
the machine and a detailed discussion of these factors is in reference [2]. 


k 


e 


T 


^pk'^pk 


^i = 


^pk 


rms 


Eqn 3-49 


Eqn 3-50 

During the balancing procedure, the machine inductances and reactances are established 

per section 3.3.5. The lengths, volumes, and masses of components and the overall generator are 

calculated using basic geometric equations and relationships similar to the first sizing method. 

Once the balancing procedure is complete, the efficiency of the machine is determined. 

Appendix E contains the detailed MATLAB code and the results are shown below. 

PM Machine Design, Version 2: Surface Magnet, Slotted Stator 
Machine Size: 


Machine Diameter = 

0.484 m 

Machine Length = 

0.979 m 

Rotor radius = 

0.147 m 

Active length = 

0.813 m 

Slot Avg Width = 

16.000 mm 

Slot Height = 

25.000 mm 

Back Iron Thick = 

39.587 mm 

Tooth Width = 

15.863 mm 

Machine Ratings: 
Power Rating = 

16000.0 kW 

Speed = 

13000 RPM 

Va (RMS) = 

2604 V 

Current = 

2064.1 A 

Ea (RMS) = 

3007 V 

Arm Resistance = 

0.00534 ohm 

Synch Reactance = 

0.723 ohm 

Synch Induct = 

0.177 mH 

Stator Cur Den = 

2064.1 A/cm2 

Tip Speed = 

200 m/s 

Efficiency = 

0.993 

Power Factor = 

1.000 

Phases = 

3 

Frequency = 

650.0 Hz 

Stator Parameters: 
Number of Slots = 

36 

Num Arm Turns = 

12 

Breadth Factor = 

0.966 

Pitch Factor = 

0.966 

Tooth Flux Den = 

1.70 T 

Back Iron = 

1.05 T 

Slots/pole/phase = 
Rotor Parameters: 

2.00 



Magnet Height = 

26.17 mm 

Magnet Angle = 

47.7 degm 

Air gap = 

4.00 mm 

Pole Pairs = 

3 

Magnet Remanence = 

1.20 T 

Aig Gap Bg = 

0.85 T 
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Magnet Factor = 
Machine Losses: 

0.958 

Skew Factor = 

0.989 

Core Loss = 

11.9kW 

Armature Loss = 

77.5 kW 

Windage Loss = 
Machine Weights: 

29.7 kW 

Rotor Loss = 

TBDkW 

Core = 

438.62 kg 

Shaft = 

424.57 kg 

Magnet = 

62.89 kg 

Armature = 

85.50 kg 

Services = 

151.74 kg 

Total = 

1163.32 kg 


3.4.4 Comparison of Methods 

Both sizing methods produce similar sized generators with the results agreeing within 1 - 
5% on most parameters. The main reasons for the differences are twofold. First, the calculations 
for the magnet dimensions are different. In the first method these parameters are input whereas 
in the second they are determined through an iterative process using material properties. Second, 
the procedure for determining the slot and tooth widths differs slightly. 

Overall, both methods underestimate the overall dimensions (length and diameter) and 
weight. This is because the sizing programs do not include calculations for portions such as 
stracture, frames, mounts, and maintenance access. These pieces are added in section 6.1 to get 
the total sizes and weights of the PM generators, specifically including: 

• 50% weight factor added to total weight 

• 0.8 m added to overall length and 0.4 m added to overall diameter 

In later chapters when waveforms are examined and machine optimization occurs, the first sizing 
method is mostly employed with the second method used for verification. 
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Chapter 4 Power Electronics and Conversion 


4.1 Background 

In conjunction with employing a PM generator, a DC bus architecture is one of the 
preferred schemes for the future [5]. DC power distribution systems can offer an advantage in 
size and a weight over high-power AC systems [3]. Therefore, a power electronics module is 
required to rectify the AC generator output and then convert it to the appropriate DC distribution 
voltage of 700 VDC. However, conversion of the high-frequency generator AC output to DC 
generates harmonic voltages and currents on the AC side and ripple on the DC side. These 
harmonics can be reduced through proper design of the conversion module. 

Since the generator is a permanent magnet machine, the field cannot be controlled like in 
a wound-rotor machine. The power electronics module (PEM) is needed not only to convert the 
machine EMF to the DC system voltage but also provide steady regulation. Also, since the 
frequency of the generator output is higher than the normal 60 Hz systems, the filter 
requirements for the PEM are less since the components can be smaller. 

Many options exist for the t 3 q)e and model of the PEM, and a combination of a rectifier 
and DC-DC buck converter are chosen. This design provides a two-step conversion from the 
generator to the DC bus and allows for good DC bus voltage regulation while reducing 
harmonics reflected back to the generator. A block diagram of the system is shown in Figure 21. 
The calculations for the rectifier and buck converter are performed for the ideal case and ignore 
parasitics such as equivalent series resistance, inductance, and capacitance. 
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4.2 Rectification 

The rectifier stage converts the AC voltage from the PM generator to a DC voltage. 

There are two main choices when designing the rectification stage: active or passive rectification 
and series or polyphase rectification. Active rectification uses controllable components, such as 
thyristors or MOSFETs, to actively rectify the AC voltage. This increases the cost and 
complexity of the rectifier but provides voltage regulation and control. Passive rectification 
consists of diode bridges to convert the AC to DC while relying on some other circuit to perform 
the voltage regulation. A passive rectification scheme is used in this paper because of its 
simplicity, potential for lower input harmonics, and lower cost. 

Series rectification is when each phase is rectified independently from the others and the 
output DC voltage is obtained by a series connection of each single phase rectifier. Each phase 
conducts during the full period resulting in high power losses and higher stresses on components. 
Also, the output DC voltage is extremely high since each rectified phase voltage is added 
together. Parallel rectification is when the lines with the most positive and negative 
instantaneous voltages provide the forward bias to turn on two diodes. The output voltage 
waveform corresponds to the instantaneous difference between two line voltages. For sinusoidal 
voltage supplies, the equation for the output DC voltage is approximately Eqn 4-1 [38]. 



where Vg = peak amplitude of phase voltage 

Eqn 4-1 

The losses are lower with parallel rectification and its output voltage provides a better 
match with the buck converter. Therefore, a parallel rectification scheme is used in this paper. 
For the initial PM generator from section 3.4.2, the rectifier is shown in Figure 22. 
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Figure 22: Basic 3-phase Rectifier 

The DC-side inductor and capacitor are necessary to obtain a good DC output voltage 
with low ripple (large C) and an acceptable AC line current waveform (large L) [39]. The goal is 
to produce maximum power factor with minimal line-current harmonics. A power factor of 0.96 
with low line harmonics can be achieved by selecting the inductor and capacitor values using 
Eqn 4-2 and Eqn 4-3 [40]. 



where Lor,=0.10 

= nominal l-n RMS source voltage 
■ref = P ref/Vref (RMS current from V rgf) 
P^gf = nominal output power 
w = nominal source frequency 



w ^ 

^ref ) 


where Cgr, = 100.0 

V^gf = nominal l-n RMS source voltage 
•ref = P ref/Vref (RMS current from V ^ef) 
P^gf = nominal output power 
f^gi = nominal source frequency 


Eqn 4-2 


Eqn 4-3 
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For the initial PM generator from section 3.4.2 where Vref = 2925 V, Pref = 16 MW, and fref = 650 
Hz, the values for the DC-side inductor and capacitor are Lr = 82.3 pH and Cr = 287.7 mF. 

These values will change as the system is optimized. 


4.3 DC-DC Conversion 

After the rectifier converts the AC voltage to DC voltage, the DC voltage must be 
stepped down to the bus voltage of 700 VDC. A high-power buck converter is ideal for this 
application. The power circuit has a basic topology to which other circuit components are added 
to perform functions such as filtering and over-voltage protection [41]. The overall goal is to 
alter electrical energy provided by an input system to that required by an output system. 


The system consists of a two-stage input filter, a converter, and a low-pass output filter as 
shown in Figure 23. Duty ratio control using a feedback controller could be used to provide DC 
bus voltage regulation during normal operation. 



Figure 23: Basic Buck Converter 

4.3.1 Buck Converter 

The buck converter is considered a high-frequency DC-DC switching converter. The 
difference between the switching frequency and the frequency of the external waveforms allows 
the use of low-pass filters in the input and output to help improve performance [41]. The buck 
converter is designed to take a high DC voltage and reduce it to a lower DC voltage under 
varying loads and over a range of input voltages. The level of voltage reduction is controlled by 
the high frequency switching. The switch opens and closes at a specific frequency with the ratio 
of on-time to the period defined as D, the duty ratio. 









For the buck converter, the voltage conversion ratio is dictated by Eqn 4-4. 


\)Ut 



= D 


^out 1 


Eqn 4-4 

It is assumed that the required output power from the generator and PEM can be as high as the 
rated 16 MW and as low as 0.1 MW. Based on this, output currents and resistive load values are 
determined and are listed in Table 12. 


Table 12: Buck Converter Load Values 





Load Resistance 

o.i MW 

700 VDC 

142.9 A 

4.90 Q 

16 MW 

700 VDC 

22,857.1 A 

0.031 


The switching frequency for the converter is selected to be 50 kHz (T = 20 psec). This is 
sufficiently fast to ensure passive components are not excessively large but also not too high 
speed so that there are not excessive switching losses. Additionally, current switching ratings for 
advanced, high-power IGBT devices are limited to the hundred kilohertz range. By selecting 50 
kHz, this provides a margin below this maximum switching frequency and ensures devices are 
presently available or will be in the near future. 


4.3.2 Output Filter 

A low-pass output filter is required to help reduce the output voltage ripple to within the 
required specification listed in Table 2. The output voltage from the converter normally 
fluctuates between 0 and Vqut, so the filter, consisting of an inductor and capacitor, is used to 
minimize this problem. It is important that the comer frequency of the filter be much lower that 
the switching frequency (fc « fsw) so that the switching frequency ripple is eliminated in the 
output voltage [42]. Normally, a factor of at least lOOx is used [43], so to ensure an adequate 
margin for parasitics and fulfill the strict DC bus voltage ripple requirement, a factor of 400 is 
selected so that fc ~ 125 Hz. For an LC-filter, Eqn 4-5 holds. 

f c = ^ 

InyfUc 


Eqn 4-5 
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Therefore, given the value selected for fc: 


L C= 1.6210 ^ 


Eqn 4-6 

With this equation the output filter inductor and capacitor are sized using a buck 
converter sizing program [44]. It is assumed that an inductor current ripple of 0.1% at maximum 
power (16 MW) is satisfactory for sizing the inductor and Eqn 4-7 is used where Vp = 1.9 V, the 
forward-bias diode voltage [44]. 
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■(y: 







I inm^ F I 
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Eqn 4-7 

For the initial generator of section 3.4.2, this produces an inductor size of 551 pH. Using this 
inductance along with Eqn 4-6, a capacitor value of 2.94 mF is calculated. 

In order to validate the filter’s performance, a rough estimate for the voltage ripple is 
calculated. It is assumed that the ripple component of inductor current flows through the 
capacitor and the average component of the output current flows through the load resistor [42]. 
Therefore, the peak-to-peak voltage ripple can be written as shown in Eqn 4-8. 

AV - - 1 - — — 

C ~ C 2 2 2 


Eqn 4-8 

The output voltage ripple is estimated to be 0.019 V which is within the specification listed in 
Table 2. This overall process is used later to design the converter output filter as the PM 
generator system is optimized. 


4.3.3 Input Filter 

The input filter is needed to prevent the converter switching fluctuations from affecting 
the voltage and current from the rectifier. The filter selected is a two-stage filter with an L-C 
filter and an R-C damping leg. From the switching frequency, the angular frequency is © = 
3.142x10*’ rad/sec. In order to examine worst case, the 16 MW, 22,857 A-output is used, and 
with the corresponding duty ratio of 0.1023 for the generator of section 3.4.2, an average current 
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of 2,338.3 A is calculated. Since the fundamental of the ripple is attenuated much less than the 
higher harmonics, the goal is to suppress the fundamental of the ripple along with some margin. 
The fundamental of the average current is calculated in Eqn 4-9. 

4 4 

iconv\ ~ ^ iconv — * (2338.3) = 29112Amps 

K K 

Eqn 4-9 

It is assumed that the buck converter input current ripple must be kept very small (< 10 A at 
maximum power) to ensure no detrimental effects on the rectifier. Therefore, at the 50 kHz 
switching frequency, the maximum gain of filter is G as calculated in Eqn 4-10. 

lO.OAmp^ 

^ ; ripple w 

-2-= 0.00168 

iconvl 2977.2 


Eqn 4-10 

In order to account for the contribution of higher harmonics and non-idealities associated with 
actual components, a 50% margin is used for G resulting in G = 2.52x10'^. 

The output impedance of the input filter is also a concern because it must be sufficiently 
small so that it does not adversely affect the converter’s performance. It is assumed that the 
output voltage ripple of the filter must be less than 10 V at maximum power. From this, a value 
of Zo at 50 kHz is calculated (Eqn 4-11). 

10.0 V 

Z „ = —-— = 0.00168 
° 2977.2 


Eqn 4-11 

Using the same margin of 50% for the same reasons, the output impedance is 1.12x10"^. 

Initially, the L-C low pass filter portion is examined to calculate values for Lf and Cf. 
Then, a well-designed filter has Cb ~ 10*Cf, and Rf is determined using reference [39]. The 
calculations are listed below. 
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-0 

* rvvv\ 


Lf 

.t 

Iin 

Lf 

♦ 


” Zo C 

--o 

^ Cf = 

= G 

L- 


Zo = 


Zo = 


Zl+ 


s-L 


1 + s LC 


G = 


Hn = 

JiL 

*out 


Zc + Zj^ 


out 




lout 


1.1210 ^ =^ 


1 + s LC 


2.52-10 ^ = 


1 + s -L-C 


1.1210 1.1210 ^-s^-L-Cr: s-L 


2.5210 ^ + 2.5210 ^-s^-L-C::! 


s-L= 0.444 
Lf = 1.415uH 


s •L-C= 395.83 

s-C= 891.51 Cf=2.84mF 


Eqn 4-12 


Cj, = n-Cf = lO-Cf = 28.4-mF 


Ro= It 


0.0223 


R 


f _ (2 + n)-(4 + 3n) 

^0 J 2-n^-(4+n) 


Rf = 0.215 -Rq = 4.79 -10 
Choose 10*Rf to be safe 


Rf = 4.79 mQ 


Eqn 4-13 

The performance of the input filter is determined by examining the filter’s transfer 
function and verifying it has sufficient attenuation at the switching frequency and also minimal 
peaking. The transfer function for the filter is given in Eqn 4-14 [45]. 
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'RfCb's) 

1 + 1 

f 3^ 

iLfCfCb-RfS j 

1 + 

Lf(Cf + C 

+ (RfCb's) - 1 - 1 


Eqn 4-14 

Figure 24 shows the Bode plot for the initial input filter and it operates satisfaetorily (Appendix 
F contains the detailed MATLAB code). This overall procedure for designing the converter 
input filter is utilized later to as the PM generator system is optimized. 



Figure 24: Bode Plot for Converter Input Filter 

4.3.4 Converter Control 

One control technique that can be used is duty ratio control utilizing a feedback 
controller. This arrangement is presented as an illustration of one control methodology for the 
power electronics module but implementation is beyond the scope of this paper. 

The control scheme involves negative feedback control using linearized, averaged state- 
space models. Figure 25 shows a block diagram of this where K represents the controller and H 
represents the buck converter. It is designed so that the closed-loop system is stable and well 
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damped for the input voltages and output loads. Voltage-mode control (duty-ratio control) is 
implemented using a Proportional Integral Derivative (PID) controller. 



Figure 25: Block Diagram of Control Scheme 


The majority of control systems used in industry applications are based on this type of 
controller. In general, proportional and integral control can be applied independently or together 
(P, I, or PI), while derivative control can only be applied in combination with one or both of the 
other two (PD, PID). Testing of a PID controller normally involves plotting the system response 
to a simple input, such as a step input, and adjusting the gains based on the output. 

The first step in developing the control scheme is to obtain a linearized, averaged model 
for the buck converter and derive the transfer function from perturbation in duty ratio to 
perturbation in output voltage. The transfer function represents H in Figure 25. Next the PID 
controlled is developed, consisting of three gains, Kp, Ki, and Kd, and this represents K in Figure 
25. The H and K are then combined to form the complete transfer function for the system. 

Initially, the state variables are assigned (it and vc) and the state-space equations are 
determined (neglecting the input filter). 

L—iL = (^l - vc)-q(t) - vc (l - q(t)) C~vc = fiL" “ f‘L” V ~ 

at at \ K J \ K J 



•(vj-qd) - Vc) 



Eqn 4-15 

Then, the circuit is averaged and it is assumed that the state variables have slow variation and 
small ripple, and d = avg(q(t)). All of the variables are now averaged quantities and Eqn 4-15 
becomes Eqn 4-16. 
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Eqn 4-16 


5‘l -’d 


d 1 
dt ^ C 


( vc^ 

Ry 


The variables are linearized resulting in the equations shown in Eqn 4-17. 


vc= Vc+ v'c 


Linearized Equations 


d = D + d' 


vi = Vi + v'j 


dt ^ L ^ L L ^ 


d , 1 1 , 

—vr = —It -Vp 

dt C ^ R-C ^ 


Eqn 4-17 

Since vc is the output voltage, utilizing the Laplace operator and solving for the transfer function 
yields Eqn 4-18. 


2 1 1 ^ 

s +-s +- 

RC LC 7 


out 


H = 


D , Vi 

-V 1 +-d 

LC ^ LC 


V, 


out LC 


d' 2 1 1 

s +-s +- 

RC LC 


For the controller the transfer function is Eqn 4-19. 

Kj where e' is 


£ 

e' 


K = — = Kp + %s + — 


Eqn 4-18 


K = 


Kj-)'S + Kp'S + Kj 


Eqn 4-19 


Overall, the total transfer function is then Eqn 4-20. 
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= HK 


^ out 
e’ 

v'out 


e' 3 1 2 1 

s +-s +-s 

RC LC 


Eqn 4-20 

Using the total system transfer function, the gains on the controller are adjusted to obtain 
the desired steady-state and transient operation. This type of system can readily be designed and 
optimized to provide the required overall performance, but this is beyond the scope of this paper. 


4.4 Conversion Losses 

An integral part of a DC-DC converter’s performance is how efficiently it operates 
during the energy conversion process. Losses that occur include winding losses due to 
conduction, switching losses, and inductor core losses. A detailed analysis of all of these losses 
is beyond the scope of this paper and they are only described qualitatively and estimated using 
analytical equations. 

Practical magnetic components, such as inductors, exhibit losses generated by winding 
conduction, hysteresis, and eddy currents. Resistive components in the converter as well as 
parasitics (equivalent series resistance) produce I^R losses. Ideal inductors and capacitors are 
employed in this paper and parasitic losses are not calculated. 

The transistors and diodes have conduction losses based on their on-state resistances. 
The equations to estimate these losses are shown in Eqn 4-2land Eqn 4-22. Lastly, the 
transistors have switching losses associated with the transistors changing states. If the voltage 
and current waveforms are assumed to rise and fall approximately linearly, the switching losses 
are determined as shown in Figure 26. 



where 1^ = average current 
D = duty cycle 

Rj ,3 = IGBT on-state resistance 

Eqn 4-21 
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•■diode = (• - ••)lL">^d 

where II = average current 
D = duty cycle 

Rj, = diode on-state resistance 

Eqn 4-22 



4.5 Component Sizes and Weights 

The PEM detailed in the previous seetions consists of IGBTs, diodes, capacitors, 
inductors, and resistors. The sizes and weights of the resistors are neglected because they are 
significantly smaller than the other components. To ascertain reasonable estimates for the other 
devices, numerous references are utilized [15], [46], [47], [48], [49], [50], [51], [52], [53], [54]. 
In addition, a 30% service fraction is added to the total component mass and volume calculations 
to account for heat sinks, cooling, etc. 

The power transistor is the main switching device in the converter circuit. In many 
applications, IGBTs are used because they have the current density and low loss of bipolar 
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transistors with the high-speed and high input impedance of MOSFETs [54]. Power diodes are 
needed to complete the buck converter circuit and provide proper operation. Table 13 lists the 
assumed device parameters for the IGBTs and diodes used. The IGBT and diode values are not 
from exact devices from particular manufacturers but rather represent data assembled from 
numerous sources. In order to achieve the high power levels, the semiconductor devices must be 
connected in parallel. Specifically, since the IGBTs and diodes are rated for 1,200A each, 19 
modules (each rated for 845 kW) are connected in parallel to achieve the required 22.5 kA. This 
can be readily accomplished using IGBTs provided the switching is properly controlled [55]. 

The inductors and capacitors are needed to perform energy storage and filtering 
functions. In high-power applications they are usually the dominant components in terms of 
mass and volume. For estimating their contribution to the converter, component weight and 
volume energy densities are utilized. Table 13 lists the assumed characteristics for the inductor 
and capacitor sizing estimations. 


Table 13: Power Electronics Module Component Characteristics 



Maximum Collector Emitter Voltage 

6500 V 

Maximum Collector Current 

1200 A 

Mass (including heat sink) 

1200 g 

Volume (including heat sink) 


Turn On / Rise Time 

450 ns 

Turn Off/Fall Time 

400 ns 

On-State Resistance 

0.06 Q 


1.9 V 

Maximum Switching Frequency 

-100 kHz 


Peak Blocking Voltage 

6500 V 

Average Forward Current 

1200 A 

Mass (including heat sink) 

450 g 

Volume (including heat sink) 


On-State Resistance (estimated) 

0.0073 n 


1.0 V 


Mass Density 


mmsmsmmmmmm 

8 J/cc 



Mass Density 

3 kg/J 

Volume Density 

0.08 J/cc 
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In addition to using the component values to size the power conversion module, a weight factor 
(l.Ox) and a volume factor (0.20x) are employed to account for portions such as structure, 
frames, mounts, and maintenance access. 
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Chapter 5 Waveforms, Models, and Machine/Module Optimization 


5.1 Initial Generator Waveforms 

The PM generator produces back EMF waveforms that are dependent on a number of 
factors as discussed in Chapter 3. The goal is to produce a voltage waveform that closely 
resembles a sinusoidal waveform with a low total harmonic distortion (THD) because this results 
in minimal harmonic content which reduces losses in the machine. THD is a measure of the 
distortion in a waveform caused by undesirable frequency components. It is calculated as shown 
in Eqn 5-1. 

THD = 

Eqn 5-1 

The back EMF waveforms are generated using Eqn 3-18, Eqn 3-19, and Eqn 3-22 which 
are shown again below. 



oo 

b( 0) = ^ Bj^ sin(np0) 


n = 1 
n odd 


where 


4 , . "P0m . 

B =-B k -sin — sin 

" nil ^ ^ \ 2 ) 

0^ = magnet physical angle 

n = harmonic number 



Eqn 5-2 


oo 

?i(0)= ^ A-n •sin(np0) 


n = l 
n odd 


where 


P 


Eqn 5-3 
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oo 


Ea= ^ Vjj Sin(np0) 


n = 1 

n odd 


where v = — A, = (On X 
at 


Eqn 5-4 

MATLAB code is developed to create the generator EMF waveforms as well as 
determine the harmonic content and THD of the generator output waveform. This code is 
contained in 0. For the initial PM generator designed in section 3.4, the air gap flux density, 
back FMF waveform, and harmonic content (line-neutral) are shown in Figure 27, Figure 28, and 
Figure 29. As a frame of reference, MIL-STD-1399 requires Type 160 Flz AC systems to have a 
total THD of less than 5 % line-line [56]. This limit does not directly apply to the PM generator 
AC bus since it is being rectified to DC before supplying loads. However, it is used as a 
guideline because it provides a reasonable limit to ensure proper rectification and performance 
on the AC side. 














' Rotor 



Figure 28: Initial Generator EMF Waveform 



Figure 29: Initial Generator Harmonic Content 






























5.2 Rotational Stress and Retaining Sleeve 

Since the PM generator is spinning at high speed, the rotor and permanent magnets are 
subjected to extremely high centrifugal forces. These forces can cause significant amounts of 
damage if the magnets and rotor are not properly restrained. The rotational components can be 
strengthened by enclosing them in a retaining sleeve/can which also increases the air gap length, 
The centrifugal force on the magnets due to the rotor spinning is ealculated in Eqn 5-5. 

M -v ^ 
p _ m ^mag 

R + h^ 

where = mass of magnets 

''mag = velocity of magnets (m/s) 


Eqn 5-5 

Using the inner surface area of the retaining sleeve, this force is converted to an outward 
pressure. Treating the retaining sleeve as a thin-walled vessel, the hoop stress felt by the sleeve 
is determined as shown in Figure 30 and Eqn 5-6. 



Figure 30: Retaining Sleeve Hoop Stress 

rn 

y Fy = -2-aj, h L+ P rLsin(0) Lde = 0 


vert 


Pr 


^^0 = 


80 


Eqn 5-6 









To provide a margin for mechanical tolerances and imperfections, a safety factor (SF) of 1.2 is 
applied to the hoop stress calculated in Eqn 5-6 to get a final SF stress for the retaining sleeve. 

The retaining sleeve can be made from many different types of materials including metal 
alloys and composites. A disadvantage of a metallic sleeve is eddy currents are induced in the 
sleeve by variations in the flux density caused by the stator slots. The advantages of a metallic 
sleeve are that it shields the magnets from most of the flux density variations and it has a high 
thermal conductivity for heat removal. 

A composite, such as carbon-fiber, provides reasonable strength while having lower 
losses since the lower conductivity reduces the eddy currents. However, the composite sleeve 
has low thermal conductivity and does not shield the magnets from the flux variations. This 
results in increased losses in the magnets themselves. In general, when the sleeve conductivity is 
low, the rotor losses are in the permanent magnets, while when the sleeve conductivity is high, 
the PM losses decrease and the sleeve losses increase. The materials considered for use as in the 
retaining sleeve are listed in Table 14 [10], [57], [58], [59], [60], [61], [62], [63], [64]. The 
numbers for the carbon fiber composite are in the middle of a range of typical values (0.5 - 18.0 
pI2-m) since wide variations exist depending on how the composite is made. 


Table 14: Retaining Sleeve Materials 





Stainless Steel 

90 

0.72 

Aluminum Alloy 

75 

0.05 

Titanium Alloy 

110 

0.78 

Inconel 

132 

0.98 

Carbon Fiber 

100 

9.25 


0 contains detailed MATLAB code to perform the retaining sleeve stress calculations. 


The initial PM generator from section 3.4.2 is analyzed with the results below. 


Retaining Sleeve Stress: 

Stress Limits: 

Stainless Steel = 90.0 ksi 

Titanium Alloy = 110.0 ksi 

Inconel = 132.0 ksi 

Actual Sleeve Stress: 

Sleeve Thickness Actual Stress 
0.50 mm 1136.9 ksi 

1.00 mm 568.4 ksi 

1.50 mm 379.0 ksi 


Aluminum Alloy = 75.0 ksi 

Carbon Fiber = 100.0 ksi 


SF Stress 
1364.2 ksi 
682.1 ksi 
454.7 ksi 
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2.00 mm 

284.2 ksi 

341.1 ksi 

2.50 mm 

227.4 ksi 

272.8 ksi 

3.00 mm 

189.5 ksi 

227.4 ksi 

3.50 mm 

162.4 ksi 

194.9 ksi 

4.00 mm 

142.1 ksi 

170.5 ksi 

4.50 mm 

126.3 ksi 

151.6 ksi 

5.00 mm 

113.7 ksi 

136.4 ksi 

5.50 mm 

103.4 ksi 

124.0 ksi 

6.00 mm 

94.7 ksi 

113.7 ksi 

6.50 mm 

87.5 ksi 

104.9 ksi 

7.00 mm 

81.2 ksi 

97.4 ksi 

7.50 mm 

75.8 ksi 

90.9 ksi 

8.00 mm 

71.1 ksi 

85.3 ksi 


It is evident that the initial 4 mm gap length is not sufficient once the retaining sleeve is 
considered (since the air gap dimension includes the retaining sleeve thickness). The machine 
must therefore be revised to allow for the retaining can. In this modification process, the 
aluminum alloy is removed as a sleeve material possibility since it has lower strength and much 
higher conductivity than the other materials which together would produce a larger air gap and 
higher eddy current losses. 

The PM generator will be redesigned taking the retaining sleeve and hoop stress limits 
into consideration. The updated PM generator will be a bigger machine with a much larger air 
gap, greater magnet height, lower Bg, lower voltage, and higher current density. 

5.3 Rotor Losses 

5.3.1 Model for Time Harmonics & Winding Space Harmonics 

The permanent magnets used in high-speed generators are electrically conductive and 
therefore support eddy currents. The retaining sleeves are sometimes made from electrically 
conductive material that also can carry eddy currents. These eddy currents are primarily caused 
by fluctuations in the magnetic flux density produced by time and space harmonics of the 
winding currents. The currents produce losses which can potentially cause excessive heating or 
demagnetization of the permanent magnets. An anal)4ical model is developed using the winding 
and current harmonics and the surface impedance to estimate the rotor losses. 
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Figure 31 shows the PM generator geometry “flattened out” into rectilinear coordinates. 
This is an accurate representation provided the dimensions are on a radial scale that is much 
smaller than the radius of the machine so that curvature is not important [65]. 



Figure 31: General Magnet Loss Model 

The direction of rotation is in the positive x-direction, the radial direction is y, and the armature 
current flows in the axial dimension, z. In addition to the arrangement of Figure 31, the 
geometry shown in Figure 32 is also utilized. 



Figure 32: Layer of Material 

The following assumptions are made in developing the rotor analytical loss model: 

• Layers of material extend to ± infinity in the ± x direction 

• Layers effectively extend to negative infinity in the negative y direction 
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• Motion/rotation is in the + x direction 

• The physical constants of the layers are homogeneous, isentropic, and linear 

• The ferromagnetic material does not saturate 

• The machine is long axially so magnetic variations in the z direction are ignored (H and 
B only vary in x, y directions) 

• All currents flow in the z direction 

• The rotor and stator are constructed of laminated steel so their conductivity in the z 
direction is negligible 

• The time and space variations are approximately sinusoidal 

• Flux density at y = infinity is zero 

• A traveling flux wave harmonic can be represented by an equivalent traveling current 
sheet (Kz) on the surface of the stator 

• The normal component of the flux density is continuous at all interfaces 

• The tangential component of the flux density is continuous at all interfaces except at the 
stator/air gap where it is increased by the current sheet density 

• The magnetic flux density crossing the air gap and magnets is perpendicular 

• The effect of magnet eddy currents on the magnetic flux density is negligible - this is 
accurate below 10 kHz [66] 

• The magnet flux density is constant over the magnet breadth 

One of Maxwell’s Law states that there is no magnetic charge and therefore flux lines 
must close on themselves, represented by Eqn 5-7. 

VB = 0 

Eqn 5-7 

Knowing this, the flux density B can be represented as the divergence of the magnetic scalar 
potential A (Eqn 5-8). 

B=VxA 

Eqn 5-8 

Substituting in poH for B, taking the gradient of Eqn 5-8, and using Ampere’s Law for current 
density yields Eqn 5-9. 


84 






jiq H = V X a 


IXqV xH=VxVxA 


HoVxH=VlV Aj-V -A 

^ 2 
^q J = -V A 

2 '^ 

|XQ a E= -V A 

Eqn 5-9 

Faraday’s Law is then employed to determine a relationship between E and A as shown below. 

V X E = —B (Faraday's Law) 

dt 

V X E = -V X — A 

dt 

j —^ 

E = — A 


Combining Eqn 5-9 and Eqn 5-10 gives Eqn 5-11. 


Eqn 5-10 


2 ' H ' 

V -A = p Q-a-—A 
^ dt 

Eqn 5-11 

Since H and B vary only in the x and y directions, A and E are in the z direction and A can be 
represented by Eqn 5-12. 


A,= Re[(A,(y)).eH“-'“)] 

where k = 2 n/X 


Eqn 5-12 


Performing the mathematical operations and solving, Eqn 5-12 becomes Eqn 5-13. 
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2 d 

V *Ay — IXn’CJ- Ay 

Z ^ dt ^ 


< J 

Since b = V x A ^nd A is in z-dir, then = — A„ and b^ = — 


X . ""Z 

dy 


dx 


^2 _ 


^ ^ 
—■ + — 2 ' + 


dx dy dz y 


,2 d^ 

= -k + —Ti 

dy 


-k^-A^Cy) + -^A^Cy) = j-co-^o-<5-Az(y) 
dy 


-^A^Cy) - (k^ + j-(o-|io-cjj-Az(y) = 0 
dy 


Let Y = Jk^'+ jm^iQ-o 
Solution of the form: 

jy , A .-yy 


Az(y) = Ap-e +An-e 
{( 


A, = Re| I A„.e™ + A„.e-”)ei 


Solving for the magnetic flux densities and magnetic fields produces Eqn 5-15. 

yy^ ^j-(cot-kx) 
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Eqn 5-13 


Eqn 5-14 
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Y-Ap -y-Ajj 

Hp = — and H„- 




**0 


«x=*{(hp' 
"y = 


+ Hjj-e' 


■7y^.gj(“t-kx) 


Y 


- YY 1 J(®t-kx) 


= Rd Hx-e 


j(a)t-kx) 


Hy = Re| 


Hy-e 


j(a)t-kx) 


Eqn 5-15 


The layer from Figure 32 can now be examined utilizing a “surface coefficient.” A surface 
coefficient is defined as the ratio of the y-directed to x-directed magnetic field amplitude (a = 
Hy/Hx). At the bottom of the layer where y = 0, the surface coefficient is given by Eqn 5-16. 

“b(y = ») = — = s- k 


H 


X Y Hp + Hj, j._. 

Y 


Hn 


+ 1 


Eqn 5-16 


At the top of the layer (y = h), the surface coefficient is shown in Eqn 5-17. 




cCtCy = h) = j—• 
Y 




- ^ 


^yh 


Hp.e' +H„e 
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«! = J— 
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_p eYh_e-Yh 


—^-e' + e ' 




J 


Eqn 5-17 

Using Eqn 5-16 to solve for the ratio Hp/Hn, a final expression for the top surface coefficient is 
determined (Eqn 5-18), and it is applicable to any uniform region. 
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Eqn 5-18 

If the region being examined is positioned on top of a ferromagnetic surface, such as the magnets 
on the steel rotor shaft, the boundary condition at the bottom of the layer (Ub infinity as Hx ^ 
0) produces Eqn 5-19. 


a ^ = j—coth( 7 h) 
Y 


Eqn 5-19 


In the case of the air gap where the conductivity is zero, Eqn 5-18 reduces to Eqn 5-20. 
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at = j- 


f j-sinh(')li) + a|j-cosh(')li)^ 
j-cosh(')!i) + a^-sinh(')h) 


Eqn 5-20 

The electric field (E) is z-directed and the magnetic flux density and field (B and H) are 
y-directed so a relationship is determined between them using Faraday’s Law. The surface 
impedance is then determined. Surface impedance (Zs) is the ratio of the z-directed electric field 
to the z-directed current. Eqn 5-21 and Eqn 5-22 show the detailed calculations. 


V X E = -5-b 
dt 


(Faraday's Law) 


V X E = -Un —H 
^ dt 

dz dx dt ■’ 


jkE2 = -jtopo-Hy 

to„ 


K 


-Ho-H. 


Eqn 5-21 




Z. = 


®surf ^ ^ 


-l^oHy 


^surf 

“n 

Zs = 


-H. 


Eqn 5-22 

Since electromagnetic power flow into the rotor is the desired quantity, Poynting’s vector 
is employed (Eqn 5-23). 

S = Ex H 


Eqn 5-23 


Dissipation in the rotor is in the negative y direction producing Eqn 5-24. 
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There is no z-directed magnetic field (H^ = 0) 

-Sy = -%H^ 

In time average form, Eqn 5-24 becomes Eqn 5-25. 

-<Sy>=~Re(Ez-Hx) 

-<5,>=--|i//-Re(—) 


Eqn 5-24 


-ReCZ) 

Eqn 5-25 

Eqn 5-25 yields the power dissipated at the stator surface. This is the correct result for 
the rotor because there is no mechanism for dissipating power between the stator and rotor. The 
power estimated by Poynting’s theorem flows directly from the stator to the rotor [65]. 

In order to calculate the rotor losses, the above analytical model is applied to the 
geometry in Figure 31. For this model, the stator is assumed to be a smooth surface without slots 
because the slot effects are considered later in section 5.3.2. The first step is calculating the 
surface coefficient at the bottom of the magnet layer. It is assumed that this is formed by the 
highly permeable rotor shaft below the magnets. This assumption allows the surface coefficient 
at the top of the magnet layer to be calculated using Eqn 5-19. 

Next, the surface coefficient at the top of the retaining sleeve is determined using the 
sleeve material’s conductivity and Eqn 5-18. Traveling across the air gap to the surface of the 
stator, the surface coefficient is estimated using Eqn 5-20. The surface impedance is then 
computed utilizing Eqn 5-22. 

Once the surface impedance is known, the rotor losses caused by the stator winding time 
and space harmonics are calculated using Eqn 5-25. Only the more significant harmonics are 
considered as shown in Figure 33. The higher order harmonic effects are ignored because the 
products of higher order time harmonics (small number) and higher order space harmonics 
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(small number) produce negligible contributions (small number * small number = very small 
number). Appendix I contains the detailed MATLAB code that performs all of the above 
calculations. 



Figure 33: Relevant Harmonics 


5.3.2 Model for Stator Slot Effects 

The stator slots cause variations in the magnetic field which produce losses in the 
retaining sleeve and magnets of the rotor. Accurate calculation of the losses in the retaining 
sleeve is extremely difficult. Several different methods have been developed and in this paper, 
the technique from reference [18] is employed. Figure 34 shows the geometry of the retaining 
can and the currents that are induced. 



Figure 34: Retaining Sleeve Induced Currents 
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As the rotor spins past the slot openings of the stator, the air gap flux density undergoes 
modulation due to the change in reluctance. The dip in Bg (shown in Figure 35) travels along the 
B-waveform which is otherwise moving synchronously with the rotor. 



Figure 35: Flux Density Variation 

The rotation of the rotor generates an E-field in the sleeve and a subsequent axial current 
density J = E/p. When this current density is integrated over the volume of the can, the average 
loss per unit area is determined (Eqn 5-26) [18]. 

w = 

Eqn 5-26 

It is evident from Eqn 5-26 that as the slot width increases, the width of the flux density dip (P) 
gets larger causing the sleeve losses to increase. The above equation only considers eddy 
currents flowing in the axial direction but there are also circumferential components. These 
portions are accounted for using a factor Ks as shown in Eqn 5-27 where the total can losses are 
determined [18]. 


3600 p y 
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Pcan =KsWA 


where a = .r.2.(R + h J-L,, 


tanh 

PLst 





P-Lst 



_2(R + l'm)_ 



Eqn 5-27 


One way to reduce the retaining can losses is to split the sleeve cylinder into separate 
rings. This reduces the effective flow path length for the eddy currents thus lowering the losses. 
There is a limit to how much this can be done since as the number of rings increases, the 
construction becomes more difficult. In this paper it is assumed that the maximum feasible 
number of sleeve sections is ten. 

When metal alloys are used as the retaining sleeve material, they effectively shield the 
magnets from the flux density variations, whereas there is little or no shielding with the carbon 
fiber composite because of its low conductivity. Therefore, with a metal alloy can, the rotor 
losses occur largely in the can compared to the carbon fiber sleeve, where the losses mostly 
occur in the magnets. These magnet losses are calculated using methods similar to Eqn 5-26 and 
Eqn 5-27 assuming that the eddy currents flow in the top 10% of the magnet volume. Appendix 
J contains the MATLAB code that performs the detailed rotor loss calculations caused by the 
stator slot effects. 


5.4 Complete System Model & Design Procedure 

The complete model for the power generation module consists of the PM generator, the 
rectifier, the filtering/energy storage components, the buck converter, and the load. The program 
PSIM (Version 6.0, Powersim, Inc.) is utilized to construct the circuit model, perform detailed 
time series and waveform analysis, and determine voltage and current values. An example of a 
7-phase version of this model system is shown in Figure 36. In each of the models, the converter 
and rectifier are simulated as one IGBT/diode module for simplicity but for actual 
implementation/construction, 19 modules are required as described in section 4.5. 


93 







CESERWOR RECTiriER IMPOT FILTER C.OBVERTER OtraPVT iPJLTEi^OtJTRUT 

Figure 36: System Model 

The procedure for designing and analyzing the power generation module consists of the 
following: 

• Design a PM generator using the methodology from sections 3.4.2 and 5.2, iterating to 
ensure the sizes, current density, magnetic saturation, and sleeve thickness/stresses are all 
satisfactory 

• Verify the PM generator has proper output waveforms and THD per section 5.1 

• Calculate the rotor losses caused by stator slot harmonics as outlined in section 5.3.2 

• Determine the rectifier circuit parameters utilizing the procedures from section 4.2 

• Design the buck converter and associated input and output filters employing the methods 
of section 4.3 

• Construct the PSIM circuit using the correct generator parameters (number of phases, 
back voltage (Ea), and synchronous inductance (Ls)), rectifier values, converter 
parameters, and filtering components 

• Simulate and analyze the PSIM model verifying the system meets the required 
specifications for power and output voltage 

• Determine the generator line harmonic current magnitudes from the PSIM model 

• Using the harmonic current values, calculate the rotor losses caused by the winding time 
and space harmonics as detailed in section 5.3.1 

• Calculate the switching losses and conduction losses of the power conversion module per 
section 4.4 

• Determine the overall losses and efficiency of the system 
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• Calculate the masses and volumes of the generator and power conversion module using 

the methodology of sections 3.4.2 and 4.5 

5.5 Optimization 

Optimizing the overall power generation module involves designing numerous PM 
generators and associated power conversion electronics to achieve an ideal result. Using the 
design procedure from section 5.4, machine parameters are varied producing different generator 
designs. For each unique generator, a power conversion module is devised so that its values are 
matched with the output of the generator. This results in a wide range of power generation 
modules. 

The main machine parameters that are varied are the number of poles, number of phases, 
and the retaining sleeve material as shown in Figure 37. This results in 80 different machines 
and associated power conversion modules. For each of these generators, the magnet thickness 
(hm), air gap (g), magnet angle (Gm), and stack length (LsO are iterated to ensure the sizes, current 
density, magnetic saturation, sleeve thickness/stresses, output waveforms, and THD are all 
satisfactory. 

Rotational speed (13,000 RPM) and magnet skew angle (10°) are held constant. The 
number of slots (N) is set to 36 for the 3 and 5-phase machines and to 72 for the 7, 9, and 11- 
phase machines to properly fit the windings and ensure a reasonably sinusoidal back voltage 
waveform. Other parameters including the number of slots/pole/phase (m), number of armature 
turns (Na), back voltage (Ea), synchronous inductance (Ls), electrical frequency (f), breadth 
factor (kb), pitch factor (kp), skew factor (ks), magnet factor (kg), and air gap flux density (Bg) all 
change as a result of the machine variations. 
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Figure 37: Machine Optimization Parameters 

Each machine is optimized to have minimal back voltage, size, and mass while ensuring 
it remains within the current density and magnetic saturation limits. The power conversion 
modules are then developed matching component ratings to those of the PM generators. The 
power eleetronics are designed to produce maximum power factor and minimal line harmonics 
as discussed in section 4.3. This yields the lowest possible losses in the PM generator thus 
maximizing efficiency and minimizing rotor heat generation. 
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Chapter 6 Results and Analysis 


6.1 General 

As discussed in section 5.5, the main machine parameters that are varied are the number 
of poles, number of phases, and the retaining sleeve material and resulting in 80 different 
machines and associated power conversion modules. These PM machines and power 
conversions modules are developed with varying results. Appendix K and Appendix L contain 
the detailed specifications for each of the generators and power conversion modules. Appendix 
M and Appendix N have the in-depth losses and weights for each design. Appendix O lists the 
mass and volume results for the rectifiers and Appendix P contains the mass and volume data for 
the converters. A summary of the range of values for the designs is contained in Table 15 and 
Table 16. 


Table 15: General Module Specifications 



Back 

Voltage Eg 
(RMS) 

Synch 

Inductance 

Ls(mH) 

Gen Volt 
THD 

Stack 

Length 

l-8t(m) 

Overall 

Machine 

Length 

L(m) 

Machine 
Diameter D 
(m) 

Air Gap Flux 
Density 
Bg(T) 

PCM Vol 
(m®) 

Minimum 

1,136.0 

0.020 

3.12% 

0.850 

1.732 

0.834 

0.61 

6.54 

Maximum 

3,380.0 

0.348 

23.05% 

1.080 

2.057 

0.900 

0.77 

7.96 


1,977.0 

0.110 

11.61% 

0.943 

1.861 

0.866 

0.68 

7.14 


Table 16: General Module Results 



Machine 

Weight 


Total 

Weight 

(kg) 

Machine 

Losses 

(kW) 

■Bl 

hSH 

Line 

Current 

THD 

Minimum 

1,388.6 

2,302.7 

3,900.1 

241.7 

1,688.3 

2,126.7 

6.80% 

Maximum 

2,302.9 

3,954.7 

6,101.8 

1,259.8 

3,293.7 

3,535.4 

24.19% 


1,754.6 

3,025.5 

4,780.2 

463.8 

2,124.0 

2,587.8 

14.56% 


In general, the PCM losses are too high due to the IGBT switching losses. This is 
because high-frequency hard switching is used in the power conversion module as opposed to 
soft switching or a lower frequency. Soft switching involves using snubber circuits or resonant 
converters to minimize device switching losses by constraining the switching of the power 
devices to time intervals when the voltage across the device or the current through it is nearly 
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zero. Soft switching can reduce switching losses as much as 50-80% [67], [68]. In addition to 
lowering switching losses, soft switching reduces EMI, permits higher switching frequencies, 
and achieves higher efficiencies. Analyzing different soft-switching topologies is beyond the 
scope of this paper but the effects soft switching has on the PCM are studied. Using a lower 
switching frequency is not examined because this would produce a larger converter since the 
passive components would have to be bigger. 

In order to select a final design, the different variants are examined to ascertain which 
one is optimal. Parameters such as weights, machine losses, PCM losses, and THD are used to 
help facilitate the selection process. 

6.2 Number of Phases 

To facilitate studying the effects the number of phases has on the power generation 
module, several parameters are held constant. The retaining sleeve material and number of poles 
are arbitrarily set to inconel and 12. This does not affect the analysis because the trends related 
to varying the number of phases are consistent across all 80 power modules as indicated in the 
detailed results contained in the appendices. 

As the number of phases varies from 3 to 11, the weights of the PM generator and the 
power conversion module change, but only slightly. Figure 38 shows this and there is a 
maximum variation of approximately 242 kg, or 5.7%. The 11-phase machine produces the 
lowest weight power module but only by a small amount. 


Weights vs. Phase Number 
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Figure 38: Weights vs. Number of Phases 
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The machine losses fluctuate significantly as the phase number changes with the 3-phase 
machine having 48.8% of the losses of the 11 phase machine (see Figure 39). The core, 
windage, and armature losses are approximately constant for the different phase generators, but 
the rotor losses vary considerably. This disparity is caused predominately by the AC line current 
harmonics from the rectifier being much greater in the higher phase machines. Their triple-n 
harmonics are higher order thus allowing lower order harmonics to have a greater effect, 
specifically the 3^'* harmonic. 

The higher rotor losses are also caused by the geometry of the machine. As the phase 
number increases, the number of slots/pole/phase (m) decreases which causes the winding 
factors to have less of an effect in reducing higher harmonic effects. The exception is the jump 
from 5-phases to 7-phases because the number of slots is doubled from 36 to 72 so that m 
increases initially. Overall, the 3-phase generator has the lowest rotor and total machine losses. 


Machine Losses vs. Phase Number 
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Figure 39: Machine Losses vs. Number of Phases 
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PCM Losses vs. Phase Number 
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Figure 40: PCM Losses vs. Number of Phases 
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The PCM losses also differ greatly depending on the number of phases (Figure 40). The 
11-phase machine has the lowest total losses due to it having the lowest IGBT switching losses. 
This is because it has a smaller rectified voltage which reduces the stresses and switching losses 
in the IGBTs (see section 4.4). However, if soft switching is implemented, the PCM losses are 
similar for the different phase machines varying by a maximum of 63.5 kW instead of the 843 
kW associated with hard switching. Therefore, the lower phase machines become more 
attractive if soft switching is utilized. 

Lastly, as the number of phases increases. Figure 41 shows that in general, the generator 
voltage and line current THD increase. This is caused by the same effects as with the rotor 
losses, the slots/pole/phase and the AC line harmonics from the rectifier. As mentioned in 
section 5.1, the specification from MIL-STD-1399 for voltage THD is normally 5% so using 
these standards as guidelines, only the 3-phase and the 7-phase power modules are satisfactory. 
The line current THD should be made as low as possible because this directly affects the amount 
of losses and heating in the rotor. The typical limit for line current THD is 10% making the 3- 
phase generator the only one with acceptable THD performance [69]. 
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THD vs. Phase Number 
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Figure 41: THD vs. Number of Phases 


6.3 Retaining Material 

Similar to the analysis of the number of phases, several parameters are held constant 
when examining the various retaining sleeve materials. Specifically, the number of phases is 
arbitrarily set to 3 and the number of poles is fixed at 6. This does not affect the analysis 
because the trends related to varying the retaining material are consistent across the continuum of 
power modules as shown in the detailed results contained in the appendices. 

With respect to the weights of the PM generator and the PCM, the type of retaining 
sleeve has little effect. Using a stainless steel retaining can produces the largest power module 
but it is not appreciably larger than the others. When the retaining sleeve is constructed from 
titanium, inconel, or carbon fiber, the difference in total weight is less than 0.11%. Figure 42 
shows the effects the retaining sleeve material has on the weights of the power generation 
module. 
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Figure 42: Weights vs. Retaining Sleeve Material 



Figure 43: Machine Losses vs. Retaining Sleeve Material 

The type of retaining sleeve material does affect the losses in the PM generator. The 
core, windage, and armature losses vary only slightly but the rotor losses fluctuate a great deal. 
Inconel has the lowest rotor losses (77.3 kW) while stainless steel has the highest (107.1 kW). 
The carbon fiber has the smallest losses in the actual can but larger losses in the magnets. This is 
because the carbon fiber provides no shielding for the magnets from flux variations and winding 








current harmonics. Overall, the inconel sleeve produces the lowest rotor and total machine 
losses as shown in Figure 43. 

The PCM losses are not influenced considerably by the kind of retaining can material 
(Figure 44). Depending on the type, the losses vary by about 157.2 kW, or 4.2%. When soft 
switching is utilized, the losses fluctuate by only 23.7 kW or 2.5%. Therefore, the PCM losses 
are not a major factor in selecting the retaining can material. 


PCM Losses vs. Retaining Sleeve Material 
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Figure 44: PCM Losses vs. Retaining Sleeve Material 

Lastly, the selection of a retaining sleeve material is not a factor in the generator voltage 
THD or the AC line current THD. The voltage THD varies by only 0.5% and the line current by 
less than 0.15% depending on the t 5 rpe of material as indicated in Figure 45. It is evident, 
however, that selected a 3-phase, 6-pole, 36-slot machine would not be wise since the voltage 
THD is above the 5% limit in all cases. 
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Figure 45; THD vs. Retaining Sleeve Material 
6.4 Number of Poles 

The number of poles in a machine is an important factor because it affects numerous 
parameters including the electrical frequency, the magnet pole pitch, and the magnetic air gap 
flux density. To facilitate studying the effects of the number of poles, several parameters are 
held constant. The retaining sleeve material and number of phases are arbitrarily set to inconel 
and 3. This does not affect the analysis because the trends related to varying the number of poles 
are consistent across the range of machines as indicated in the detailed results in the appendices. 

First, the weights of the generator and PCM are examined to ascertain the effects of the 
number of poles. As Figure 46 indicates, as the number of poles increases the weights of both 
the machine and the PCM decrease. The generator weight goes down because assuming the 
same flux density and circumferential arc, a greater number of poles produces the same radial 
flux but requires less stator core back iron (see Figure 5). The PCM weight decreases because 
the higher frequencies and lower rectified voltage associated with the increased pole numbers 
reduce the size and energy requirements of the passive components, especially the capacitors. 
Therefore, the total power generation module weight decreases as the number of poles increases 
making the 12-pole machine the most attractive in terms of weight. 
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Weights vs. Pole Number 
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Figure 46: Weights vs. Number of Poles 

Next, looking at Figure 47, as the number of poles goes up, the losses in the machine 
increase, some more than others. The core losses increase because the higher frequencies cause 
the eddy currents and hysteresis to have a greater effect. However, the range in core losses is 
only 7 kW so this increase is not a major contributor. The windage losses are not significantly 
affected by the number of poles differing by only 3 kW across the pole variation. 

Machine Losses vs. Pole Number 
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Figure 47: Machine Losses vs. Number of Poles 

The armature losses do fluctuate some as the pole number changes with the 6-pole 
machine having the lowest (117.6 kW) and the 12-pole machine having the highest (137.9 kW). 
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This is due to the higher pole machines having higher stator current densities as the machine 
performance is optimized. The rotor losses also increase as the number of poles gets higher 
mainly because of the winding time and space harmonics. As the pole number increases, the 
number of slots/pole/phase (m) decreases which causes the winding factors to have less of an 
influence in reducing higher harmonic effects. 

Overall, all of the different pole machines have less than 160 kW of rotor losses with the 
6-pole machine having the lowest. However, as mentioned in section 6.3, the generator voltage 
THD is more significant in this machine causing the higher pole machines to become viable 
options. 

The number of poles also produces sizeable effects on the PCM losses. As the pole 
number increases, the PCM losses decrease due to much lower IGBT switching losses. This is 
because the higher pole machines have lower rectified voltages which reduce the stresses and 
switching losses in the IGBTs (see section 4.4). The difference between having 6-poles and 12- 
poles is 754.3 kW in PCM losses, a 22.9% reduction. If soft switching is implemented, the PCM 
losses are still lower for the higher pole machines, but the difference is now 98.1 kW, a 10.3% 
reduction. Overall, as shown in Figure 48, the 12-pole machine is the best option when 
considering PCM losses. 
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Figure 48: PCM Losses vs. Number of Poles 

The generator voltage THD and the AC line current THD are both affected by the 
number of poles. Examining Figure 49, the 10-pole and 12-pole generators are the only variants 
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that produce voltage waveforms that have acceptable THD levels, 3.74% and 4.95% 
respectively. The AC line current THD values are satisfactory for all of the pole numbers with 
the 12-pole machine having the lowest value (6.82%). If the THD is examined in conjunction 
with the machine losses, the 12-pole machine is the best option. It has the top overall THD 
performance, has 25.8 kW lower rotor losses compared to the 10-pole generator, and has 14.6 
kW lower total machine losses compared to the 10-pole design. 
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Figure 49: THD vs. Number of Poles 
6.5 Final Power Module 

6.5.1 PM Generator 

Based on the previous analysis and the goal of producing an optimal power generation 
module, a final design is selected. The sizing codes and analysis from sections 3.4.1,3.4.2, 

3.4.3, and 3.4.4 are also utilized in this process. The final PM generator is a 3-phase, 12-pole 
machine with an inconel retaining sleeve. The final design is not a fully optimized in terms of all 
parameters but represents the best selection using the methodology of this paper. Detailed 
characteristics for the PM generator are listed in Table 17. 
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Table 17: PM Generator Final Design Parameters 



Total Machine Diameter 

0.848 m 

Total Machine Length 

1.868 m 

Rotor Radius 

0.147 m 

Active Length 

1.000 m 

Slot Average Width 

16.973 mm 

Slot Height 

25.000 mm 

Back Iron Thickness 

17.150 mm 

Tooth Width 

15.882 mm 




Power Rating 

16,000 kW 

Speed 

ImMioMM 

Va(RMS) 

1,739 V 

Phase Current 

3,169.3 A 

Ea (RMS) 

2,380 V 

Armature Resistance 

0.00457 Q 

Synchronous Reactance 

0.377 Q 

Synchronous Inductance 

0.046 mH 


\WJE!SKStSfB!!!M 

Tip Speed 

200 m/s 

Efficiency 

0.979 

Power Factor 

0.960 

Phases 

3 

Frequency 

1,300 Hz 





Number of Slots 

36 

Number of Armature Turns 

12 

Fund Breadth Factor 

1.000 

Fund Pitch Factor 

0.866 

Tooth Flux Density 

1.30 T 

Back Iron Flux Density 

0.93 T 

Slots/pole/phase (m) 

1.00 




iffuSilRS ’ 


29.00 mm 

Magnet Angle 

20 deg 

Air Gap 

5.(X) mm 

Pole Paris 

6 

Magnet Remanence 

1.20 T 

Air Gap Bg 

0.69 T 

Magnet Factor 

0.961 

Skew Factor 

0.955 





Core Loss 

19.9 kW 

Armature Loss 

137.9 kW 

Windage Loss 

34.4 kW 

Rotor Loss 

158.9 kW 



•v"' *■ 


Core 

292.49 kg 

Shaft 

522.73 kg 

Magnet 

72.59 kg 

Armature 

82.48 kg 

Services 

145.54 kg 

Structure 

557.92 kg 

Total 

1,673.75 kg 




A diagram of the PM generator is shown in Figure 50 with the stator winding layout 
given in Table 18. The numbers in the winding layout table refer to the slot number locations of 
Figure 50. 
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Figure 50: Diagram of PM Generator Final Design 


Table 18: Winding Layout 
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6.5.2 Power Electronics Module 

In conjunction with selecting the final PM generator design, the associated power 
electronics module is determined. Figure 51 shows the complete power module, including 
individual component values. 
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Figure 51: Power Module Final Design Diagram 

Along with Figure 51, additional characteristics for the final power module are given in Table 19 
with detailed data contained in the appendices. 

Table 19: PCM Final Design Parameters 




IGBT Switching Losses 

2,103.6 kW 

IGBT Conduction Losses 

267.6 kW 

Efficiency 

0.841 

iSiHE 


IGBT Switching Losses 

420.7 kW 

IGBT Conduction Losses 

267.6 kW 

Efficiency 

0.946 


Diode Conduction Losses 


Total PCM Losses 


Diode Conduction Losses 


Total PCM Losses 


168.2 kW 


2,539.4 kW 


168.2 kW 


856.5 kW 




Capacitors 


IGBTs 


Services 


Total 



396.2 k 


22.8 k 


267.5 k 


2,318.4 k 





2.0 m 


2.0 m 


•iinensi 


Width 


Volume 




1.66 m 


6.6 



Max Rectified Voltage 


Duty Cycle 
(@16MW load) 


Duty Cycle 
(@100kWload) 


5,567 V 

Output Voltage 

700.0 V 

0.1622 

Output Voltage Ripple 
(@16MW load) 

0.09 V 

0.1277 

Output Voltage Ripple 
(@100 kW load) 

1.40 V 
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6.5.3 PerformanceAVaveforms 


As discussed in section 5.4, all 80 power modules are constructed and simulated using 
PSIM. In order to conserve space, only the simulation results for the final design are presented 
here. The module is simulated operating at the maximum load of 16 MW and at the mi nimum 
load of 100 kW. 

For the 16 MW loading. Figure 52 shows the generator voltage waveforms. Since the 
THD is low, the curves fairly closely resemble sine waves. 



Figure 52: PM Generator Voltage Waveforms 


Figure 53 through Figure 55 show plots for various parameters for the complete module loaded 
at 16 MW. The output voltage fluctuates approximately 0.09 V (0.0013%) which is within 
specification, the output current varies by less than 22 A (0.0096%), and the AC line currents 
have low harmonic content. 
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Figure 54: Output Current (16 MW) 
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Figure 55: AC Line Current (16 MW) 


Figure 56 through Figure 58 display the module performance at the minimum load of 100 
kW. The output voltage ripple is 1.4 V, higher than at full load but still not significant. The 
output current also has some variation (about 25 A) but this is acceptable at the lower power 
levels. The AC line current has much higher harmonic content because of notches that occur in 
the waveforms caused by the DC link capacitor charging and discharging at the low load. 
However, the magnitudes of the harmonics are insignificant since the magnitude of the 
fundamental is low and therefore the rotor losses and heating are not adversely affected. 
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Figure 56: Output Voltage (100 kW) 
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Figure 57: Output Current (100 kW) 
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Figure 58: AC Line Current (100 kW) 


6.6 Comparison 

The premise for designing and employing a high-speed PM generator and power 
conversion module is to offer a light weight, power dense alternative to existing wound rotor 
machines. In Table 1 in section 1.3.2, parameters for some typical wound rotor machines are 
presented. In order to compare them with the new design power module, these machines are 
scaled to 16 MW using the general relationship of Eqn 6-1. 

4 

Pwr°^VoP 


Eqn 6-1 


Table 20 then provides an evaluation of the different machines and power modules. 


Table 20: Comparison of Machines/Modules 


fedchine 

Phases 






Volume 


1 

3 

16 

3,600 

4.4 

3.7 

3.2 

52.0 

40.8 

2 

3 

16 

3,600 

5.0 

3.6 

3.1 

54.3 

34.3 

3 

3 

16 

3,600 

5.5 

3.2 

4.0 

69.0 

47.3 

4 

3 

16 

3,600 

4.6 

2.8 

3.7 

47.6 

41.0 

PM Gen 

3 

16 

13,000 

1.87 

0.85 

0.85 

1.4 

1.67 

PCM 

N/A 

16 

N/A 

2.0 

1.66 

2.0 

6.6 

2.32 

Total PM 

3 

16 

13,000 

3.87 

Varies 

Varies 

8.0 

3.99 
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Examining Table 20, the new design offers significant reductions in both weight and 
volume. Specifically, the new total power module is 11.6 - 16.8% the volume and 8.5 - 11.7% 
the weight of the other machines. This translates to about a 7x reduction in volume and a lOx 
reduction in weight. These reductions can provide flexibility to naval architects since power, 
weight, and volume are integral parts of the design and construction processes. 
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Chapter 7 Conclusions and Recommendations 


7.1 Design Lessons Learned 

Numerous lessons are learned from both the generator and power electronics design 
processes. The entire 16 MW power module is a complex system with many interrelated 
components that affect the overall performance. 

7.1.1 PM Generator 

The final generator is the product of examining multiple parameters to determine an 
optimal design. Various lessons are learned relating to different portions of the machine 
including the number of slots, the number of poles, the number of phases, the magnetic 
geometry, the windings and output voltage, and optimization methods. 

7.1.1.1 Number of Slots 

Initially, it was thought that the number of slots could be held constant over the range of 
machine designs. However, this is not possible because as the number of phases increases, 
additional slots are necessary to fit the greater number of windings. One of the most important 
parameters that can be used to help reduce the harmonics in a rotating machine is the number of 
slots/pole/phase (m). Therefore, the number of slots can be used to help optimize this parameter 
and improve the shape of the output voltage waveform and THD. 

The number of slots also influences the electrical loading of the machine. The number of 
slots directly affects the slot size which helps determine the current density in the machine. This 
is usually one of the limits on the machine’s performance because of cooling. Overall, the 
number of slots is an important parameter that must be carefully selecting when designing a 
machine. 

7.1.1.2 Number of Poles 

Traditionally, the number of poles is driven by the system frequency and prime mover 
speed. With a high-speed generator rectified to a DC distribution system, the number of poles 
becomes an important design variable. Therefore, throughout the design progression, careful 
consideration should be given to the number of poles. This proved to be true throughout this 
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paper since the number of poles affects the electrical, magnetic, and structural performance, 
including the electrical frequency, the voltage waveforms, the magnetic flux, the magnet volume, 
the magnet pole pitch, the air gap size, and the stator back iron thickness. 

It is recommended that initially the designer select a range of acceptable output electrical 
frequencies which puts bounds on the allowable number of poles. From there, the number of 
poles can be varied to help determine the size, waveforms, and performance of the generator. In 
general, if the rotational speed is held constant as the number of poles increases: 

• The number of slots/pole/phase (m) changes affecting the output waveforms and THD 

• The weight of the machine decreases 

• The electrical frequency increases 

7.1.1.3 Number of Phases 

An odd number of phases is usually the best option for rotating machines because it 
produces balanced operation with limited mechanical fluctuations. Initially, it was assumed that 
a higher number of phases (greater than 3) would produce a more optimal power module. It is 
discovered that this is not necessarily the case. As the number of phases increases, m decreases 
which can adversely affect the voltage THD. Also, the AC line current harmonics are more 
substantial in increased phase machines because the triple-n harmonics are higher order as the 
phase number increases. However, a big advantage of having a higher number of phases is that 
it reduces the requirements for the PCM because of lower rectifier voltage ripple and lower 
phase voltages and currents. 

It is concluded that the benefits of 3-phases and a higher number of phases can be 
simultaneously realized through a modified design. Multiple 3-phase sets of windings can be 
located equidistant around the stator and then connected to the PCM to achieve the desired 
results. An example of this would be three 3-phase sets of windings, one set at oP, 120°, and 
240°; one at 45°, 165°, and 285°; and one at 90°, 210°, and 330°. 

7.1.1.4 Magnetic Geometry 

Designing the magnetic geometry in a PM generator is an intricate process. Careful 
consideration should be given to the type of retaining material, the magnet material, and the 
desired voltage. It is determined that the magnet height and air gap length greatly affect the 
overall machine performance, and if they are not properly selected, the generator does not 
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operated properly. A typical starting point in a design is to have the magnet height be 
approximately 5 - lOx the length of the air gap [17], [18], [30]. This ensures that the magnetic 
flux is adequate while not saturating the core material. 

7.1.1.5 Windings and Output Voltage 

As discussed in section 3.2.7, wye series connected windings are used in this paper. This 
produces a generator output voltage that is too high compared to the desired PCM output voltage 
of 700 V. It also unnecessarily stresses the PCM devices and produces a lower than optimal duty 
cycle. It is concluded that connecting some of the generator windings in parallel or using a lower 
number of turns to produce a lower output voltage might be viable options to alleviate this 
problem. However, care must be taken to account for possible circulating currents in the parallel 
windings in the event of a fault. 

7.1.1.6 Optimization Methods 

When designing the PM generator, the goal was to produce an optimal design that met 
the required performance criteria while minimizing the rotor losses, the size, and the weight of 
the machine. Optimizing the design involves trying to simultaneously maximize both the 
electrical loading (current density) and magnetic loading (air gap flux density). It is concluded 
that it is extremely difficult to do both at the same time since there are competing factors such as 
saturation limits, cooling limits, and rotational stress limits. 

In this paper, each of the 80 machines is optimized by trying to maximize the current 
density (within the limits) while ensuring an adequate air gap flux density. This is done because 
it seems to produce the smallest, high-performance machines. It is therefore recommended in 
any PM generator design process that careful consideration be given to both the electrical and 
magnetic loading. 

7.1.2 Power Electronics 

The power electronics module (PEM) designed in conjunction with the PM generator 
uses one of many possible topologies. The combination of the passive rectifier and buck 
converter is a stable, well-proven design. It is optimized to match the generator by selecting the 
passive components properly. However, the design does have limitations because it uses hard- 
switching and large passive components. 
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7.1.2.1 Rectification 

Passive rectification is utilized in the PEM because of its simplicity and potential for 
lower input harmonics. It is concluded that the rectification scheme is satisfactory and couples 
well with the PM generator. However, it is recognized that a more detailed analysis of the 
rectification scheme could produce a better topology. 

7.1.2.2 Conversion and Switching 

As discussed in Chapter 6, one of the main problems with the power electronics module 
is that the losses are too high due to the IGBT switching losses. This is because hard switching 
is used in the power conversion module as opposed to soft switching. Soft switching minimizes 
device switching losses, reduces EMI, permits higher switching frequencies, and achieves higher 
efficiencies. 

Analyzing different soft-switching topologies is beyond the scope of this paper but the 
effects soft switching has on the PCM are studied. It is concluded that if soft switching is 
implemented in the PCM, the efficiency increases by 10 - 20% with the losses reduced by 
overall 1 MW. Given the power requirements of the power generation module, soft switching 
should be employed. It is therefore recommended that in any high-power system that utilizes 
power conversion modules, soft switching should be examined to ascertain its usefulness. 

7.1.2.3 Control 

One control scheme for providing voltage regulation for the power generation module is 
discussed in section 4.3.4. This control methodology utilizing duty ratio control with a PCD 
controller is relatively simple to implement. This ease of control is possible because the buck 
converter is a straightforward topology. It is concluded that if a more complex, soft-switching 
design is implemented, more extensive control schemes would be required. 

Even though it is not discussed in detail, the PM generator provides for easy control 
because it is a synchronous machine. This is a distinct advantage over an induction generation 
which requires extremely complicated control algorithms. 

7.2 Power Generation Module 

Overall, 80 different 16 MW machines and power conversion modules are sized, 
designed, and analyzed with a final design selected. The permanent magnet generation module. 
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consisting of the generator and associated power electronics, provides an excellent alternative to 
traditional wound rotor machines for naval applications. It is concluded that the new design 
offers significant reductions in both weight and volume. Specifically, it is estimated that the PM 
generation module has a 7x reduction in volume and a lOx reduction in weight compared to 
similarly rated wound rotor systems. These reductions can provide flexibility to naval architects 
since power, weight, and volume are integral parts of the design and construction processes. 
However, further study is necessary to verify the 16 MW PM generation module’s thermal, 
structural, and mechanical performance. 

7.3 Recommendations/Further Study 

This paper discussed the electrical and magnetic design and analysis of a permanent 
magnet generation module including sizing, detailed machine design and analysis, analytical 
models for different rotor loss mechanisms as well as with other machine losses, power 
electronics and conversion for connecting the high-speed generator to a DC distribution system, 
and in depth simulation of the complete system. Numerous design issues are addressed and 
several issues are raised about the potential improvements a PM generation system can offer. A 
proposed PM generation module design is presented along with a detailed design methodology. 

However, this paper does not full cover all aspects associated with the design and 
analysis of a PM generation module. In order to fully ascertain whether the proposed design is 
technically sound, further study is required in both the machine and power electronics areas. 
Specifically, with the respect to the generator, the following is recommended: 

• Using the calculated losses, perform a detailed thermal analysis of the machine including 
development of the required cooling systems 

• Conduct an in-depth structural and mechanical study of the PM generator, including 
vibration and acoustic analysis 

For the power electronics portions, the further examination is recommended, including: 

• Examine different rectification schemes to see what effects they would have on the 
whole system 

• Explore soft switching topologies and perform detailed analysis on how to implement 
them to ensure acceptable switching losses while matching performance with the 
generator 
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For the complete PM generation system, the following is recommended: 

• Using the developed PM generator and power electronics module, implement a control 
scheme and conduct in-depth analysis with the system connected to an entire distribution 
system under changing load conditions 
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Winding Pitch 

Bottom Surface Coefficient 

Magnet Pitch Coverage Coefficient 

Top Surface Coefficient 

Flux Density Dip Width 

Flux Density Exponent 

Frequency Exponent 

Winding Electrical Angle 

Flux Linkage or Wavelength (depending) 

Flux Density Dip Interval 

Slot Fill Fraction 

Efficiency 

Kinematic Viscosity of Air at 20° C (m^/s) 

Power Factor Angle 

Resistivity 

Air Density at 20° C (kg/m^) 

Steel Density (kg/m^) 

Magnet Density (kg/m^) 

Conductor Density (kg/m^) 

Stator Conductivity (S/m) 

Magnet physical angle 

Magnet skew angle 

Air Gap Shear Stress 

Stator Tooth Width + Stator Slot Width 

Free Space Permeability 

Recoil Permeability 

Eleetrical Frequency (rad/sec) 

Mechanical Frequency (rad/sec) 



Magnetic Vector Potential 
Armature Conductor Area 
Air Gap Area 
Magnet Area 
Slot Cross-sectional Area 


Bo 

Bi 

Bb 

Bd 


Br 

B. 


Base Flux Density for Core Losses 
Fundamental Flux Density 
Back Iron Flux Density 
Flux Density Dip 
Ideal Radial Flux Density 
Air Gap Flux Density 
Remnant Flux Density 
Tooth Flux Density 
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specific Heat Capacity of Air (J/kg*C) 
Flux Concentration Factor 
Friction Coefficient 

Buck Controller Average Duty Cycle 
Stator Core Back Iron Thickness 
Duty Ratio 

Overall Machine Diameter 

Buck Controller Error Signal 
Back Voltage (EMF) 

Electrical Frequency (Hz) 

Base Frequency for Core Losses 

Air Gap 

Effective Air Gap 

Slot Depression Depth 
Magnet Height/Thickness 
Slot Depth 

Phase Current 

Average Converter Output Current 
Current Density 
Wave Number 

Winding Breadth/Distribution Factor 

Electrical Power Waveform Factor 

Magnet Factor 

Current Waveform Factor 

Winding Pitch Factor 

Winding Skew Factor 

Winding Factor 

Carter Coefficient 

Leakage Factor 

Reluctance Factor 

Retaining Can Loss Factor 

Surface Current Density 

Overall Machine Length 
Armature Conductor Length 
Air Gap Inductance 
Self Slot Leakage Inductance 



Lmach 


Lslot 

Ls 

Lst 


Mutual Slot Leakage Inductance 
End Turn Inductance 
Overall Machine Length 
Slot Leakage Inductance 
Synchronous Inductance 
Rotor Stack Length 


m 

Mac 

Me 

Mcb 

Met 

M^ 

Ms 

Mtot 


Slots per Pole per Phase 
Mass of Armature Conductor 
Mass of Core 
Mass of Core Back Iron 
Mass of Core Teeth 
Mass of Magnets 
Mass of Shaft 
Total Mass 


N 

Na 

Nc 

Ns 

Nset 

Nsfp 

Nsp 


Rotational Speed (RPM) 

Number of Armature Turns 
Turns per Coil 
Number of Stator Slots 
Actual Coil Throw (in slots) 
Full-pitch Coil Throw (in slots) 
Number of Stator Slots Short Pitched 


P 

Pf 

Po 

Pc 

Peb 


Pin 


Pwind 


PC 

PCM 

PEM 


Perm 

Pwr 


Pole Pairs 
Power Factor 

Base Power for Core Losses 
Total Core Loss 
Core Iron Loss 
Retaining Can Loss 
Tooth Loss 
Input Power 

Friction & Windage Losses 
Permeance Coefficient 
Power Conversion Module 
Power Electronics Module 
Permeance 
Machine Power 


q 


Number of phases 


R 

Ri 

R2 

Ra 


Rotor Radius 
Inner Bound of Magnet 
Outer Bound of Magnet 
Armature Resistance 
Core Inside Radius 
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Rco 

Rd 

Rds 

Ri 

Rs 

Rey 

syrat 

t 

tfrac 

THD 

Vmag 

Vtip 

Va 


Wd 

Ws 

Wsb 

Wst 

Wt 

Xs 

Zs 


Core Outside Radius 
Diode On-State Resistance 
IGBT On-State Resistance 
Inner Magnetic Boundary 
Outer Magnetic Boundary 
Reynold’s Coefficient 

Stator Back Iron Ratio (Yoke Thickness/Rotor Radius) 

Retaining Sleeve Thickness 
Peripheral Tooth Fraction 
Total Harmonic Distortion 

Magnet Speed (m/s) 

Rotor Tip Speed (m/s) 

Terminal Voltage 

Retaining Can Loss per Unit Area 

Slot Depression Width 

Slot Average Width 

Slot Bottom Width 

Slot Top Width (closest to rotor) 

Tooth Width 

Synchronous Reactance 

Surface impedance 
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Appendix A. Detailed Power Requirements 


In order to properly size the PM generator being designed, a typical load list is developed 
for an IPS naval ship. For the propulsion portion, it is assumed that the ship contains a single 
shaft requiring 30,000 SHP. Shaft horsepower (SHP) is the power measured in the shafting 
within the ship by a torsion meter as close as possible to the propeller or stem tube [10]. This 
SHP then translates to 32,552.1 BHP (24,275 kW) assuming a 96% efficient electric propulsion 
motor and a 96% efficient power control module for the motor. Brake horsepower (BHP) is the 
power required at the engines/generators. It is also assumed that the ship contains two equal¬ 
sized generators. Overall, Table 21 lists the systems and their associated connected (maximum) 
electrical power requirements. 


Table 21: Ship Connected Loads 





III Engineering Systems || 

Hydraulics 

150 

Engine Room Fresh Water 

50 

Main Sea Water 

300 

Electronics Fresh Water 

25 

Aux Sea Water 

80 

Reactor Cooling 

2,000 

Air Conditioning 

300 

Reactor Fresh Water 

50 

Chill Water 

50 

Reactor Heating 

350 

Distilling 

20 

Steam Feed System 

300 

Lube Oil 

50 

Steam Condensing System 

80 

Miscellaneous 

25 



II Auxiliary Systems III 

Atmosphere Control 

150 

Drain System 

100 

Refrigeration 

10 

Ventilation 

250 

High Pressure Air 

100 

Damage Control 

20 

Low Pressure Air 

25 

Sanitary System 

10 

Trim System 

100 



III Combat Systems III 

jlSonar 

250 

Radio/ESM 

40 II 

IlCombat Control 

250 

Navigation 

20 II 

III Hotel 

Loads II 

IlLighting 

60 


50 

IlLaundry 

20 

Galley 

70 



III Propulsion ||| 

Main Motor 

24,275 



Auxiliary Propulsion 

300 
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Classic U.S. Navy design practices (DDS 310-1) determine generator sizes by first 
estimating the expected load during the preliminary design phases. Propulsion load is a function 
of the hull form, maximum speed, and drive technology. Ship service loads are a function of the 
ship’s mission, crew size, and ship’s size. Once this is complete, the estimates are adjusted for 
anticipated construction and life cycle growth. Propulsion load is not expected to change but 
ship service loads are given a 20% margin for construction growth and a 20% margin for life 
cycle growth [11]. The generators are then sized using Eqn A-1, where Pgen is the power 
provided by one generator. 

_ (1.21.2 MaxLoad ) 

S®" ” [0.9 (n - 1)] 

where n = number of generators 

Eqn A-1 

Applying this method to the total electrical load from Table 21 produces two generators that are 
each capable of 47.9 MW. 

Since IPS ships require electrical power for ship loads and propulsion, efficient power 
management should be implemented. This means that power'is shared among all systems so that 
at top speed, more power is diverted to the propulsion motor, and at lower speeds power is 
available for other uses. It is assumed that rated ship service power and full propulsion power 
are possible at the same time only if both PM generators are available. If one generator is lost, 
the rated ship service load is still provided with the remaining power sent to the propulsion 
motor. In addition, full power for both the main motor and auxiliary propulsion is not required 
simultaneously. Similar to the classical method, a 20% margin for construction growth and a 
20% margin for life cycle growth are allotted. Based on these ideas, Eqn A- 2 is developed to 
size the generators, where Ptotai is the total generation capacity (both generators). 

^ (1.21.2-MaxNonPropulsionLoad) + MainMotorLoad 

” otai - ;;t1 

where n = number of generators 

Eqn A- 2 

Utilizing this process produces a total generation capacity of 31.99 MW or two generators that 
are each capable of 16 MW. 


136 


A standard U.S. Navy ship spends 95% of its time at less than full speed [12]. Because 
of this, the classical method unnecessarily over-estimates the power requirements for the ship 
and results in generator sets that are too large. If the average demand for each system is also 
considered, then the disparity becomes even worse. Table 22 lists each system along with a 
demand factor which represents the average amount of power required for normal operation of 
the systems. Based on the 9.8 MW total average loads, the classical method generators would 
each be loaded to only 10.2 % capacity which is extremely inefficient. 

The improved generator sizing method does a better job of ensuring efficient operation 
while not excessively limiting the ship if one generator is damaged or down for maintenance. If 
only one of the 16 MW generators is available, the connected ship service load (5.4 MW) can be 
supplied while still having 10.6 MW for propulsion, giving 44 % propulsion capability which is 
acceptable when operating in reduced status. 
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Table 22: Ship Average Loads 





i,A>v'<f,^^'ii'».>>^^:.r!'i!r.«.;w.fco Ass'iv'siwSi .s.> 




Engineering Systems 



Hydraulics 


Main Sea Water 


Aux Sea Water 


Air Conditionin 


Chill Water 


Distillin 


Lube Oil 


Engine Room Fresh Water 


Electronics Fresh Water 


Reactor Coolin 


Reactor Fresh Water 


Reactor Heatin 


Steam Feed System 


Steam Condensing System 


Miscellaneous 


Atmosphere Control 


Refrigeration 


High Pressure Air 


Low Pressure Air 


Trim System 


Dram System 


Ventilation 


Damage Control 


Sanitary System 


Sonar 


Combat Control 


Radio/ESM 


Navigation 


150 

1.00 

150 

300 

0.20 

60 

80 

0.20 

16 

300 

0.75 

225 

50 

1.00 

50 

20 

0.50 

10 

50 

1.00 

50 

50 

1.00 

50 

25 

1.00 

25 

2,000 

0.25 

500 

50 

1.00 

50 

350 

0.25 

88 

300 

0.75 

225 

80 

1.00 

80 

25 

1.00 

25 


Lightin 


Laund 


Water Heatin 


Galley 


Auxiliary Systems 


150 


10 


100 


25 


100 


100 


250 


20 


10 


Combat Systems 


250 


250 


40 


20 


Hotel Loads 


60 





II Propulsion I 

IlMain Motor 

24,275 

0.30 

7,283 

1 Auxiliary Propulsion 

300 

0.05 

15 
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Appendix B. MATLAB Code: Basic Sizing Method 


Jonathan Rucker, MIT Thesis 
%May 2005 
% Program: pmlbasic 

% Program performs basic sizing and parameter calculations 
% for generators. 


% Definition of variables 
% Name 

% General variables 
% Pwr 
% rpm 
% psi 
%f 


% omega 
% vtip 
% LovD 
% stress 

% Rotor variables 


Variable 


Required power 
Speed (RPM) 

Power factor angle 
Electrical fi*equency (Hz) 
Electrical frequency (rad/sec) 
Tip speed (m/s) 

L/D ratio 

Gap shear stress (psi) 


%R 

%D 

%Lst 

%p 

%Bg 

% Stator variables 
%Kz 
% Jz 
% hs 


Rotor radius (m) 

Rotor diameter (m) 

Rotor stack length (m) 

Number of pole pairs 
Expected air gap flux density (T) 

Surface current density (A/m) 
Current density (A/m2) 

Slot height (m) 


%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 


clear; 

% Constants & conversion factors 


hs = .015; 
lams = 0.5; 
convl = 9.81; 
conv2 = 703.0696; 


) Assume slot depth of 15 mm 
) Assume slot fill fraction 
) 9.81 W per Nm/s 
> 703.0696 N/m2 per psi 


%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 

% INPUTS 

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 

Pwr = 16e6; % Required power 

vtip = 200; % Max tip speed (m/s) 

LovD = 2.851; % Wound rotor usually 0.5-1.0, PM 1.0-3.0 ^ 

% Shear stress usually 1-10 psi small machines, 10-20 large liquid 
% liquid cooled machines 
stress =15; 

p = 3; % Pole pairs 

Bg = 0.8; % Tesla 


%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 

% Calculations 
% Size 
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% Initially use Pwr = 2*pi*R*Lst*stress*vtip 
% Lst = 2*LovD*R 
hscm = hs*100; 

R = sqrt(Pwr/(2*pi*(LovD*2)*vtip*stress*convl*conv2)); 

D = 2*R; 

Lst = LovD*D; 

% Speed 

omega = (p*vtip)/R; 
f = omeg^(2*pi); 
rpm = (60*f)/p; 

% Current densities 

Kz = (stress*conv2)/(Bg*100); 

Ja = 10*Kz/(hscm*lams); 

% Output 

fprintfCBasic Machine DesignXn'); 
fprintf('Input Parameters:\n'); 

fprintf('Power = %10.1f kW Shear Stress = %10.1f psi\n',Pwr/le3,stress); 
]^rintf(L/D Ratio =%10.2f Tip Speed = n[i/s\n',LovD,vtip); 

fi)rintf(’Pole Pairs = %10.1f Air Gap Bg = T\n',p,Bg); 

fprintfCOutputAn'); 

fprintf(’Rotor Radius = %10.3f m Stack Length = %10.3f m\n',R,Lst); 
fprintfCSpeed =%10.0fRPM Frequency =% 10.1 fHzW,rpm,f); 
fprintfCJa = %10.2f A/cm2\n’,Ja); 
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Appendix C. PM Machine Database 




Power factor 


Rotational Speed (rom 


Rotor radius 




Maanet thickness fmm 


Active lenoth (m 


Pole pairs 


Number of phases 


Number of Slots 


Slots per pole per phase 


Slots short pitched 


Slot depth (mm 


Slot depression depth (mm 


Slot depression width 


heral Tooth Fract 


Back iron thickness 


Turns per coil 




Slot fill fraction 


Maqnet remant flux density (T 


Maqnet electrical anqie 


Base frequenc 
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Rotor radius (cm 





Number of phases 


Number of Slots 


Slots per pole per phase 


Slots short pitched 


Slot depth (mm 


Slot depression depth (mm 


Back iron thickness (mm 


Turns per coll 


Coils per phase 


Slot fill fraction 


Maqnet remant flux density (T 
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Appendix D. MATLAB Code: Sizing Method 1 


% Jonathan Rucker, MIT Thesis 
% May 2005 
% Program: pm 1 input 
% Program used as input file for pmlcalc 
% All necessary input parameters entered here. 


clear; 

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 

% Definition & Entry of variables 
% General variables 

Pwr = 16e6; % Required power (W) 

rpm = 13000; % Speed (RPM) 

psi = 0; % Power factor angle 

% Rotor variables 


R = 0.147; 
hm = 0.025; 

Lst = 0.838; 
p = 3; 

Br= 1.2; 
thm = 50; 
thsk =10; 

% Stator variables 
q = 3; 

Ns = 36; 

Nsp= 1; 
g = .004; 
tfrac = 0.5; 
hs = .025; 
hd = .0005; 
wd = le-6; 
syrat = 0.7; 

Nc = l; 
lams = 0.5; 
sigst = 6.0e+7; 

% Densities 
rhos = 7700; 
rhom = 7400; 
rhoc = 8900; 


% Rotor radius (m) 

% Magnet thickness (m) 

% Rotor stack length (m) 

% Number of pole pairs 
% Magnet remnant flux density (T) 

% Magnet physical angle (deg) 

% Magnet skew angle (actual deg) 

% Number of phases 
% Number of slots 
% Number of slots short pitched 
% Air gap (m) 

% Peripheral tooth fraction 
% Slot depth (m) 

% Slot depression depth (m) 

% Slot depression width (m) 

% Stator back iron ratio (yoke thick/rotor radius) 
% Turns per coil 
% Slot fill fraction 
% Stator winding conductivity 

% Steel density (kg/m3) 

% Magnet density (kg/m3) 

% Conductor density (kg/m3) 


% Jonathan Rucker, MIT Thesis 
% May 2005 
% Program: pmlcalc 

% Program performs sizing and parameter calculations 
% for permanent magnet machines with surface magnets and 
% slotted stators. 

% Program developed from J.L. Kirtley script with permission 
% MUST RUN pmlinput PRIOR TO RUNNING pmlcalc 
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% Definition of variables 
% Name 

% General variables 
% Pwr 
% rpm 
% psi 
%f 


% omega 
% vtip 
% lambda 
%Ea 

% Rotor variables 


Variable 


Required power (W) 

Speed (RPM) 

Power factor angle 
Electrical frequency (Hz) 
Electrical frequency (rad/sec) 
Tip speed (m/s) 

Flux linkage 

RMS Internal voltage (V) 


% R Rotor radius (m) 

% hm Magnet thickness (m) 

% Lst Rotor stack length (m) 

% p Number of pole pairs 

% Br Magnet remnant flux density (T) 

% thm Magnet physical angle (deg) 

% thsk Magnet skew angle (actual deg) 

% Stator variables 


% q 
% m 
%Ns 
% Nsp 

%g 
% ge 
% tfrac 
%hs 
%hd 
% wd 
% syrat 
%Nc 
% lams 
% sigst 
%Kc 

% Loss Models 
% PO Base power for core losses 

% FO Base frequency for core loss 

% BO Base flux density 

% epsb Flux density exponent 

% epsf Frequency exponent 

% rhos Steel density 

% rhom Magnet density 

% rhoc Conductor density 

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 


Number of phases 

Slots per pole per phase 

Number of slots 

Number of slots short pitched 

Air gap (m) 

Effective air gap (m) 

Peripheral tooth fraction 
Slot depth (m) 

Slot depression depth (m) 

Slot depression width (m) 

Stator back iron ratio (yoke thick/rotor radius) 

Turns per coil 

Slot fill fraction 

Stator conductivity 

Carter coefficient 


% Constants to be used 
mu0 = 4*pi*le--7; 
tol = le-2; 
cpair= 1005.7; 
rhoair = L205 ; 
nuair= l.Se-S; 

P0 = 36.79; 

F0=1000; 

BO =1.0; 


% Free space permeability 
% Tolerance factor 

% Specific heat capacity of air (J/kg*C) 

% Density of air at 20 C (kg/m3) 

% Kinematic viscosity of air at 20 C (m2/s) 
% Base Power Losss, W/lb 
% Base freuency, 60 Hz 
% Base flux density, 1.0 T 


144 







epsb = 2.12; 
epsf=1.68; 

% Generate geometry of machine 

% Number of slots/pole/phase 
m = Ns/(2*p*q); 

% Number of armature turns (each slot has 2 half coils) 

Na = 2*p*m*Nc; 

% Tooth width 

wt= 2*pi*(R+g+hm+hd)*tfrac/Ns; 

% Slot top width (at air gap) 

wst = 2*pi*(R+g+hm+hd)*(l-tfrac)/Ns; 

% Slot bottom width 

wsb = wst*(R+g+hd+hs)/(R+g+hm+hd); 

% Stator core back iron depth (as p increases, dc decreases) 
dc = syrat*R/p; 

% Full-pitch coil throw 
Nsfp = floor(Ns/(2*p)); 

% Actual coil throw 
Nsct = Nsfp - Nsp; 

% Estimate end turn length 
% End turn travel (one end) 
laz = pi*(R+g+hm+hd+0.5*hs)*Nsct/Ns; 

% End length (half coil) 
le2 = pi*laz; 

% End length (axial direction) 
lei = 2*le2/(2*pi); 

% Calculate electrical frequency & surface speed 
f = p*rpm/60; 
omega = 2*pi*f; 
vtip = R*omega/p; 

% Winding & skew factors 

gama = 2*pi*p/Ns; 

alfa = pi*Nsct/Nsfp; 

kp = sin(pi/2)*sin(alfa/2); 

kb = sin(m*gama/2)/(m*sin(gama/2)); 

kw = kp*kb; 

ths = ((p*thsk)+le-6)*(pi/180); % skew angle (elec rad) 

ks = sin(ths/2)/(ths/2); 

% Calculate magnetic gap factor 
Rs = R+hm+g; 

Ri = R; 

R1 = R; 

R2 = R+hm; 

kg = ((Ri^(p-l))/(Rs^(2*p)-Ri^(2*p)))*((p/(p+l))*(R2^(p+l)-Rl^(p+l))... 
+(p*Rs^(2*p)/(p-l))*(Rl^(l-p)-R2^(l-p))); 

% Calculate air gap magnetic flux density 

% Account for slots, reluctance, and leakage 
ws = (wst+wsb)/2; % Average slot width 

taus = ws + wt; % Width of slot and tooth 
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Kc = l/(l-(l/((taus/ws)*((5*g/ws)+l)))); 
ge = Kc*g; 

Cphi = (p*thm)/180; % Flux concentration factor 

K1 = 0.95; % Leakage factor 

Kr = 1.05; % Reluctance factor 

murec =1.05; % Recoil permeability 

PC = hm/(ge*Cphi); % Permeance coefficient 

Bg = ((Kl*Cphi)/(l+(Kr*murec/PC)))*Br; 

% Calculate magnetic flux and internal voltage 
thmrad = thm*(pi/180); 

B1 = (4/pi)*Bg*kg*sin(p*thmrad/2); 
lambda = 2*Rs*Lst*Na*kw*ks*B 1/p; 

Ea = omega*lambda/sqrt(2); % RMS back voltage 

% Calculation of inductances/reactances 

% Air-gap inductance 

Lag = (q/2)*(4/pi)*(muO*Na'^2*kw'^2*Lst*Rs)/(p^2*(g+hm)); 

% Slot leakage inductance 

perm = mu0*((l/3)*(hs/wst) + hd/wst); 

Las = 2*p*Lst*perm*(4*Nc'^2*(m-Nsp)+2*Nsp*Nc^2); 

Lam = 2*p*Lst*Nsp*Nc'^2*perm; 
if q == 3 

Lslot = Las + 2*Lam*cos(2*piyq); % 3 phase equation 

else 

Lslot = Las - 2*Lam*cos(2*pi/q); % multiple phases 

end 

% End-turn inductance (Hanselman) 

As = ws*hs; % Slot area 

Le = ((Nc*muO*(taus)*Na'^2)/2)*log(wt*sqrt(pi)/sqrt(2*As)); 

% Total inductance and reactance 
Ls = Lag+Lslot+Le; 

Xs = omega*Ls; 

% Lengths, Volumes, and Weights 

% Armature conductor length 
Lac = 2*Na*(Lst+2*le2); 

% Armature conductor area (assumes form wound) 

Aac = As*lams/(2*Nc); 

% Mass of armature conductor 
Mac = q*Lac*Aac*rhoc; 

% Overall machine length 
Lmach = Lst+2*lel; 

% Core inside radius 
Rci = R+hm+g+hd+hs; 

% Core outside radius 
Rco = Rci+dc; 

% Overall diameter 
Dmach = 2*Rco; 

% Core mass 

Mcb = rhos*pi*(Rco^2-Rci'^2)*Lst; % Back iron 

Met = rhos*Lst*(Ns*wt*hs+2*pi*R*hd-Ns*hd*wd); % Teeth 

Me = Mcb + Met; 

% Magnet mass 
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Mm = 0.5*(p*thmrad)*((R+hm)'^2“R^2)*Lst*rhom; 

% Shaft mass 

Ms = pi*R^2*Lst*rhos; 

% 15% service fraction 

Mser = 0.15*(Mc+Ms+Mm+Mac); 

% Total mass 

Mtot = Mser+Mc+Ms+Mm+Mac; 

% Stator resistance 
Ra = Lac/(sigst*Aac); 

% Core Loss Calculations 

% Tooth Flux Density 
Bt = Bg/tfrac; 

% Back iron flux density (Hanselman) 

Bb = Bg*R/(p*dc); 

% Core back iron loss 

Pcb = Mcb*PO*abs(Bb/BO)^epsb*abs(f/FO)''epsf; 

% Teeth Loss 

Pet = Mct*PO*abs(Bt/BO)''epsb*abs(f/FO)^epsf; 

% Total core loss 
Pc = Pcb + Pet; 

% Start loop to determine terminal voltage and current 
notdone = 1; 
i = 0; 

la = Pwr/(q*Ea); 
while notdone ==1 
i = i + 1; 
xa = Xs*Ia/Ea; 

% Conductor losses 
Pa = q*Ia^2*Ra; 


% Gap friction losses 
% Reynold's number in air gap 
omegam = omega/p; 

Rey = omegam*R*g/nuair; 

% Friction coefficient 
Cf = .0725/Rey^.2; 

% Windage losses 

Pwind = Cf*pi*rhoair*omegam'^3*RM*Lst; 

% Get terminal voltage 

Va = sqrt(Ea^2-((Xs+Ra)*Ia*cos(psi))^2)-(Xs+Ra)*Ia*sin(psi); 

Ptemp = q*Va*Ia*cos(psi)-Pwind; 
error = Pwr/Ptemp; 
err(i) = error; 
if abs(error-l) < tol 
notdone = 0; 
else 

la = Ia*error; 
end 
end 
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% Remaining performance parameters 

% Current density 
Ja = la/Aac; 

% Power and efficiency 
Pin = Pwr+Pc+Pa+Pwind; 
eff = Pwr/Pin; 
pf = cos(psi); 

fprintfCpmlcalc complete: Ready.\n'); 


% Jonathan Rucker, MIT Thesis 
% May 2005 
% Program: pm 1 output 
% Program outputs values from pmlcalc. 

% Program developed from J.L. Kirtley script with permission 

% MUST RUN pmlinput and pmlcalc PRIOR TO RUNNING pmloutput 


% Variables for output display 
Pout = Pwr/le3; 

Jao = Ja/le4; 

Pco = Pc/le3; 

Pwindo = Pwind/le3; 

Pao = Pa/le3; 
wso = ws*1000; 
hso = hs*1000; 
wto = wt*1000; 
dco = dc*1000; 

Lso = Ls^lOOO; 
hmo = hm*1000; 
go = g*1000; 


% Output Section: 

fprintf('\nPM Machine Design, Version 1: Surface Magnet, Slotted Stator\n’); 


fprintf(’Machine Size:\n'); 
fprintf(’Machine Diameter = %8.3f m 
fprintf(’Rotor radius = %8.3f m 

fprintf(‘Slot Avg Width = %8.3f mm 
fprintf('Back Iron Thick = %8.3f mm 

fprintf(Machine Ratings:\n'); 
fprintf(Power Rating = %8.1f kW 

fprintf( Va (RMS) = %8.0f V 

fprintfCEa (RMS) = %8.0f V 

fprintf(’Synch Reactance = %8.3f ohm 
fprintf('Stator Cur Den = %8.1f A/cm2 

f^rintfCEfficiency = %8.3f 

fprintf(’Phases = %8.0f 

fprintfCStator Parameters:\n'); 
fi)rintf(’Number of Slots = %8.0f 


Machine Length = %8.3f m\n',Dmach,Lmach); 
Active length = %8.3f m\n',R,Lst); 

Slot Height = %8.3f mm\n',wso,hso); 

Tooth Width = %8.3f mm\n',dco,wto); 


Speed = %8.0f RPM\n', Pout,rpm); 

Current = %8.1f A\n', Va,Ia); 

Arm Resistance = %8.5f ohm\n*,Ea,Ra); 
Synch Induct = %8.3f mH\n',Xs,Lso); 
Tip Speed = %8.0f m/sW, Jao,vtip); 
Power Factor = %8.3f\n', eff,pf); 
Frequency = %8. If Hz\n’,q,f); 


Num Arm Turns = %8.0f \n',Ns,Na); 
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fprintfCBreadth Factor = %8.3f 
fprintf(Tooth Flux Den = %8,2f T 
fprintf(’Slots/pole/phase = %8,2f\n‘,m); 

fprintf(’Rotor Parameters:\n'); 
fi)rintf('Magnet Height = %8.2f mm 

fprintf('Air gap = %8.2f mm 

fprintf('Magnet Remanence = %8.2fT 
fprintf('Magnet Factor = %8.3f 

fprintf('Machine Losses:\n'); 
fi)rintf('Core Loss = %8.1f kW 

fprintf(Windage Loss = %8.If kW 

fprintf('Machine Weights :\n’); 
fprintf('Core = %8.2fkg 

fprintf(’Magnet = %8.2f kg 

fprintf(’Services = %8.2f kg 


Pitch Factor = %8.3f \n', kb,kp); 
Back Iron - %8.2f T\n’, Bt,Bb); 


Magnet Angle = %8.1f degm\n',hmo,thm); 
Pole Pairs = %8.0f \n',go,p); 

Aig Gap Bg = %8.2f T\n',Br,Bg); 

Skew Factor = %8.3f \n’,kg,ks); 


Armature Loss = %8.1f kW\n', Pco,Pao); 
Rotor Loss = TBD kWVn', Pwindo); 


Shaft = %8.2f kg\n’,Mc,Ms); 

Armature = %8.2f kg\n’,Mm,Mac); 

Total = %8.2f kg\n',Mser,Mtot); 
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Appendix E. MATLAB Code: Sizing Method 2 


% Jonathan Rucker, MIT Thesis 
% May 2005 
% Program: pm2input 
% Program used as input file for pm2calc 
% All necessary input parameters entered here. 

clear; 


% Definition & Entry of variables 
% General variables 

Pwr = 16e6; % Required power (W) 

rpm = 13000; % Speed (RPM) 

psi = 0; % Power factor angle 

Bsat = 1.65; % Stator saturation flux density 

% Rotor variables 

vtip = 200; % Tip speed limit (m/s) 

p = 3; % Number of pole pairs 

Br = 1.2; % Magnet remnant flux density (T) 

thsk =10; % Magnet skew angle (elec deg) 

PC = 5.74; % Permeance coefficient for magnets 


% Stator variables 
Ja = 2200; 
q = 3; 
m = 2; 

Nsp = 1; 
g = .004; 
hs = .025; 
hd = .0005; 
wd = le-6; 
ws = .016; 

Nc= 1; 
lams = 0.5; 
sigst = 6.0e+7; 

% Densities 
rhos = 7700; 
rhom = 7400; 
rhoc = 8900; 


% Initial current density (A/cm2) 
% Number of phases 
% Slots/pole/phase 
% Number of slots short pitched 
% Air gap (m) 

% Slot depth (m) 

% Slot depression depth (m) 

% Slot depression width (m) 

% Avg slot width (m) 

% Turns per coil 
% Slot fill fraction 
% Stator winding conductivity 

% Steel density (kg/m3) 

% Magnet density (kg/m3) 

% Conductor density (kg/m3) 


% Jonathan Rucker, MIT Thesis 
% May 2005 
% Program: pm2calc 

% Program performs sizing and parameter calculations 
% for permanent magnet machines with surface magnets and 
% slotted stators. 

% MUST RUN pm2input PRIOR TO RUNNING pm2calc 


%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 

% Definition of variables 
% Name Variable 
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% General variables 


% Pwr 

Required power (W) 

% rpm 

Speed (RPM) 

% psi 

Power factor angle 

%f 

Electrical frequency (Hz) 

% omega 

Electrical frequency (rad/sec) 

% vtip 

Tip speed (m/s) 

% lambda 

Flux linkage 

%Ea 

% Rotor variables 

RMS Internal voltage (V) 

%R 

Rotor radius (m) 

% hm 

Magnet thickness (m) 

%Lst 

Rotor stack length (m) 

%p 

Number of pole pairs 

%Br 

Magnet remnant flux density (T) 

% thm 

Magnet physical angle (deg) 

% thsk 

% Stator variables 

Magnet skew angle (actual deg) 

%q 

Number of phases 

% m 

Slots per pole per phase 

%Ns 

Number of slots 

%Nsp 

Number of slots short pitched 

%g 

Air gap (m) 

% ge 

Effective air gap (m) 

% tfrac 

Peripheral tooth fraction 

% hs 

Slot depth (m) 

%hd 

Slot depression depth (m) 

% wd 

Slot depression width (m) 

% syrat 

Stator back iron ratio (yoke thick/rotor radius) 

%Nc 

Turns per coil 

% lams 

Slot fill fraction 

% sigst 

Stator conductivity 

%Kc 

% Loss Models 

Carter coefficient 

%P0 

Base power for core losses 

%F0 

Base frequency for core loss 

%B0 

Base flux density 

% epsb 

Flux density exponent 

% epsf 

Frequency exponent 

% rhos 

Steel density 

% rhom 

Magnet density 

% rhoc 

Conductor density 

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 

% Constants to be used 

mu0 = 4*pi*le-7; 

% Free space permeability 

tol = le-2; 

% Tolerance factor 

cpair = 1005.7; 

% Specific heat capacity of air (J/kg*C) 

rhoair = 1.205 ; 

% Density of air at 20 C (kg/m3) 

nuair = 1.5e-5; 

% Kinematic viscosity of air at 20 C (m2/s) 

P0 = 36.79; 

% Base Power Losss, W/lb 

FO = 1000; 

% Base freuency, 60 Hz 

BO =1.0; 
epsb = 2.12; 
epsf = 1.68; 

% Base flux density, 1.0 T 
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% Calculate electrical frequency & rotor radius 
f = p*rpni/60; 
omega = 2*pi*f; 

R = p*vtip/omega; 

% Winding & skew factors 

Ns = floor(2*q*p*m); % Number of slots 

gama = 2*pi*p/Ns; 

Nsfp = floor(Ns/(2*p)); 

Nsct = Nsfp - Nsp; 

alfa = pi*Nsct/Nsfp; 

kp = sin(pl/2)*sin(alfa/2); 

kb = sin(m*gama/2)/(m*sin(gama/2)); 

kw = kp*kb; 

ths = ((p*thsk)+le-6)*(pi/180); , % skew angle (elec rad) 

ks = sin(ths/2)/(ths/2); 

% Calculate magnet dimensions, tooth width, & air gap flux density 
thme =1; % Initial Magnet angle (deg e) 

notdone = 1; 

ge = g; % Initial effective air gap 

while notdone == 1 

alpham = thme/180; % Pitch coverage coefficient 

Cphi = (2*alpham)/(l+alpham); % Flux concentration factor 
hm = ge*Cphi*PC; % Magnet height 

Ds = 2*(R+hm+g); % Inner stator/air gap diameter 

K1 = 0.95; % Leakage factor 

Kr = 1.05; % Reluctance factor 

murec =1.05; % Recoil permeability 

Bg = ((Kl*Cphi)/(l+(Kr*murec/PC)))*Br; 
wt = ((pi*Ds)/Ns)*(Bg/Bsat); % Tooth width 
taus = ws + wt; % Width of slot and tooth 

Kc= l/(l-(l/((taus/ws)*((5*g/ws)+l)))); % Carter's coefficient 
ge = Kc*g; 
eratio = ws/wt; 
if abs(eratio - 1) < tol 
notdone = 0; 
else 

thme = thme + 1; 
end 
end 

% Set final values 
thm = thme/p; 
thmrad = thm*(pi/180); 
hm = ge*Cphi*PC; 

Ds = 2*(R+hm+g); 

% Generate geometry of machine 

% Peripheral tooth fraction 
tfrac = wt/(wt+ws); 

% Slot top width (at air gap) 
wst= 2*pi*(R+g+hm+hd)*tfrac/Ns; 

% Slot bottom width 

wsb = wst*(R+g+hm+hd+hs)/(R4-g+hm+hd); 


% Magnet physical angle 

% Magnet height 
% Inner stator/air gap diameter 
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% Stator core back iron depth 

dc = (pi*Ds*thmrad/(4*p))*(Bg/Bsat); 

% Core inside radius 
Rci = R+hm+g+hd+hs; 

% Core outside radius 
Rco = Rci+dc; 

% Slot area 
As = ws*hs; 

% Estimate end turn length 
% End turn travel (one end) 
laz = pi*(R+g+hm+hd+0.5*hs)*Nsct/Ns; 

% End length (half coil) 
le2 = pi*laz; 

% End length (axial direction) 
lei = 2*le2/(2*pi); 

% Calculate magnetic gap factor 
Rs = R+hm+g; 

Ri = R; 

R1=R; 

R2 = R+hm; 

kg = ((Ri^p4))/(Rs'^(2*p)-Ri'^(2*p)))*((p/(p+l))*(R2^(p+l)-Rl^(p+l))... 
+(p*Rs^(2*p)/(p-l))*(Rl'^(l-p)-R2^(l-p))); 

% Core loss calculations (per length) 

% Core mass per length 
McbperL = rhos*pi*(Rco^2-Rci'^2); 

MctperL = rhos*(Ns*wt*hs+2*pi*R*hd-Ns*hd*wd); 

McperL = McbperL + MctperL; 

% Tooth Flux Density 
Bt = Bg/tfrac; 

% Back iron flux density (Hanselman) 

Bb = Bg*R/(p*dc); 

% Core back iron loss per length 

PcbperL = McbperL*PO*abs(Bb/BO)^epsb*abs(£/FO)^epsf; 

% Teeth Loss per length 

PctperL = MctperL *PO*abs(Bt/BO)^epsb*abs(f/FO)^epsf; 

% Total core loss per length 
PcperL = PcbperL + PctperL; 

% Current and surface current density 
% Armature turns (each slot has 2 half coils) 

Na = 2*p*m*Nc; 

% Arm cond area (assumes form wound) 

Aac = (As*lams)/(2*Nc); 

% Power & Current waveform factors (Lipo) 
ke = 0.52; 
ki = sqrt(2); 

% Initial terminal current 
la = Ns*lams*As*Ja*le4/(2*q*Na); 
notfin = 1; 

Lst = 0.1; % Initial stack length 

i=l; 


% Back iron 
% Teeth 
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% Start loop to determine Lst, Ea, Va, and la 
notdone = 1; 
k = 0; 

while notdone == 1 
k = k+l; 

% Surface current density 
A = 2*q*Na*Ia/(pi*Ds); 


% Calculate stack length of machine 
% Loop to get stack length 
while notfin == 1 
% Gap power 

Pgap = 4*pi*ke*ki*kw*ks*kg*sin(thmrad)*(f/p)*A*Bg*(Ds^2)*Lst; 
% Length of conductor 
Lac = 2*Na*(Lst+2*le2); 

% Stator resistance 
Ra = Lac/(sigst*Aac); 

% Copper Loss 
Pa = q*Ia^2*Ra; 

% Core losses 
Pc = PcperL*Lst; 

% Iterate to get length 
Ptempl = Pgap-Pa-Pc; 
error = P wr/Ptemp 1; 
err(i) = error; 
if abs(error-l) < tol 
notfin = 0; 
else 

Lst = Lst*error; 
i = i+ 1; 
end 
end 

% Calculate magnetic flux and internal voltage 
thmrad = thm*(pi/180); 

B1 = (4/pi)*Bg*kg*sin(p*thmrad/2); 
lambda = 2*Rs*Lst*Na*kw*ks*Bl/p; 

Ea = omega*lambda/sqrt(2); % RMS back voltage 

% Calculation of inductances/reactances 

% Air-gap inductance 

Lag = (q/2)*(4/pi)*(muO*Na''2*kw^2*Lst*Rs)/(p^2*(g+hm)); 

% Slot leakage inductance 

perm = mu0*((l/3)*(hs/wst) + hd/wst); 

Las = 2*p*Lst*perm*(4*Nc^2*(m-Nsp)+2*Nsp*Nc^2); 

Lam = 2*p*Lst*Nsp*Nc^2*perm; 
if q == 3 

Lslot = Las + 2*Lam*cos(2*pi/q); % 3 phase equation 

else 

Lslot = Las - 2*Lam*cos(2*pi/q); % multiple phases 

end 

% End-turn inductance (Hanselman) 

taus = ws + wt; % Width of slot and tooth 

Le = ((Nc*muO*(taus)*Na^2)/2)*log(wt*sqrt(pi)/sqrt(2*As)); 
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% Total inductance and reactance 
Ls = Lag+Lslot+Le; 

Xs = omega*Ls; 

% Lengths, Volumes, and Weights 

% Armature conductor length 
Lac = 2*Na*(Lst+2*le2); 

% Mass of armature conductor 
Mac = q*Lac*Aac*rhoc; 

% Overall machine length 
Lmach = Lst+2*lel; 

% Overall diameter 
Dmach = 2*Rco; 

% Core mass 
Me = McperL*Lst; 

% Magnet mass 

Mm = 0.5*(p*thmrad)*((R+hm)'^2-R^2)*Lst*rhom; 

% Shaft mass 

Ms = pi*R^2*Lst*rhos; 

% 15% service fraction 

Mser = 0.15*(Mc+Ms+Mm+Mac); 

% Total mass 

Mtot = Mser+Mc+Ms+Mm+Mac; 

% Gap friction losses 
% Reynold's number in air gap 
omegam = omega/p; 

Rey = omegam*R*g/nuair; 

% Friction coefficient 
Cf = m25/Rcy\2\ 

% Windage losses 

Pwind = Cf*pi*rhoair*omegam'^3*RM*Lst; 

% Get terminal voltage 
xa = Xs*Ia/Ea; 

Va = sqrt(Ea^2-((Xs+Ra)*Ia*cos(psi))^2)-(Xs+Ra)*Ia*sin(psi); 

Ptemp = q*Va*Ia*cos(psi)“Pwind; 

Perror = Pwr/Ptemp; 

Perr(k) = Perror; 
if abs(Perror--l) < tol 
notdone = 0; 
else 

la = Ia*Perror; 
end 
end 

% Remaining performance parameters 

% Current density 
Ja = la/Aac; 

% Power and efficiency 
Pin = Pwr+Pc+Pa+Pwind; 
eff = Pwr/Pin; 
pf = cos(psi); 
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fprintf('pm2calc complete: Ready.\n'); 


% Jonathan Rucker, MIT Thesis 
% May 2005 
% Program: pm2output 
% Program outputs values from pm2calc. 

% Program developed from J.L. Kirtley script with permission 

% MUST RUN pm2input and pm2calc PRIOR TO RUNNING pm2output 

% Variables for output display 
Pout = Pwr/le3; 

Jao = Ja/le4; 

Pco = Pc/le3; 

Pwindo = Pwind/le3; 

Pao = Pa/le3; 
wso = ws*1000; 
hso = hs*1000; 
wto = wt*1000; 
dco = dc*1000; 

Lso = Ls*1000; 
hmo = hm*1000; 
go = g*1000; 

% Output Section: 

fprintfOVnPM Machine Design, Version 2: Surface Magnet, Slotted Stator\n'); 
fprintf(Machine Size:\n'); 

fprintf(Machine Diameter = %8.3f m Machine Length = %8.3f m\n’,Dmach,Lmach); 

fprintf('Rotor radius = %8.3f m Active length = %8,3f mVn',R,Lst); 

fprintf(’Slot Avg Width = %8.3f mm Slot Height = %8.3f mm\n’,wso,hso); 

fprintf('Back Iron Thick = %8.3f mm Tooth Width = %8.3f mm\n',dco,wto); 

fprintf(Machine Ratings:\n'); 

f^rintf(Power Rating = %8.If kW Speed = %8.0f RPM\n', Pout,rpm); 

fprintf( Va (RMS) = %8.0f V Current = %8. If A\n*, Va,Ia); 

fprintfCEa (RMS) = %8.0f V Arm Resistance = %8.5f ohm\n',Ea,Ra); 

fprintf('Synch Reactance = %8.3f ohm Synch Induct = %8.3f mH\n',Xs,Lso); 

fprintf(’Stator Cur Den = %8.1fA/cm2 Tip Speed = %8.0f m/sW, Jao,vtip); 

fprintf(’Efficiency = %8.3f Power Factor = %8.3f\n', eff,pf); 

fprintf(’Phases = %8.0f Frequency = %8.1f Hz\n’,q,f); 

fprintf('Stator Parameters:\n'); 

fprintf('Number of Slots = %8.0f Num Arm Turns = %8.0f \n',Ns,Na); 

fprintf(’Breadth Factor = %8.3f Pitch Factor = %8.3f \n', kb,kp); 

fprintfCTooth Flux Den = %8.2fT Back Iron = %8.2f T\n', Bt,Bb); 

fprintf(’Slots/pole/phase = %8.2f^',m); 

fprintf('Rotor Parameters:\n‘); 

fprintf(’Magnet Height = %8.2f mm Magnet Angle = %8.1f degm\n',hmo,thm); 

fi)rintf(’Air gap = %8.2fmm Pole Pairs = %8.0f \n’,go,p); 

fprintf(Magnet Remanence = %8.2fT AigGapBg= %8.2f T\n’,Br,Bg); 
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fprintf('Magnet Factor = %8.3f 


Skew Factor = %8.3f \n*,kg,ks); 


fprintf('Machine Losses:\n'); 

fprintf('Core Loss = %8.1fkW Armature Loss = %8.1f kW\n', Pco,Pao); 

fprintf('Windage Loss = %8.1f kW Rotor Loss = TBD kW\n', Pwindo); 

fprintf(’Machine Weights:\n'); 


fprintf(’Core = 

%8.2fkg 

Shaft = 

%8.2fkg\n',Mc,Ms); 

fprintf('Magnet = 

%8.2fkg 

Armature = 

%8.2f kg\n',Mm,Mac); 

fprintf(’Services = 

%8.2fkg 

Total = 

%8.2f kg\n',Mser,Mtot); 
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Appendix F. MATLAB Code: Bode Plot 


% Jonathan Rucker, MIT Thesis 
% May 2005 
% Program: Buckfilter 

% Program calculates transfer function and outputs 
% Bode plot for buck converter input filter 

clear; 

% Input parameters 
R = 4.79e-3; 

CfF2.84e-3; 

Cb = 28.4e-3; 

Lf=1.415e-6; 

% Set up transfer function 
num=[R*Cb 1]; 

den = [Lf*R*Cf*Cb Lf^^CCf+Cb) R*Cb 1]; 

H = tf(num,den); 
bode(H) 
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Appendix G. 

MATLAB Code: PM Generator Waveforms 

% Jonathan Rucker, MIT Thesis 

% May 2005 


% Program: pmwave 


% Program calculates and outputs different waveforms, 

% calculates THD, and computes the harmonic content 

% for permanent magnet machines with surface magnets and 

% slotted stators. 


% MUST RUN pmlinput and pmlcalc PRIOR TO RUNNING pmwave 

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 

% Definition of variables 


% Name 

Variable 

% General variables 


% Pwr 

Required power (W) 

% rpm 

Speed (RPM) 

% psi 

Power factor angle 

%f 

Electrical frequency (Hz) 

% omega 

Electrical frequency (rad/sec) 

% vtip 

Tip speed (m/s) 

% lambda 

Flux linkage 

%Ea 

RMS Internal voltage (V) 

% Rotor variables 


%R 

Rotor radius (m) 

% hm 

Magnet thickness (m) 

% Lst 

Rotor stack length (m) 

%p 

Number of pole pairs 

%Br 

Magnet remnant flux density (T) 

% thm 

Magnet physical angle (deg) 

% thsk 

Magnet skew angle (actual deg) 

% Stator variables 


%q 

Number of phases 

% m 

Slots per pole per phase 

%Ns 

Number of slots 

% Nsp 

Number of slots short pitched 

%g 

Air gap (m) 

% ge 

Effective air gap (m) 

% tfrac 

Peripheral tooth fraction 

%hs 

Slot depth (m) 

%hd 

Slot depression depth (m) 

% wd 

Slot depression width (m) 

% syrat 

Stator back iron ratio (yoke thick/rotor radius) 

%Nc 

Turns per coil 

% lams 

Slot fill fraction 

% sigst 

Stator conductivity 

%Kc 

Carter coefficient 

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 

% Constants to be used 


inu0 = 4*pi*le-7; 

% Free space permeability 

tol = le-2; 

% Tolerance factor 
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% Harmonics to be evaluated 
n = 1:2:35; 

np = p .* n; % Use in kgn equation 

w = n .* omega; % Harmonic angular frequencies 

freq = w ./ (2*pi); % Harmonic frequencies 

% Harmonic winding and skew factors 

gama = 2*pi*p/Ns; 

alfa = pi*Nsct/Nsfp; 

kpn = sin(n .* pi/2) .* sin(n .* alfa/2); 

kbn = sin(n .* m*gama/2) J (m*sin(n .* gama/2)); 

kwn = kpn .* kbn; 

ths = ((p*thsk)+le-6)*(pi/180); % skew angle (elec rad) 

ksn = sin(n .* ths/2)./ (n .* ths/2); 

% Calculate magnetic gap factor 
Rs = R+hm+g; 

Ri = R; 

R1=R; 

R2 = R+hm; 

kgn = ((Ri.^(np4))./(Rs.^(2.*np)~Ri.^(2.*np))) *((np./(np+l))... 
.*(R2.^(np+l)^Rl.^(np+l))+(np*Rs.^(2.*np)y(np-l))... 
.*(RL^(Unp)-R2.'^(Unp))); 

% Calculate magnetic flux and internal voltage 
thmrad = thm*(pi/180); 
thmerad = p*thmrad; 

Bn = Bg.*((4/pi)7n).*kgn.*sin(n.*thmerad/2).*sin(n.*pi/2); 
lambdan = ((2*Rs*Lst*Na).*kwn.*ksn.*Bn)./p; 

Ban = (omega. *lambdan); % Peak back voltage 

% Normalized values for plotting 
Eanorm = abs(Ean)./ Ean(l); 

% Voltage THD 
Eah = 0; 

for r = 2:length(n) 

Eah = Eah + Ean(r)^2; 
end 

THD= 100*sqrt(Eah/(Ean(l)^2)); 

% Generate waveforms 

% Rotor physical angle goes from 0 to 2*pi - electrical to 2*p*pi 
ang = 0:pi/100:2*pi; 
angp = p*ang; 

Bout = zeros(size(angp)); 

Eaout = zeros(size(angp)); 
for i = 1 :length(n) 

Bout = Bout + Bn(i).*sin(n(i).*angp); 

Eaout = Eaout + Ean(i).*sin(n(i).*angp); 
end 

% Plot waveforms 

figure(l) 

plot(ang,Bout); 

title(’PM Generator: Flux Density'); 
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ylabeK’B (Tesla)'); 
xlabeK'Rotor Angle (rad)'); 

figure(2) 

plot(ang,Eaout); 

title(’PM Generator: Back EMF'); 
ylabeK’Peak Voltage (V)'); 
xlabeK'Rotor Angle (rad)'); 

figure(3) 
hold on 

title(['PM Generator: EMF Harmonics, THD = ',num2str(THD),' %']); 
ylabeK’Normalized Back EMF'); 
xlabeK'Harmonic Number'); 

text(20,0.7,'Dark: Above 10% of Fundamental','FontSize', 10); 
text(20,0.65,'Light: Below 10% of Fundamental',EontSize', 10); 
for z = l:length(Eanorm) 
if Eanorm(z) < 0.10 
bar(n(z),Eanorm(z),'c'); 
else 

bar(n(z),Eanorm(z),'b'); 

end 

end 

hold off 
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Appendix H. MATLAB Code: Retaining Sleeve Stress Calculations 


% Jonathan Rucker, MIT Thesis 
% May 2005 
% Program: pmcanstress 

% Program calculates and outputs retaining can stress 
% for permanent magnet machines with surface magnets and 
% slotted stators. 

% MUST RUN pmlinput, pmlcalc, and pmwave PRIOR TO RUNNING pmcanstress 


%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 
% Calculate retaining sleeve stresses 

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 
% Conversion 

Patopsi = 1.45038e-4; % psi per Pa 

% Material yield stresses (ksi) 

Stain_str = 90; 

Alum_str = 75; 

Titan_str = 110; 

CarFib_str= 100; 

InconeLstr = 132; 

% Safety factor 
SF=:1.2; 

% Force on magnets/sleeves is centrifugal force 
% Magnet tangential velocity 
vmag = ((R+hm)*omega)/p; 

% Centrifugal force 

Fm = (Mm*vmag^2)/(R+hm); 

% Outward pressure 

Phoop = Fm^(2*pi*(R+hm)*Lst); 

% Hoop Stress (in general, str = P*R/t) 
stop = 22; 
for i = listop 

t(i) = i*.0005; % sleeve thickness t 

slev(i) = t(i)*1000; 

Sthoop(i) = (Phoop*(R+hm)/t(i))*Patopsi/1000; 

SFHoop(i) = Sthoop(i)*SF; 
end 

% Output results 

fjprintfORetaining Sleeve Stress:\n'); 

^rintf(’Stress Limits:\n'); 

fprintf(’Stainless Steel = %6.1f ksi Aluminum Alloy = %6.1f ksi\n',Stain_str,Alum_str); 
fprintf(Titanium Alloy = %6.1fksi Carbon Fiber = %6.1f ksi\n',Titan_str,CarFib_str); 
fprintf('Inconel = %6.1 f ksi\n',InconeLstr); 

fprintf(Actual Sleeve Stress:\n‘); 

fprintf('Sleeve Thickness Actual Stress SF Stress\n'); 
for i = l:stop 

fprintf(’ %5.2fmm %6.1fksi %6.1fksi\n',... 
slev(i),Sthoop(i),SFHoop(i)); 
end 
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Appendix 1. MATLAB Code: Rotor Losses from Winding Time and 
Space Harmonics 

% Jonathan Rucker, MIT Thesis 
% May 2005 . 

% Program: pmharmloss 

% Program performs rotor loss calculations caused by 
% winding time and space harmonics for permanent magnet 
% machines with surface magnets and slotted stators. 

% MUST RUN pmlinput, pmlcalc, pmloutput, and get harmonic 
% current data from PSIM prior to running pmharmloss 

% Constants to be used 

muO = 4*pi*le-7; % Free space permeability 

tol = le-2; % Tolerance factor 

% Retaining sleeve/magnet material resistivity (ohm-m) 

Stain_res = 0.72e-6; 

Titan_res = 0.78e“6; 

CarFib_res = 9.25e-6; 

InconeLres = 0.98e>6; 

Magnet_res = 1.43e-6; 

% Retaining sleeve thickness set at 0.5mm less than air gap 
h_sl = g - 0.0005; % Sleeve thickness 

g_act = g - h_sl; % Actual air gap 

% Retaining sleeve conductivities (S/m) 

cond„s = 1/Stain_res; 

cond_t = 1/Titan_res; 

cond_c = l/CarFib_res; 

cond_i = l/Inconel_res; 

% Magnet & actual sleeve cond (S/m) 
cond_m = 1/Magnet_res; 
cond^sl = cond_s; 

% Input time harmonic peak currents from PSIM 
11 =2895; 

13 = 0; 

15 = 209; 

17 = 89.2; 

19 = 0; 

111=39.2; 

113 = 27.6; 

115=0; 

117 = 17.0; 

119=12.8; 

121 = 0 ; 

123 = 7.3; 

125 = 6.3; 

127 = 0; 
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129 = 4.8; 

131=3.9; 

% Put currents in array 

Iharm = [II 13 15 17 19II1 113 115 117 119121 123 125... 

127 129131]; 

% Calculate current THD 
lah = 0; 

for r = 2:length(Iharm) 
lah = lah + Iharm(r)^2; 
end 

THDi = 100*sqrt(Iah/(Iharm(l)^2)); 

% Calculate current densities 
Iz = (l/sqrt(2)).*Iharm; 

Kz = ((q/2)*(Na/(2*pi*Rs))).*Iz; 

Iz^l = (l/sqrt(2)).*Il; 

Kz^l = ((q/2)*(Na/(2*pi*Rs))).*Iz_l; 

% Harmonics to be evaluated 
n= 1:2:31; ' 
w = n .* omega; 
freq = w ./ (2*pi); 
lam = (2*(2*pi/(2*p)))./n; 
k = (2*pi)./lam; 

% Eta values 
eta_m = sqrt((j*muO*cond_m).*w + (k.'^2)); 
eta_s = sqrt((j*muO*cond_sl).*w + (k.'^2)); 

% Surface coefficient at top of magnet layer 
alpha_m = j. *(k./eta_m). *coth(eta_m. *hm); 

% Surface coefficient at top of retaining sleeve 

topi = (j.*(k./eta_s).*sinh(eta„s.*h_sl)) + (alpha_m.*cosh(eta_s.*h_sl)); 
botl = (j.*(k./eta_s).*cosh(eta_s.*h_sl)) + (alpha_m.*sinh(eta_s.*h_sl)); 
alpha_s = j. *(k./eta_s). *(top 1 ./bot 1); 

% Surface coefficient at surface of stator 
top2 = (j.*sinh(k.*g_act)) + (alpha_s.*cosh(k.*g_act)); 
bot2 = (j.*cosh(k.*g_act)) + (alpha_s.*sinh(k.*g_act)); 
alpha_f = j.*(top2./bot2); 

% Surface impedance 
Zs = (muO.*w./k).*alpha_f; 

% Calculate losses due to time harmonics 
% Use only fundamental space harmonic factors 
Kz_t = kw.*Kz; 

Syt = 0; 

fori= l:length(n) 

Sy^t(i) = 0.5*(abs(Kz_t(i))''2)*real(Zs(l)); 

Syt = Syt + Sy_t(i); 
end 


% Harmonic angular frequencies 
% Harmonic frequencies 

% Wavenumbers 


% Fundamental RMS current 
% Fundamental current density 
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% Calculate losses due to space harmonics 
% Use only fundamental time harmonic current 
kpn = sin(n .* pi/2) .* sin(n alfa/2); 
kbn = sin(n .* m*gama/2) 7 (m.*sin(n .* gama/2)); 
kwn = kpn .* kbn; 

Kz_s = kwn .* Kz_l 7 n; 

Sys = 0; 

for i = l:length(n) 

Sy_s(i) = 0.5*(abs(Kz_s(i))^2)*real(Zs(i)); 

Sys = Sys + Sy_s(i); 
end 

fprintf(’\nRotor Losses Caused by Harmonics:\n'); 
fprintf(Time Harmonic Losses = %6.1f kW\n',Syt/1000); 
fJ)rintf('Space Harmonic Losses = %6.1f kW\n',Sys/1000); 
fprintf(Total Losses = %6. If kW\n',(Syt+Sys)/1000); 

fprintf(Current THD = %6.2f %%\n’,THDi); 
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Appendix J. MATLAB Code: Rotor Losses from Slot Effects 


% Jonathan Rucker, MIT Thesis 
% May 2005 
% Program: pmcanloss 

% Program calculates and outputs rotor losses caused by 
% stator slot effects for permanent magnet machines 
% with surface magnets and slotted stators. 

% MUST RUN pmlinput, pmlcalc, pmwave, and pmcanstress 
% PRIOR TO RUNNING pmcanloss 


%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 
% Calculate retaining sleeve losses 

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 


% Retaining sleeve/magnet material resistivity (ohm-m) 

Stain_res = 0.72e-6; 

Titan_res = 0.78e“6; 

CarFib_res = 9.25e-6; 

Inconel_res = 0.98e-6; 

Magnet_res= 1.43e-6; 

% Calculate Bd as function of wst and wt (max 10% of Bg) 

Bd = (wst/wt)*0.1*Bg; 

% Calculate flux variation parameters 
beta = (wst/(2*pi*Rs))*2*pi; 
lamB = 2*pi/Ns; 

B = (Bd/sqrt(2))*sqrt(beta/lamB); 

% Calculate geometry and can loss factor for different rings 

% k is number of rings 

fork=l:10 

A(k) = pi*2*(R+hm)*Lst/k; 

Ks(k) = 1 - ((tanh(p*Lst/(k*2*(R+hm))))/(p*Lst/(k*2*(R+hm)))); 
end 

% Input Stainless Steel sleeve thickness based on stress results 
for i = l:stop 

if SFHoop(i) <= Stain_str 
t_Stain = t(i); 
break 

elseif t(stop) > Stain_str 

fprintf(’Hoop Stress too high for Stainless Steel.Xn'); 
else 

dummy = t(i); 
end 
end 

% Input Titanium sleeve thickness based on stress results 
for i = 1 :stop 

if SFHoop(i) <= Titan_str 
t_Titan = t(i); 
break 

elseif t(stop) > Titan_str 
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fprintfCHoop Stress too high for Titanium.Xn’); 
else 

dummy = t(i); 
end 
end 

% Input Carbon Fiber sleeve thickness based on stress results 
fori= l:stop 

if SFHoop(i) <= CarFib_str 
LCarFib = t(i); 
break 

elseif t(stop) > CarFib_str 
fprintf('Hoop Stress too high for Carbon Fiber.\n'); 
else 

dummy = t(i); 
end 
end 

% Input Inconel sleeve thickness based on stress results 
for i = l:stop 

if SFHoop(i) <= InconeLstr 
tjnconel = t(i); 
break 

elseif t(stop) > InconeLstr 
fprintf('Hoop Stress too high for Inconel.Vn'); 
else 

dummy = t(i); 
end 
end 

% Calculate can losses 

w_Stain = (pi'^2/3600)*((B*rpm*(R+hm))'^2*t_Stain)/Stain„res; 
w_Titan = (pi^2/3600)*((B*rpm*(R+hm))'^2*LTitan)/Titan_res; 
w_CarFib = (pi^2/3600)*((B*rpm*(R+hm))^2*t„CarFib)/CarFib„res; 
wjnconel = (pi'^2/3600)*((B*rpm*(R+hm))^2*t_Inconel)/Inconel_res; 

fork=l:10 

P^Stain(k) = k*w_Stain*Ks(k)*A(k)/1000; 

P_Titan(k) = k*w^Titan*Ks(k)*A(k)/1000; 

P^CarFib(k) = k*w^CarFib*Ks(k)*A(k)/1000; 

P_Inconel(k) = k*w Jnconel*Ks(k)*A(k)/1000; 
end 

% Calculate magnet losses (only with carbon steel) 

% Calculate geometry and can loss factor 
Am = pi*2*R*Lst; 

Ksm = 1 - ((tanh(p*Lst/(2*R)))/(p*Lst/(2*R))); 

% Calculate magnet losses 

% Assumes only 10% of magnet thickness sees effects 
w_Magnet = (pi'^2/3600)*((B*rpm*R)^2*0.1*hm)/MagneLres; 
P_Magnet = w_Magnet*Ksm* Am/1000; 

% Output Results 
z = [15 10]; 

fprintf('Retaining Can LossesAn'); 
for i=l:3 
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k = z(i); 

fprintf(’% 1 .Of Rings:\n',k) 

fprintf(’Material Thickness Can Loss\n'); 

fprintfCStainless Steel %5.2f mm %6.1f kW\n’,t_Stain*1000,P_Stain(k)); 
fprintf(Titanium %5.2f mm %6.1f kW\n',t_Titan*1000,P_Titan(k)); 
fprintf('Carbon Fiber %5.2f mm %6,lf kW\n’,t_CarFib*1000,P^CarFib(k)); 
fprintf(’ Associated Magnet Loss %6.1f kW\n',P_Magnet); 
fprintf(’Inconel %5.2f mm %6.1f kW\n\n',tJnconel*1000,P_Inconel(k)); 

end 
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Appendix K. Results for PM Machine Variants 


Main Parameters 




Machine Parameter Data 


0ii(d«g) 0s(<*®g) 



_ _UiHa ■ ■rfdtltl ■.-jiliM TilM ■RaPM ■! »>■ WX!^ glPHia KtliWJ gngna J 

■ !«■<«■ mJB KWcicl gliSIM ■cjHtiM T^TiM gtaBM ■ WPTffB Btl;!;! ■ ■^»kgl Bif :l;yJ 

MEM BBM E*l^M gEfai i Kiltkicj ■>;JiM MigPM gtliMtl MBIga ■!>;!:?J Ktliiga ■iMEl ■«)[ 

Ki!«m k^m gygsgi ■iii?.; ■ ■ WciiiiM MnnM ■♦-n lm Btiiiigi «aiM Miiti ■ KiM mw^fum Bti:!;! ■ KKtiiiai mti^m 

' »jri!iM 1| MllimiimillHB'IBM|i||iHII'lHIIHB^P 

ggrna ntiiM ^wsm ^yna mm^ ■ .g:W>^ MnpgM ■iihkgg B O Eiigi 

■m>T ~ ".. .. .. .. _; 

fc KMM WBSM ■tUMJM tiatM 


Mc!;M MkM gc lt!:Iil ■!> T»IiM ■ VkWM ■»!;!:« »I»M MfrJtiM ■ K<kkM J ■*!;!;■ ■iliM 

_ Bjt^ KqniaB ; ! y3 aB;gMB [.y« igE^ B*n aEEaM .-i.B iB ElJ» K :gM ii»Eag «»ryH iBEBiB .i--M: iaaa 

_ WK:VyM Mif<!rf>JI i|ti!;ic1Rl t WWISM ■il!:^il »H*M MiIM ■ K:VJ:M Ktlif;!J KtlRgEl 

McM mwism M um ■< Ic!:M giUigra l MriU g ■!!!].-« gglM m KtMil ■ grna ■ WtVjm MirigM 

■hE M gCTnsaB tm ii m^imm u r m j mmm 

IjWtItM BIPHB ga»M Tiy ■»I;n>il mVKM ■;IiM ■*!♦!>■ BiliVa BlgCTl ■iK.-M ■ ■JigM 

_ _ _Bjg WiM ■il>I.-<;a KTi!^ ■< ■ [>!»■ WT^WM gel »«■ fHEM K i;g f Bil!!;!;! Btl:!;; J gliTgtl ■>!.-! J 

ill 11 XII "'HP' BJP I MWI'i MI_|I.ViH'm ll.ii.iiBW I'MBBiUBB iill ilUKBi|i|i|l l|j.'_ll H N iB.8 '.1 


■ fcM;! gypa ■agigil ■.iJiM TiIiM i V£ WA ■iltiiiiB ■ I WiM ■ KiStM ■iI:!.^M ■VWZiB 

■cj.-ltM »^EM KcyiM BEVa IHlca M K>!tgM BF^l M*!:!:!--! ■*X;H1 

MUtyj;'! BEga iyt!:jcici grgtiM TiiB ■if.-r.-iJMcitga ■;!#■ m*v^ irrgB ■»>»?■ KtinTia KiiiitM 

■ K»x-]:1 BEHiB ■ Kt!i!il»j TiliM IMt ■:!#■ Mgitga KiIrtfcB M BtitUl KiliCT MiliVi ■ Mtlii ■ 

■ K-'icZ ■ mniai gHifti*! B-ititM miiEM miyita »!« m:ttm KtiiMta Miiiitf MtikPM 

■fcai BBSM EBiBbEB liMiai BEl>l'lB±MliBA>-MBE?a K^?-aEE^BM«!l WfmMVmMmfkm 


Kgkl-l Bfiggl H.liiclcl ■>atMMI.I.» metga MPyiMg.IigaMilTRlB.I;!.-; ■nr?M 

KieWJ Byg3Miib!i!iM gSCTMlil.M Il^gEa ■>t»VlBgg5R?a 

■jliaM^l«EE«5EM»I>EB IHgSaM*^liaB;5Ea ■MiMBPISai KffMMqRaBa-pa Blt=K^Bmi 


Kg;f;S« Bfg3aB:!.-!dha ■EgfMMliI.M BETEl WPWSM 

KJfa K»ii&SllK»:>>ElEMa[M Bfri^MBlPiaiKMa MfUiM MPFH WfUM 


nE3imEgaiE5H»^»ECTiHtyjy^Bfr^iag^«MiqEa»iP?^gmiwp^ 


gr.i.M ■<««■ mhem kk.WiM mr^ mwm 

BJUi;tJBESl PtM«M BriEB Bi^^ Bt!»E^ BglHI BiJt?! P»)l;Vitl BITiIgl ■ tEgiqi 

B|g!cj<lBEgaiiK»^»-^*l Ba*B Bg^ glgltl ■»»:?;!J ■»I:I«5M WTf^ 

■it!tyjlgE)3niggt!i!»Iil gt.-f:!cV/ailt)bl.-t.-l ■ggglMfc JTM **_»_♦* 


Et:!tCTg!gyrag;!.-ftiil Mj<jU TiEB V>tV:M ■»<■ Blhna ■Jtlilil ■tfiMil gtlihgl ■tglM »»/■ 

BU El M'MM BckittM M r»EB g*I:I>l gHM » MiB ■iM IM mPCT ■i)l:yi;l BlEina gtliria 

_ gjt tcM B» l»yJ i]r<!:!c|c| BidliM ■ [iltB HiWM ■tfiWtl gTHM ■.♦»!■ K*M«1 ■g;!ig ■ ■>»;!;?■ ■■>!<!!■ BgHil 

■«!:!Jl H M ■ HiliMBEBaIgtti!*til■HatJTEM WiKVm KtlilA-l■gnCTI■:»!■ ■*!!«fini^■g>!iVlKliWJ■ifrftB 
- »>■ ■.-tiXiM mutm Fwm u^hbum mnm wmim wimm Kfcrrfj ■.i.iigi gma gin?y mtum 

aatfel »5 KiM TE» lEKUta ■il!l;I>l ggiEM giMM Kiigir^ ■il;!Al ■tl;l«yj ■tIiSE l MtW 

rnt MiJikM BXI hm ■ g>!;!clcl MlHitm ■ t>It» g?f:!:5M MUIiai MritM giljilra WnTIiB MtfMM ■ itWtM ■nttVi.-l ■iIiH;l m*f^M 

_ BEina |Pc!*l>EJ ■[if.B wsmi Bg^lcM Blpnagg^ 

KEBWEEW*].-kl BEEl BI:};agif;Ka BlRilB ■>!!»■ 

BjU-kkl BIEM BcidUg^EW g*lM WngB Bg:!;k^gilil:y^ B»iy^M 

■jK;!»a BH^ ■gi]:kfc| MgatipijEM KKiTa ■tI:M^MK«IiM giEHia ■>K!k!;l 

Ijlkya BEHgiJ^W*MBirt1tM 1 <IiM MltEf tm-l EHtiM m-giMBtltiba ■K;kl:lMtI:gH|MKi!«!tl Bigigy^ 

By:glB BEai MRiM IjEB iinF^4 BJIM Kiliba Btlilia KiKi&g ■ glEnB Btlzlika ■»>icM 

ESM wmm ga*B jgiEggBEMij ggiM mitm bmi BMjja KiEga ■»««■ 

_ _ Mr ji BnCT g »i!:!cM mwim m riu gdgrsM ■tg.-r.-fti ■' >'»*■ m;Iim Ktiti im ■ kiIimj ■ g»M>i giiroa BiinTii ■»!;! j 

■KM ■ii>M ■Ji^ BEiCT wggiira «ot« ■r.iTiM icjtbiai g.iii.ij m mm mam Bi t-ktiM ■K.T«ki j rntf^vm m*t^m 


By pa g»^a MMthM B.jiliM TjEB ■il;i;t:l KEM »!>».■■ ■ g>^ gg^Fl KiIsIiM KtlEnB ■««■ 

MiSa B3KFI1BMAK ggEB mbhm KHiJiKI BEMJ BiifM BX^Sa BI^ Bl^itJ MEEl BMjjl BlPg» 


■»kt«M ■K»’:gin K#?Wi^ ■in'^il IWiTIM KiIijiM ■HmXM MiTiM 

■gkVK Ktlga ggc!i!tltl »»iM ■mfM IHIf.-tM ■i»!lA-K»ggi» gj J Klit-M MWnm ■i»;g{;K»g»IiaMiliIia»iI;H<l 

Bgpa BlF^a yyUfki^^ ■*!:agI*lM:l.a Blhyi«K BlP!ia BI;?iaK<f:IiaBlPPTBK»g;t!lK fiMB 

KB#i Bn^a gaa BEijam-n>a i*M<a bem Bi^iaBEagEiaFgiisi mmi 






■K-i»iiMSEai»K«i»y^l!EEaiM«ffwaMai 


wm nmm mwrmmnmmmwr yMUimm gmi 


BI#1 BX^lKtlv^EI gaaia ■g[^M»EK:aK^aM!8i£IdlKE551 BJP^ Ktl-ka mirm 
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EEIB 




■PM 




■iltba »Kia KigifaMgtNa MiibWM Ktitga Mtiira 


M iMi Biipa ■■ :ka giag»Bg . ! .a m * i rm 

■r^a BygaHTtligltl gKi:yKBI;t;aBgtI«kK KtliHiK Bl!!*!!! Mif 

Bjgy^K KiHgailtWtM BjiEBMMtM ■ jtllMKim-a B^CT aHhaKit-MJ mWiiimMthZJm ggiSa KiIITia ■*»;■ 





m 

IS 

m 

0s (deg) 

1^ 




■Pm 






v^^im ■«i:Ka mtm 


■K}.t^Biilk«Kr^!?g Mni»M KUnti^M Bf?!;!;! »>j:kty| ■HMtca ■K;i.i»Kii:MiB ■mn>f ■iuum KiUHtg 

Bf».»:l gI «glK»»!Alcl Bj*B MPfm Kf ^ KiliMtl ^3CMBl:!.a »KnaBI;Wcl liTigi KiggiK BTPa Mil;? J 

MHJBEaKtmq g?JM ■t>X»a EFGEa WiKm MM KFgia Bf-Ka BEEl 


Mcr#^B]tPa»;:aarJ McauM ||pi»x«M ■hhhm m^Am WBSM Maiai P>iHigM mW?jyM ■iKnt.-a ■qHM 


DEsamHig 








=51KE 


Bma ■«K»yi ■aa 


ganaiKlK**/^ 
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M[=!ll.'l?>ia iHMJ Bji EM KliEB iid»v^ K^Ea m^mnm BlPHa BKica BK»MJ 

FiiJCnaa iM KEEb Ea-aai BKiiiil gE^a BEPEM ■aaaBE^ 


BcaiM mhtm racma ■*!;?.■ mnm MiM Bm« K^r^giffa Bm<mETa 

■KJ;aK»ltU:l ■gt!;lcM U'l^iMMlilttM B;gi!«ya ■ilil.-^iK BJHtiB ■aM B»;y;il KTira Bl;l»a ■npai 
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Appendix L. Results for Power Conversion Module Variants 




Phases Poles R«cW1«i L,(uH) C,{mF) Q, Z, »»LC L,{uH) C,(mF) Cb (mF) n,(mn) LC (uH) (mF 

can (HMb; I •e®*"' I 



KEM EE:iaaia!a3^ gEM ■KiatMi lEagg ■sMsai 

EEggggM^gia »^nggM laaiaa ggEHEaiiEaaia wKit-^ MEga ■saaa ■tgaw ^aasaJ iteigai 


EgSaEMIEEM E^M Et^H EB33Eg ■JBasaiCN*]^'?^ glEiaMaa!^ fcifliM tiaglaaa teritcy 

EE^EEEElEgMEgM EB50EKBlE51i;HiB«g™igMjifcjm3igEEa^itMa»Ki*MtE35a»a£aJW£E!ia^ 


EEE5MKtiMdlBMEEEM»aHPiliafaaaigfcimjtilEiaaii£E£B£g»RU3i»;&ia«EEB^aiaiEaaa!ateEEiEiEM 

EBigiE^3BEaMEMIEEHEEBE51iI*M^»aMiiataMDElBiEb]aMEft£aEEMiBM^IEDHiE^ 

EEEaE^CTEaaEaaEmgliaiaa!iai;fe:ia»aEm^»Kiiategi£t«EEMiKlflJIg^^gaf!gj»a!i!M 

EEEMEB^ME^EMcJfei;«;gii»iiH^glEaii^!3iat>a^EEafcRiEMEnM»agiii^tifcgaiHlMgEa»aiirf»a!aM 

gtKfeMBEEBlEgMEfeMEBgiEl^l^l^mdUiaEiaKiEtMEiaaiMEEaUKdsaiagikaMEEia 

■gBigMKi ic|ciidiMjtfeM Maiaiiiaai;^BKiaaaingi^5ifcga3ijaMtit;il»a!iiMEgiM^KtJifeaa!idit-i;LiMlliEE^ 

■i:kUtH mtti ta 


Phase. "•^'"'"9 Poles 

Can (RMS) 




L, (uH) 

C, (mF) 


A n-»^aa"rt* ii i i mm i ii i mi 

G» 

z, 

s*LC 

Lf(uH) 

C, (mF) 

Cb (mF) 


miisyjnm pcmititgi HKHas^yj 

mi^vm g i?kkl;M ■■.-« JiM gar#j J<g| iHMdafig ■g8le-Mgg:ti;B glHbgM Mtlricl 
^^a Mte1:>ikM Kl Ki^L^J|Ayiai£I EBgaggIHiEBMEEM E£MIM-a iBaEaM 
^ANie-a in Egnra KijHsgiici m .-vj i wuMafig »gnai ■.-f.4.-»a »cgi«M tiHjgira MticiiiM 

_ _«eCTMEEHOlEn3jggnR^miEE53!!g »glfea gEMgEMBgaiIKS35a^EMgE^ 


rniMfy^m fcaagai WESEM Kimgtjti i:«;i an-i Biiiil3l!gj MDEI^ waa^ ■feataail W-fEfcyiq!t«I^ taiisdi 

MEM EggQ MEM ^^lMti<=saHi mng?^g »ai;ia!iirari=n»/4fcKuai«gifeMBEijauai»i;HMl^^teiiiiiM■£&£!■ 
ME^ EpgnM mfxm EGM latn>*i ■Kiiaai gi»gja»a toiisM EBaa gaaiM utE^sm H^asa ■aaiJ ■asEHl 

MKiM EraaEaMEM EgM KfcjigtMEgg^gl MMjgiig mKum Msa^ BEB^B iiaaatia teig-MEBESi 

»g WcM M M fcfj • 1 i^ 1 r>ji;« ■ gi>] aigl K^iKW gifl WM M a 5 t >3 fcil Ktm fcf iWcM MEEaM ■gmi.'M 


Phases Poles 

Can (RMS) 


L,(uH) Cr(mF) 


Of Z, 8 “LC Lt (uH) C, (mF) Cb (mF) R, (mO) LC Lo„, (uH) Co„, (mF 



■nagaji^^M MiiLaM laiag-i tcjEgata meiu ■aasaw Ma!M n^^E MagMi^ ■aatlii 

MECTEn^igaMEgMg^MEtEadii;i£lM!!JIK»liEElElEH!!aMlME;BaiEEgiMaMt^^^jafeJteEMB 
EjEJIEEaa^^MEkMide^^AiiTO^EE^B^aEaMElIJiMsiaEIMMEMtEagaHEEniMiliai 
bagaataaiaaigi^^ig*^ c>tiiii=n*MMfclltaig^MtHt«y WifeadiMsftaa 

E>hagBBaiM5BME»kMEBgEiia^tia!aiiw^^'ia&aia2!aElEB»aajiMteitoiiMEIMigaaa*iafcHiaaE)ga 
MEMEBBaMElMggMEGB^KSiEaiME^ Effi ^TlPglWSip-aMfeEaigEEaiMftMiB-aaffijEjBlggai^ 
EslMEGEEaEMEgMEEEBlfcasSaiklfefeE^EEHaEgMlIdiMIfeMEEMMlMIfr-tf^^ERgiBHa-a 
M»IMEEgM«tMteMaEEEngiaiaa!aKKda»£i&aai!iaEM»£giawfcca»MMEMlEaa!ailE^£iaigEn»a 
mjmm ■n(»i:i>jKaMbEMi»aaKa isyaaw itgggE3»Bt?4:atjti la^.aTOi'^ Mma m^isM EfcU ■aaaw I^a*B ■aagJ gsi-aa 


»greME^n;CTaKnn» »ga Mici5Tn3ragjigi:fliViinE!gggigff!3g!«aMgnj»KcIsi£MF‘»‘»t^-*^*^™nE?^^^ g »aiflfli»a!iaa 
EEgM EBa-M gSaSW tss^ssU ■r«=fa>y]i Cfcaaiia teaa« ^stM IK-IM EEEa 

E<lME^gMIMiEMi£Kl!lJlfiISW>ailGHg>aEIfea!!aEiaiB§IMP;feIM»aMnEBga»aa«EEM 
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Appendix M. Results for Power Module Losses 
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Appendix N. Results for Power Module Weights 


Power Conversion Module 
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Appendix O. Rectifier/Input Filter Mass and Volume Calculations 
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I ■tH ■i»i!i7Ji gilgi^a I 

■ OtHPiO grliliVJ gl«!i;a I 


1 TOgHl p^naniM u .lil ilttliliin I 

I g»TRr^ EfiT^rnFl Mc'it’Jj ftl.-jcVAM UtMiM |_ 

I g»»f?Mo itggiin Eif=ngi EBggg mesem o*hmo ggpi^ga i 

I B:ir;«l giMryitl KTcTT^ MEIIPW Id^Jiklcl WIB3i ■ 

I Eg »:M lEgOB IIJi.Vl g1 111 1 B glTig^ MBiTiM ■iH l^M ■ilrliV ■tliliini | 

■ lajgEi giinncgi j M11.! L ■epi.-m vj^bm »■ ■ ^ ■iIKO ■iitiNv.a giitiinia | 


I iarM:»ti Ehtiia ■_ 

■ M*Vi !«!?■ ItIji SEl tM k!:» ■ ■<!;« ■ 


I Millpni KiKilM I 

■ Mil«!:t:a BiM-WtE ■ 


I ’flM IjM MrlErM PBcTifeM 


■ ■iMrVj mWm 


I Min:*vm » wmi ■igiM ■fc«:!.-«i 

I g«I?M ■e1P4dH 


Retaining „ . ^ p ■^®“ ^ ■^®^' Capadtor Auorago , J®*®* ^ ,J,®“ Inductor Dk>d. Inductor Capacitor Dlodo J®"*®* 

Can I W*^"*** ^ S**"*^ Cap^ltanco ^podt^ lnduc^o( Volumo Volumo Volumo Volumo Volumo ''®‘“f * 


Voltago I Cydo I Voltago 


uH) Enorgy(J) 


I Mifma ■a:M:ia WBUg^B MSSU WtEMEtE^I 


I Kil [;^0 ■»lt«!iVOI glg^ I 

■ gnTiiaia I 

■ OrurrM M*r*r*im g«n<!:iro ■ 

i lanmni ■.ir!riM liKua ■ 


I Eg!i!>H Dra:l-J M.-tr^l EfCSFl ■ 




I ia*T:?«Icl icWtVZM MEEIEM Iflf EI;jtiEl i 


I ■rmo ■tIr!iVJi OiUfi!^ | 




■ ■littM MtfMm gp^ing I 

■ Oryr^MrTgni Oifgiga I 


I EpTiia EilgiTiaM UKEI■tys^MierF;^ ■ 

■ Eir.-f;»;i wnt^mm wrnma ifiiirMH ^<f.-}.-»itrgnEi ■ 

I Eir^^.rti giiRiHa ifcmiiEi iPiiiJii o>»ri^o giXt!«v-a gi«!?ga | 

■ Wtgr;;!! IifJV!;! KlinnCT MSgCIM ga;gT;ia Mr«!IR.^ i< I*!!iI?lMFnicM glgpoi ■ 

■ E »A i *j g miti fcyA^iE i EnzEi MiWM OETTTi Mh g r M m *%uz m m * f M m o « t « t =ho ■ 

I ETOtiK;! iilMsIiM iaRiSin HCTM gfimO ■ 

■ Er»gtKi ■iiRnii’i Bnrg^ «iiM ft/;v>/;i Kicigini me^em oijikvo MrigiM giBinEI I 

■ ■«»[■« giTg^TKl I 

i ji«f i{;t ■ I 


I ■«>#;, !■ MnrigM HA-M ■ «» 

I MrAriiM tjj^W BfelRillll 


I m* t *:*m msvwm ■? !.« MritM 

■ »JHai M-a waiM 


I OtKrVJ MrltM ■«>iOi Mc?W=a 

■ Orfi^if MiV^^ MCM ttgiKM 


I pjflitTO ErW'Z \rM ^E^Oi * MiHim * 


■ ■«]&>?■■ mdiM 


I Eg«!:lH(tl >1 EgT^gm ■|Qpi|OI|E|{I^E|l 


| OMtil3| 




Ratalnina ““ ***““' ^®‘®' ^®‘“* Copodior Avorooo ^®*®* 

Phaaaa ® Poles Roetifiod Duty Roctlfiod Copocltone* Copodtor 7 *^*^ lnducfanco( 

Voltago Cydo Voltago /«n ..m 


I g]ngTaa idr;W.-El | 


lOtM Moauio I UiiMaKAM I lOUU « , _ 

Inductor Dlodo I I Inductor Capacitor Dlodo volumo Volumo 

Vainm. I of I Val«^. Vnlu„. Valnm. ''®'“‘"® ''®'‘'"'* 


Volumo Volumo 
(m*) (m*) 


mj«H:giirf wR!iaFii 


I ETiT:«gl ig^Sl.itl i*»iW-hpff M 0/.-«;Mg *].Oi OTPg^l 

■ EHhM gqggpni E»;g!^ | 

I IeH EtltVhM Er*!CTl UMiTTM |jjtr?[:y«ll 


II 


iBjiiiaai 


i EWOTil lilCTigl Big^TiEl WESEm OTinaa ■ 


IBiit'TigaEgnaOrggiiWiir 


IDE^HI 


I E^-ggM lilESa WlFMm BIEHgi ■ 

■ iarnr^ iraninH | 


I laiiVifJglErltWM igr«:?a<| ^ESEW EiggtEO 1 


■ El<Ti!^ | 


I BP-Z^AJ ■»»c«ci»>j I 

I EftHJH gtlRiga BKgig I 

■ EKfr;;i:i Ei»{!»yAJ WFTTiifni ■ _ 

■ iti:>^?)tM EifcgHi iariT;Hn mseEWI g;»>tVM WKWkmz^^^^ I 


i lj«>t|r3a| _ 


I MW Mil Erfc! [«a RfiWilil ■ 


lEmiPii 


■ WiMry^ EiT*!^ | 

I MrBIcM Krlrlryj glPnia | 
I ■»»![>■ gtr*™ giigiEi I 
1 »rT«MM gTOITni | 

l o«m-o»tiyiMgyi!iia I 
I OtiiP^ J oagiiai I 

I ■iKKO ■«lt»!i^ ■ glp!^ I 
I ■»!<:?!;■ ■«yi!iSBEilP!^ I 
I ■«»r^-0 ■«Ii!«VJ EilgHfl;! | 

■ ■«if«!:a MtiiNs ■ griiTnini | 
I grliMO MiIiMM ggililEl ■ 

■ ■rir^M MifrlrV-M ■tTiT^mi j 
I ■iKkfcM Mifirryji KEiiiii ■ 

■ ■iHtiMMiiyrMiiniiini l 

I ■ilr'rt ■ ■ntt^mi 1 


I g]E:*Iiro#^ ■K-f [MM 


I OrKr^ liM MilvUM !■{»;'!« 




IMOPIMM Mi-M WMl 


■ r . fi Wtw mt bim mmm 

I MiM* WeXSM WKMm mWB 


I MEM Mini M 

■ m-f 


■ m n ^m miWFm WA-MM« l lr;i M MR ' IrtM 
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Appendix P. Converter/Output Filter Mass and Volume Calculations 


fHwV-wiT 1 iVyn i>rT M —_ 


Phasos Poles 

Can 


Output 

Voltage 

Total 

Capacitance 

(mF) 

Total 
Capacitor 
Energy (J) 

Total 
Capacitor 
Mass (kg) 

Average 

Output 

Current 

Total 

Inductance 

(uH) 

Total 
Inductor 
Energy (J) 

Total 
Inductor 
Mass (kg) 

Module 
IQBT Mass 
(kg) 

Module 
Diode 
Mass (kg) 

Number 

of 

Modules 

Total 
Inductor 
Mass (kg) 

Total 
Capacitor 
Mass (kg) 

Total IQBT 
Mass (kg) 

Total 
Diode 
Mass (kg) 



■iiiiif I* Mtnwm ■EWiiil» w:g ^ 

MEglEMi Mt MESatil I 

1 ■aacMI 


mmm i 

■ ill B M B11 MM I Ml 11! I tl ^'^-^'t-M ^ EiqpaTECT I 

■/i I i I . B ^gr ;by^ M .W M i I iyi M M^ignai tlCT ^kl;M I 

K55CT W H 11 WW !B?»M«»«>'i<a^MI3jRiPi^'gl:JMKggl>MgFglM I 
■iiiyzM ■^>jr:ga ■ktofm ■ e i w mcbcm i 



m PJtlcM MiltM ■ ■igPI g gfe’Kai 

g 8MU Mig g g MM 

gki.-f;■ ^ltty« »g?ggii 1«g gTigJHl 

mmm 

gEkfPM 

_l iEMlHi 

■OEHUgtyim 

m kkKtm Miityz ■ ■igai i^nai 

g>3«« ■riIiCTM»->4:M ■«lEI.-»HFggE» 

g kU-fcB ■iliyiM MEaM ■*» kM Mnigyg 



Phases Poles 

Can 



Totat Total Total Avoiego Total Total Total Modulo Modulo INumoer Total Total TotaMQBT 

Capacitenco Capacitor Capacitor Output Irtductanco Inductor Inductor IQBT Masa INodo of Inductor Capacitor ^ ^ DIodi 

<mF) Enorgy(J) Masa (kg) Currant (uH) Enorgy(J) Mass (kg) (kg) Mass (kg) Modules Mass (kg) Mass (kg) ^ Mass (I 



■i»T>Ti» Mcfiisngi MifriHM gtTtVJxMP^gnjg M.I^.I.M gk!;gB MiPUg l 
■kutn^ Mcreag gl tW7:M ■iX>yi:M BEEP^™ g tck]@i I 

MlIiI.!:M griWEM gtliVi^g ggngB MIMM | 

■/»T>yj i^ikM miwfm ■tity#g HBfp™ mFom mmm \ 

B/ 11 IM HI f Hi MM B N Ml IBI il ^l»!»»MgigqHii»gaiM l 
■TmrgiMt^x-EM ■.»:!■ gif ii:igg^i»^«g;gqMgR5iag^gM i 
■ATiU M Hi|' WMijij/ Ml 11! I iiM '' I BBil I i BM \\ f 

MES5M ^EM g«inyAM K^iif gegggg 

B^||||iBBii|! ii MM'II I'MM BI ''Ml i » PX>tagt!igM 

MEEBM »***/i^ ^E5EM Mlti:iir■ K* t^»M 

^aEEMi gifiyAM gdtix-yj MjI m.^ ■ kknj i 


■ gnaiM gtityAM gtie \m m k!:gM 

I gtKtrif g]nig m WPIJI 

I EM3 MEHM HEE9 KSSil 





■ MiirkiM BiMix-ya^'uiaiMi taatM gckkgiM 


■ g{«!tiiB gtityiiM wKJ'Xm ■roiggBrenai 

■ gliMPg grUtyAM ■*gkg g ikiia 

■ g»yj:lg M ■Eligy 


I Cspsctlsnco I Capacitor I Capacitor I Output I Inductanco I Inductor I Inductor I IQBT Masai DIodo 


Enorgy (J}l Mass (kg) I Currant 


I Energy (J)l Mass (kg) 




MTIliBBEBE^Ml _ _ 

■Al.lia McTi gEM M.]^4W BEai ■ 

mmm gKiri.gM ■>w ■ 

■Vi I i l . B «ggM ifrAM g * T .VA« gE k»g m^SOSM ■ 

■AiiinM MEEESM gnmeM M.L klf mti Qj I 

miiitjiM igigiM 1 

■iTlll^■MBEgH _ 

lirAiulrMl BEESSMi ^E^3M W ' 

* n I BM H MMiWcMgtI>yAM^k<:U-yj»igg|gMlggCTiBggaM I 
■AMiB ^EEEMI ■« W giligg gElWg gkkJkM I 

WEES3M giliUgg^XJikyj gEkt-’B ^ky#:M i 

B11 MM Hi MM I i/ Bl 111 Ii M EkJkMgEkJkB»kt;jr^ i 
■AI«l.BMEgM M«>i£M BEESB^aiBEnMKnia»heM I 

»n»Kg MtftCCTB » El*lgB |f-kJtjpg I 



gkfclB giltiKcM gWOEl 

gkklcB g»l>yzg g^XlPM Mtg kg g E! Ht.-B 


g ky#Hi T«ftyai wmssKim ■tgk» gwii;! 
gkWg«»«yAg »?IiI:MM«gkMg|.-EIEB 
gkMI:B ■*l<lrAM gPXJrri gigHg ggRKB 


g kkJkB gtf«v;ai ■Pang wngeg graarai 
gpj Jfg giTtBEg I 

g'Kkijai g«K«yz g I 

gky#:agilPHg|_ _ 

gki;»g giiiy^g mtfiXim ■fgbg guicitB 

gpj^gg giKiykg Mrg kg g E lrJ?B 

gktii;^g*g>yzg »»>:»MM_*gai[ggEig 


Ph...s "•“"^'"8 Pol.. 


I I I I I TTotal I Total I Total Modulo Modulo I Number I Total I Total iToialiaBrl 

V I Capacitanco I Capacitor I Capacitor I Output I lr>ductanco I Inductor I ItKluctor IGBT Mass DIodo I I Inductor I Capacitor | ».«• n(n\ I 


I Energy (J) | Mass (kg) I Currant 


Energy (J) Mass (kg) 


Mass (kg) I Modules [Mass (kg)I Masa (I 



gAi»T>» BEEESBI m*%rkm gii»yAg ■feiiUia MjE kP-m mu K?g kt.TriH i 
■AT»T>BMgT»I;nggf*!Kggigtf»yAgl I 11 >JTnrWgkT:lfg-»>■■■! 


_ ifcL.iM MilfiTrag gJ-XEg MSgaCM glcTJ^ gcTrT.-fjg I 

■AYigai BEEESM BEKM M*n*ykM WJrM^vtm ^EE^B fc^j B' ' 
B M11 B BrmB pi i^Ag g<g<yAg BJjTgg MacrntM gk^CT g irigtg i 
BEDEW ■»nyi?M ^ESBB fc iikiiW i 

■ATtlg MtWEag gifiV/g 1 I I I H glckt ■ g- lil«K ■ I 

IIIK^Sgi^BEBM^«*^^^*^*y^^^^-<^^^L_ __I~~ ~ 

■AIM MPETM mjriM mmm ms^m i 

■mMMcgTMMiiBSMgqriHiggarngBCgMigiciMKECT I 
■MMHl^^W»W:BiB»X«yABEE^M»iK%MBEtMBEIEB I 
gAM. B g t I ryk M ■gamaaigKnEM g EEgg »kgrg 

■gggBMBEMWaJMiXaiaEE^ 



■ gkkfrg gtgtyAg ■«iEkBI gEMiB 

I miiMM rotigig MigM gEaai 

I g k E g.g mmm B ga f M igk-g angaoi 

I ■cklii.g gtgtyjig ■<>-»:■ ■qnig g !!*■ 

■ g tr;g«B g«I«yAg gigEg gkl«KB 

■ gTiIMgii:iCTgig^inMM<lEi.-^gktlEB 

■ g(«I»»B grliyiig ■rgk^l gkklEB 

I mw^m gitrVAg »>>JrMMti!gEg gEfc»:B 
I gnaEB grgtyAg »»>»;■ Mnnsg gEkgg 
I ghblg grgtyjg ■>gkg gkkJg 

IgkEliB 

■ gRggg gifryAg ■OEM gEWHi 

■ gf. T HE B g . i t iig 



_ gkk»g 

■ckyjcBEiliMI »i^gnMg«1gHggtfM 

■yAJIggrlrkkg 
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Totol 

PC,., c^L ^ ‘"c:z 


I ilSS^ HBH EESn Kl^a I 



Number 

of 

Modules 

Total 

Inductor 

Volumo 

(m*) 

Total 

Capacitor 

Volumo 

(m*) 

Total IQBT 
Volumo 
(m*) 

Total 

Oiodo 

Volumo 

<m*) 

Sorvico 

Volumo 

(m*) 

tdtal 

Volumo 

|m*> 


■ ■»XtHi:^Miiilh’;i l 


I >/»:«»!■ ip<ii>j<M miwzm ■»x»;i!»si mrnKm ■ii«SKrHi | 

I MrlgiliM ■gCjM IcTHFa » HrZH gtltl [»! J giXtj U-M | 

1 »zir«i»M ^ji’iiiM gx»!»>ai^.kX:iyj i 

I TigiliM MilBcM latirilBljjt:!?! ■.tyjcfKEItCT mtMUim Ktltlk!;! ■ 

I M il'-M BEfZM ■«X«:V»nigt^4;l.-»J TiirM ■>ItfP!agtXill<!;M I 

■ «»:il«« m*Wim Mi>i:i:i« ■At.-f :■ gETiCTI M ■ililiril gilil k!:l ■ 

i ■ritrtXta ipjiM:i miwpm ■.-g»:« ■ mcpt%m m ■tiii [«!:i KiXiiRr^ | 

I >^i:i»i»iipti«i>t M mfwzm aiii:i;«Ti MtXtW*:;! ■♦riiRii | 

I waaaw MUmm gX»:i!iIB ■4c!*iB KETiU m IFai MiftWM giXtS b!;l I 

I Mir*x«B M B • UTM MgBM ■>][«>!*rB ^kiigg m*hi m;1 ■ 

I BP 111 w. M mtiki ■ gx»;t:»m ■ENga Kiiinra w*h\ k! j i 

■ MI h! i!. M m*wim ■iXi>:iiB tmwzm bixo Ka j ■ 


i ■fXi!»>iB MiiRM ■iiK;M Mrii?« 


■ m.-M Kjiixini MgRiiM ftTM m*Ti^.m 

I Mv^ 


mai « «x»: i™i _ _^ 

I mwim gx»:»;«t ■ MORtM ■»! bep^ 

I gXt>;»i« MilRM ■ilKM m*Tli'tmmX^‘W 

■ ■ IBM BXtXtiM m*pxm ■ix^iM mFmm 




I KiTgiiga mt^m mr^ 

I m Wim gX«TiT«TB MiHiM ■<>>; J MCTn>» 

■ m9fim m»ir« wfx^ 


PH„.. Pol- 


Total 

Vohim. ^ 

* Currant | uH) 


TotiM Modulo Modulo Numba- Total 

Inductor IQBT DIoda . Inducb 

Volumo Vohimo Velum# ®’ Volum 

(m*) (m*) (m*) Modules 


lnS“r C^tor 5S1 

Vdum. vTm. Volum. 

(m*) (m*) (m*) 


l lil" 11 i MIEE31MilBTrTYB KtfJl-.l&M m&X’fiM M MM ‘ W MirnTTiM g»lilhr;« ■ 

I »itrix»B mmm aiXi;*:*?! ■»>i:i.-vj mt ul-i; ■ m w&m ■tXii [«!;■ giXiitc!;! | 

I ■ti:i ■ gntiTCT mpjiVr^ m fom m grirs uim giiiyK!;i | 

I |iPcl«!it?Mi MilM gilTrliVa M.-WiXiM MI<!;1! ■ » WMM KiXtS WrM KtTri WM | 

I uptn^-Mi miwsm gtitririT j m ■ kVKm » tL^m ma*t t<:;i g«T*yR!ii | 

I ■iltHi! JgtTgTtia ■ 

I —i I i.i M MiiM gx«ytm if MR;i ^ ■iftirra gii tSKia | 

I ■rg!iXt» ipcii KTM milfrm ■gri!«iii ■ gtcrnra mwz^ bim km giiii tt!j i 

I EggiM McWlIM m*TlrMM EttMiil mhlhM BpIR3i KiMb!J ■ 

i MPil>J< !iM TWiM ■atr>!*ni TcnOi m ■ gtXri Uim giXrn?;! 1 

I gir«!r»ii ijjtiliM y ^ MR.iM »«xtii*:j gtwgra i 

■ >/»:»Xia Mn»k*M TtigM git«!»r*ni KHiiai gtXiiK j gitijc!^ ■ 

■ !?■ ■«]lr2M gXi!»!«« MliTWm 1Btf il t«;a 1 

■ «»!*>»■ mutimmnitVAm mryji^jm^iiVJkmmpj-jm MriiH*:;! Kiiiik!;i ■ 


■ MiM CG55a wBEm 

i ■ KtXtWii ■ mtit'im 


I m^JM HFiaM 


■ M WitM fciKtIO ^ BERil mfUniM w 

I mfAm gxgtiiiB m*H'Xm ■tX;;^ ■bkim 

■ »»;!<■ ■«x«!i!«l ■ MilRiM ■<!>;;■ Mti:;V Mlfim 

I ■ frZ ■ gXiltltiB ■iygiM MtKJ WiX^ MFItTnB 

■ BXi!i:iTB MilKiM ■«!»;■ WtlM mWim 

■ m^m mm mmm 

■ MIT.-IM BX»r*!«;BMqRI«M M«Xj J ■tXiaiiKIBiM 

I mwim gxi:>:«^B miK^mii^mmm 

I gx»!»:ti« ■tjR^ m»j>:imm*v^mwB’jfiM 


_ . , _ Total . _ , Total Modulo Modulo M.imKor Total Total 

„ Output „ f. Capacitor Total ^*** Inductor IQBT Oiodo Inductor Capacitor 

P”!** Vol^. V«l™. Vi™.. V.I.... Vlum. “* Vl.m. Vlun.. 

|n.F) B-WW) C.™,! uH) B,.™(J, Modute. 


Sorvico Total 

Volumo Volumo. 
(m*) <m*) 


I »i>:»XtM ^Hi3M m*iikm ■*Xt:t;iy »Pcf« ■ giXti wm gxtsitTa | 

■ Ml 1111 ■ ■ilZM gXiyill P>>Xgg ■-•gMiM KCnMIKiM gXtijtM 1 

■ ■zgajMtiHMgiiziiKgtytiiijjEaiMsggx^Mfti^aMEMgqpngagiXiTttTa i 
I ■riiigiM 1^^531 Mrw*x?r«y m’>j:u-a mb m gnCTRjigtititcrii | 
I m • n r M mtwsm KtXt!*:*; ■ ■fcj.-iii ■ me^twm » K:it^ mtftiit'.m gx«it£!-i | 

I »l«!iX»B^JI;l:!:MiMKcMgX«!>:iTBB'>4:kaiMjgkXf Mni!aMfr#^iM*XilL*! JMrSk’J | 

M ^ ■ M WEEHM m ■tx^trim g'>X:i.-»i Kg j«»M WKk W **xit Kr;’! | 

I >z^>x»a ^4.ij'i mtwim g«'x»y«iBB»Jj^iia m-~^^^*x:m ■»xi»«;ii Ktirik!:! | 

i ■r/i:>x»M MiVtii tim MEzwm m w^i»m t«!a gx»i him ■ 

I «CT>» MjlZ»MiX»!i>iM H-4:J.ia|ip.-g[iXtMMgnMm^ ■tX*iI«U gtgB^ 1 

I ■ii^lXi ■ M H I M M W MWi W gEiCT ■iXilRa gXtlfcUl 1 

I »i»:»x«M IP»j 11 mtwim ■»Xi!t:im gggia ■ gixti [i: j TOtiitc!:! a 




■ ■»X«H»!.^MiXitH;li 


i HrXtrtlrf BEygtM ■■X-JMMrXsy^ MFXItiM 

I ■ wim ■tXt.yti ■ MtiRiM m*Kim ■«x:{;m mmm 
I m K;kM K»X»y«IM MiX^M ■«!«J MRIM 


■ ■ i:M ■tXi.ViI ■ Mlt^M ■»X.^.;M 


I ■ ■»Xt!»!trB ■raFM MtlWM ■«X;f?M MRIiRH 


■ ■!«;■ gtXiyjIBMPXg iM Mtiwai »13;a 


Output Co^SStor 

Cap^^c ^padl« ^tput l„du^nco( 

^ (mF) Energy (J) Cuirint uH) 


^.T rZi* Numbar T.™ t.™ toWiobt "»i 

S. v^!^. ot v.l.m. “I* viui. vwvn.. 


Voluraa Volumo Volumo 

(m*) (m*) (m*) 


Volumo Volumo 

Modules 


I MtiwiM KtXi!#!*?! Mcinga | 

i M ■ 


■ MrXri [«!:■ gXiIltfil | 

■ BrXti t«!;l gXitKil I 

■ ■«T»irt!;lB»Xi'Ck?;H 

■ KilrSfiM KiTrlK!;! | 


I I'iXTOTB BtfeXtUrJI ■; kVJ’M 1 

_ _ ■■?XrriTa^>JI;kfrlMa-^!c^M»g8K!B| 

I ■ ill 111 MESgM M>ii;>!im M gxitrriigx*; j<T;i | 

I mK*t*m IP m M mvkim Miyim giXtiria gxtstt;;! ■ 

I ■rZ»!iX*M Ml H 3M TiNMl ggTtgM ■KM J ■ 

I ■r;»!tx«a bee^M ■ii'M gxr*r>y b^^xhi ■iXffTriagiT»ifcr;i i 

I >i»itXta WEEUi M«Xrai gx«y*iB ■.-jcvjna ptc?iXiB ■ofzm ■■Xriya gXtih! ^ | 

I ■rgaj McUrM MiT*y*'Sm grXrlKJ Kilrlk!.-! I 

i >i*!«Xia WElEMi g1K!i!CT WkkI:B M f fdBI gXilTT;! ■iltlfcT.-l | 

i >ii>XiM MEEEai IrfrAM gxroni ■.^XW:M t spa 1 frfM »*X?n»a gXrXk! J I 

I ^KfXil ME 3 aWW*^'M WgXiTT^ | 

■ ■r7«TiXiB iptip^aa m!Wim gX'tV«M»>-I:).ia gRt^ ■lXl»l^^^^Xl^^f;l | 


I ■ WSm gX«!iriT MijRM MrMM ■»!;; M 


■ m wA U gxty^ ■ iwa n ■ « ? « mnmm 

i MiXiIiNni ■roRM M«X>T;M MtXiM HEgTWi 

i mtiigjia mMm mga Kftiaa 


I M iTiiil fciiKHiM BEIRiM 

I giXi^IitaMEFiTiM 

I gXiVti ■ mtpjm MiXi^^ wkpx^ 

I ■R^Ji gE55IMM»K«JI 


IKi^J 


Ph..«i "•“i’'"" Pol., jss; ‘=^““ “to." 

vottaga voiumo 


I gtXil t«!.^ IPXffg^l ■ 


I ■KM gT«yCT M.)K« M ■ ■ X .-yM BUgJPmH 

I mwitm gx«!iM ■ m*Pim ■■x.-rit hfim 

■ KE55I1M«K<»M>K 


Total Modulo Modulo ajurnber 

Inductor IQBT Oiodo . Inductor Capacitor 

Volumo Volumo Volumo Volumo V^umo 

(m*) (m*) (m*) Modules 


I ■nrgj| M.x»yin jijj Pg : .M ■jxgangxMiiii^^ 


Sorvico TotM 

Volumo Volumo 
(m*> (m*) 


I >^!«X»B MEEI3W Pf>Ki1 gX»!«MB M friiBi I 

I »ii>x»l WEIEMI M!Wim grXtTrTiIB ■.-ftlsXiM ■iXtirirJBXtStcriM | 

I »;i:»X«M1JJ-l.MTiXMg«X>yiniMa:].^l »t^^g ■iXi»>!JgiXiifc!;l I 

■ «!!;?«■ M JillM MrlLirM gRiMl ■«X«iriM gilfiT jtM ■ 

I Kgx;a M«Xt!:<:M TOEM gT^yril »cT:I>B MBcM ^«X«^F«^^ ■ff.IkTJ ■ 

■ MrANXia BEB3Bi MtXiiM gxrrCT ■>>J:i.-vi f 11 gw»*B Ml j B ■«t«tt^ gt«xi<r.-i | 

I >i»>XiB MEEEM »»X:feW EgrgW »laitM 1 ri'ii gXrirc!^ I 

I >i»>X>M |P< l«!.-M MrlEM KrXOirrni ■ WcITCT U*h\ lt!;l gXtIK;l | 

■ ■rit!«Xia ^lEi »!»;:■ ffxranig>>X:U-»i ■HMgwitgu m rtJfnBTCTgrtmc?!! i 

I >i»:«XiB W **Tf»Trr=wg*Xiti<?;i ■ 

I ■ i I 111M MEESM gfiyCT l:i g-’M ■njjf :■ 1 P-y Brlit t«!?l BtXiIk?;! j 

I ■»Xiy»TBB>>j^ai KHM Mfr>-»gXiY[>a gX»IH.^ ■ 

I ■rZgiXrM iptfc ^4<^^ gTiyrm ■.■frrX;M tc!tx.-a m km grTri uim k«Xi; uim i 

I ifr^«!iX»a BEBUM M«X:?*M gliyCT »l; t Ui KVJtlM WflM ■tXffTta gx«l k!;l 1 

■ mri*:*hm Egyira ■>>i:uvj^ i:iB»M»paCT BiXtnra ■ 


■ grXtlrMM MPIRM Miyjai MrX.-^ MFIEM 

I ■ gf ly ri ■ ■«xtMi rntPim m«x.^ ■ bexi!m 




I mw ii mKfiTimm iPxm 

I M H tJl M WEEM Wif:^ M BBd(;TMI 

I mwsm Bx«:i!i7B MgRiM m*'9rhm ■iX^iM mrim 


mimESSSl 




i MOl KESiTB MHFM M .Xj;M M « X ;iM »gai 

i MMm BEsga iiBBaB 


J iKaaM ^ 


I M Bij ^ HfiTiTiTiI M Mt]K^ 

- 
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